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SUMMARY 

Full  details  are  given  of  the  untruncated  second-order  orbital  theory,  of 
which  a  r6sum6  (including  the  formulae  finally  derived)  was  published  in  1983. 
The  analysis  only  takes  account  of  the  first  two  zonal  harmonics,  J2  and  J3 
(geopotential  assumed) ,  but  the  extension  to  an  arbitrary  zonal  harmonic  will  be 
covered  in  a  later  Report  (Part  2) . 

The  principal  merit  of  the  theory  derives  from  the  extreme  compactness  of 
the  results  obtained  for  short-period  perturbations  in  position  and  velocity, 
these  perturbations  being  expressed  relative  to  a  system  of  spherical-polar 
coordinates  based  on  a  rotating  'mean  orbital  plane'  The  use  of  osculating 
elements  is  thereby  avoided,  and  mean  elements  are  defined  so  as  to  give  the 
coordinate  perturbations  their  simplest  possible  form.  To  make  the  theory  both 
complete  and  compact,  an  intermediate  set  of  elements,  described  as  'semi-mean', 
is  also  required. 

A  numerical  assessment  of  the  theory  will  be  included  in  another  Report 
(Part  3)  . 
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1  ZHTRODUCTXOM 

The  'main  problem'  in  modelling  the  motion  of  an  artificial  Earth  satellite 
has  always  been  recognized  as  the  representation  of  the  perturbations  due  to  he 
low-degree  zonal  harmonics  of  the  geopotential.  A  considerable  literature  has 
resulted,  much  of  which  was  surveyed  by  the  present  author  in  earlier  papets^"^. 
The  last  two  of  these  papers  gave  a  second-order  theory*  for  orbits  of  low  eccen¬ 
tricity:  Ref  2  was  confined  to  the  perturbations  due  to  J2  ,  but  Ref  3  covered 
the  effects  (for  low  e  )  of  the  general  (tesseral)  harmonic,  ;  as  in 

Ref  1,  the  formulae  associated  with  were  shown  to  be  applicable  to  luni- 

solar  perturbations  when  f  and  m  are  taken  as  negative  and  zero  respectively. 
The  novelty  in  Refs  2  and  3  lay  in  the  compactness  that  could  be  achieved  by 
expressing  the  short-period  perturbations  in  terms  ot  a  particular  system  of 
cylindrical-polar  coordinates;  however,  when  the  theory  was  extended  to  an 
arbitrary  elliptic  orbit  (with  formulae  therefore  untruncated  in  e  > ,  the 
corresponding  system  of  spherical-polar  coordinates  was  found  to  be  more 
appropriate. 

The  outcome  of  the  extension  just  referred  to  has  been  a  complete  second- 
order  theory  for  Jj  and  J3  ,  and  a  summarized  account  of  the  new  theory  was 
presented  to  the  1982  lAF  Congress^.  The  present  Report  gives  the  full  details 
of  the  purely  analytical  part  of  the  theory,  including  (in  particular)  the 
lengthy  expressions  for  the  perturbations  in  the  osculating  elements  as  well  as 
the  strikingly  compact  expressions  for  perturbations  in  the  coordinates.  The 
original  intention  was  to  include  details  of  the  hybrid  (semi-analytical)  compon¬ 
ent®  of  the  orbital  model,  the  component  that  models  the  long-term  evolution  of 
the  chosen  set  of  mean  orbital  elements.  However,  the  resulting  paper  would  have 
been  unduly  long,  so  most  of  this  material  has  been  held  over  for  a  later  Report. 
Further,  formulae  have  recently  been  derived  that  generalize  the  J3-results  (and 
the  first-order  J2-results)  of  the  present  paper  to  an  arbitrary  zonal  harmonic 
(Jf);  these  formulae  are  unexpectedly  compact  and  merit  priority  in  publication. 

Thus  this  Report  constitutes  Part  1  of  an  intended  trilogy.  Part  2 
(Technical  Report  89022)  will  cover  the  purely  analytical  theory  of  the  general¬ 
ization  to  ;  since  each  ,  for  t  >  2  ,  only  needs  be  covered  to  first 
order,  the  results  for  different  t  may  be  superposed  linearly  to  provide,  in 
combination  with  the  J2^-re3ulta  from  Part  1,  a  complete  second-order  theory  for 


*  Since  all  other  geopotential  coefficients  are  of  order  J2^  ,  a  'second-order 
theory'  is  second-order  in  J2  ,  but  only  has  to  be  first-order  in  any  other 
covered. 
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the  zonal  haztnonlcs.  Finally,  Part  3  triil  cover  the  evolution  of  mean  elementa, 
the  topic  dropped  from  Part  1;  it  is  a  consequence  of  the  revision  in  intention 
for  the  present  paper  that  a  large  number  of  references  to  'Part  3'  will  be  found. 

Section  2  of  the  Report  provides  sufficient  background  to  make  the  paper 
essentially  Independent  of  the  author's  previous  contributions^'®  to  the  subject. 
Some  of  the  more  inportant  notation  is  introduced  in  section  2.1,  but  in  view  of 
the  complexity  of  the  notation  used  throughout  the  paper,  a  Li st  of  Symbols  is 
included  at  the  end  of  the  Report. 

Section  3  is  devoted  to  a  number  of  important  issues,  and  in  particular  to 
introducing  the  notion  of  'semi-mean  elements'.  Unli..e  mean  elements,  these  are 
not  free  of  short-period  variation;  their  utility  arises  from  the  fact  that  the 
independent  variable  of  Lagrange's  planetary  equations  has  to  be  transformed  from 
time  to  true  anomaly,  to  permit  an  untruncated  integration,  sc  each  semi-mean 
element  differs  from  the  corresponding  mean  element  by  a  quantity  proportional  to 
the  difference  between  true  and  mean  anomaly.  Considerable  attention  is  devoted 
to  the  way  in  which  position  and  velocity  are  derived  (via  semi-mean  elements) 
from  mean  elements,  since  the  whole  ethos  of  the  paper  is  to  stress  the  import¬ 
ance  of  an  algorithm  in  which  (short-period)  perturbations  are  applied  directly 
to  a  conceptual  'mean'  position  and  velocity.  It  is  essential,  of  course,  that 
no  accuracy  should  be  lost  when  the  algorithm  is  inverted,  so  that  mean  elements 
can  be  obtained  from  position  and  velocity.  This  requirement  is  of  particular 
inportance  at  epoch,  and  is  easier  to  implement  than  is  sometimes  supposed;  a 
method  of  inverting  the  algorithm,  iteratively,  was  described  in  Ref  5.  Two 
aspects  of  orbital-element  singularity  are  discussed  to  complete  section  3: 
first,  a  difficulty  that  arises  in  velocity  computation  but  not  in  position 
computation;  and  secondly,  difficulties  in  the  long-term  propagation  of  mean 
elements.  The  latter  discussion  is  in  the  nature  of  an  introduction  to  Part  3. 

Specific  expressions  for  the  perturbations  due  to  J2  and  J3  ,  constitu¬ 
ting  the  meat  of  the  Report,  are  given  in  sections  4  to  7.  Section  4  covers  the 
first-order  perturbations  due  to  J2  ,  in  both  (osculating)  elements  and  coordin¬ 
ates.  Section  5  derives  the  second-order  perturbations,  due  to  J2  alone,  in 
the  elements,  and  section  6  derives  the  resulting  perturbations  in  ccordinates. 
Finally,  section  7  covers  (for  both  e.’-ements  and  coordinates)  the  first-order 
perturbations  due  to  J3  alone;  as  already  remarked,  these  are  regarded  as 
contributing  to  ,  he  second-order  solution  of  the  pro):i.em  with  both  C2  and  J3 
present. 


TR  8BO68 


7 


At  the  end  of  the  concluding  section,  section  8/  the  most  iirportant 
formulae  derived  in  the  Report  are  referenced  for  convenience. 

2  BACKGROUND 

2 . 1  Osculating  •leoisnts  and  related  quantities 

Our  starting  point  is  a  standard  set  of  osculating  elements  for  an 
elliptical  orbit,  viz  a  (semi -major  axis),  e  (eccentricity),  i  (inclination), 
n  (right  ascension  of  the  ascending  node),  (0  (argument  of  perigee)  and  M  (mean 
anomaly) ,  all  being  functions  of  t  (time) .  The  mean  motion,  n  ,  is  directly 
related  to  semi -major  axis  by  Kepler's  third  law, 

n^  a^  =  H  ,  (1) 

where  )l  is  the  (Earth's)  gravitational  constant. 

We  also  recfuire,  as  an  intermediary  for  perturbations  in  the  'fast' 
element,  M  ,  the  quantity  a  ,  )cnown  as  the  modified  mean  anomaly  at  epoch 
(where  t  -  0);  this  is  defined  such  that  (with  t  a  dummy  variable  for  t) 


t 

M  =  O  +  /  n  dT  , 
0 


(2) 


and  we  rewrite  (2)  simply  as 


M  =  a  +  i  ,  <3) 

since  this  shorthand  use  of  J  will  prove  very  convenient.  Despite  its  .lame  and 
definition,  O  is  a  current  variable,  determined  from  M  ,  in  principle,  by 
ta)cing  the  integral  of  n  bac)cward3  from  t  to  epoch.  It  is  clearly  less 
accessible  than  the  'unmodified  mean  anomaly  at  epoch',  which  is  determined  (also 
as  a  current  variable)  just  by  M  -  nt  ,  but  a  has  the  advantage  that  we  have, 
direct  from  (2) , 


M  -  ct  +  n  .  (4) 

An  important  quantity  is  the  true  anomaly,  v  ,  which  is  required  both  in 
its  own  right  and  (in  a  particular  'mean'  form)  as  an  alternative  independent 
variable  to  t  .  The  derivation  of  v  from  M  is  given  in  section  3.3,  and  two 
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other  important  quantities  may  be  derived  at  once  when  v  is  known:  first,  r  , 
the  satellite's  geocentric  distance,  is  given  by 

p/r  -  1  +  e  cos  V  , 

where  p  is  the  serai-latus  rectum  or  'parameter*  of  the  ellipse,  defined  by 

p  -  a(l  -  e^)  ;  (6) 

second,  u  ,  the  argument  of  latitude,  is  given  (and  defined)  by 

u  -  <0  +  V  .  (7) 

We  shall  so  often  require  to  use  the  difference  between  the  mean  and  true 
anomalies  as  a  distinctive  quantity  that  it  is  useful  to  have  a  shorthand  for  it;, 
the  quantity  is  traditionally  known  as  'the  equation  of  the  centre'  and  we  denote 
it  here  by  m  ,  so  that 


m  “  V  -  M  . 


(8) 


The  elements  Q,  co  and  M  ,  together  with  some  of  the  other  quantities, 
suffer  certain  well-known  indeterminacies,  or  'singularities',  for  orbits  that 
are  close  to  being  circular,  equatorial  or  both.  Dealing  properly  with  these 
singularities  has  been  a  major  consideration  of  the  study  (see  section  3.5  and 
Part  3),  but  for  the  moment  we  simply  introduce  some  further  notation.  As  in 
Refs  2  and  3,  we  use  4  af’d  T)  for  the  quantities  defined  by 


and 


5  -  e  cos  0) 


O) 


T)  -  e  sin  0)  , 


(10) 


these  being  well-defined  for  a  near-circular  orbit  that  is  not  also  near- 
equatorial.  However,  §  and  T)  are  not  as  significant  now  as  in  the  earlier 
papers,  which  were  largely  written  on  the  hypothesis  of  low  eccentricity:  it  was 
as  a  corollary  of  the  low-e  assumption  that  U  ,  defined  as  M  +  CO  ,  was  also 
appropriately  used  ss  a  fundamental  element,  with  t]  and  U  replacing  e,  (o 
and  M  ,  but  the  use  of  U  is  not  nearly  so  effective  when  the  analysis  is  not 
truncated  in  respect  of  e  .  In  considering  near-equatorial  orbits,  it  will 
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sometimes  be  convenient,  in  the  same  way,  to  use  the  triple  (§,  JJ,  0  to  refer  to 
a  set  of  quantities  (sin  i  sin  Cl,  -  sin  i  cos  Si,  cos  i)  that  are  always  well 

defined  (and  in  fact  constitute  the  direction  cosines  of  the  orbital  momentum 

vector  in  the  assumed  axis  system) .  Some  other  quantities  -  y,  p  and  L  ~  are 

conceptually  useful  in  the  context  of  the  circular  and  equatorial  singularities, 
but  they  are  only  fully  defined  at  the  differential  level  and  their  formal  intro¬ 
duction  is  postponed  until  section  2.2. 

As  an  aid  to  the  concise  presentation  of  formulae, 

define  seven  quantities,  as  follows: 

it  is  convenient  now  to 

c 

“  COS  i  , 

(11) 

3 

*  sin  i  , 

(12) 

f 

) 

“  3*  , 

(13) 

g 

“  ^  T  ^  ' 

4 

(14) 

h 

1 

fl 

(15) 

q 

(16) 

(so  that  p  -  aq^)  and 

W 

"  q’^  (p/r)^  . 

(17) 

It  is  worth  remarking  that  e  and  q  have  the  same  relation  to  the  so-called 
'angle  of  eccentricity',  ^  ,  as  s  and  c  have  to  the  inclination,  i  ,  so 
that  ^  (defined  as  sin"^  e)  might  have  been  taken  as  a  basic  orbital  element 
instead  of  e  .  Though  there  is  indeed  an  important  parallel  between,  in  par¬ 
ticular,  e  and  s  (of  which  there  will  be  a  good  example  in  Part  3),  ’■“ere  is 
an  essential  distinction.  Thus,  e  can  exceed  unity  (hyperbolic  orbits)  and  so 
remains  a  satisfactory  parameter  when  ^  is  no  longer  real  (universal  para¬ 
meters,  valid  for  all  types  of  orbit,  are  considered  in  Ref  7) .  It  is  quite 
different  for  s  ,  since  as  i  increases  from  hn  to  n  (retrograde  orbits) 
s  returns  to  zero,  so  that  it  is  an  inherently  ambiguous  quantity;  further,  s 
fails  to  distinguish  satisfactorily  between  distinct  near-polar  orbits,  since  its 
derivative  with  respect  to  i  is  c  ,  which  is  then  close  to  zero. 

The  goal  of  concise  presentation  is  further  facilitated  by  the  introduc¬ 
tion  of  the  eight  families  of  quantities  defined,  in  principle,  by: 
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“  003(jv 

+ 

U)  , 

Sj  " 

sinjjv  + 

03)  , 

(18) 

Cj 

-  cos ( jv 

+ 

2m)  , 

Sj  - 

sin(jv  + 

2(0)  , 

(19) 

Yj 

-  cos(3v 

+ 

30))  , 

®j  “ 

sin(jv 

3(0)  , 

(20) 

rj 

”  cos(jv 

+ 

403)  , 

h  " 

ain(jv  + 

40)  , 

(21) 

where  j  is  any  (positive,  zero  or  negative)  integer.  It  is  to  be  noted, 
however,  that  the  arguments  of  the  cj  etc  will  normally  not  be  the  osculating 
V  and  0)  themselves,  but  the  corresponding  'semi -mean'  v  and  5  ,  to  be 
introduced  in  section  3.2.  Thus  the  notation  Cj  etc  world  normally  be 
appropriate,  but  for  simp.icity  the  tildes  will  be  omitted  unless  (as  in  parts  of 
section  5)  the  distinction  becomes  important. 

Final.'.y,  it  is  convenient  to  introduce  here  a  pai  r  of  notation  conventions 
that  will  be  widely  used  in  the  long  algebraic  expressions  that  are  required 
en  route  to  tlie  compact  expressions  of  the  final  solution.  The  conventions 
involve  the  use  of  round  and  square  brackets,  respectively,  to  abbreviate  poly- 
nominals  in  f  and  e^  .  Thus  (31,32,  js)  stands  (in  the  appropriate  context)  for 
jl  jz  ^  +  33  Cjl'jzl  stands  for  ji  +  •'2  ;  and,  sines  the  conventions 

will  be  combined,  [(ji,j2),  j3(j4»js>l  stands  for  3i  +  jz  f  +  j3  «'^(34  +  js  f )  • 

The  j's  will  always  be  integers  such  that,  in  the  last  example,  ji  and  32 
are  co-prime,  and  likewise  and  js  , 


2.2  Lagrange's  planetary  equations 


Rates  of  change  of  osculating  elements  may  be  expressed  exactly  by 
Lagrange's  planetary  equations  developed  in  Refs  8  and  9,  for  example,  or  (with 
more  detail;  in  other  text-books^®'^^;  these  assume  the  existence  of  a  disturbing 
function,  U  ,  that  expresses  the  perturbation  potential.  The  arguments  of  7 
are  taken  to  be  a,  e,  i,  fl,  o)  and  M  (with  t  also  in  principle,  but  not  in 
the  present  study) ;  then  the  six  planetary  equations  are 


_2_  du 


1 

2 

na^e 


1  /  du 

na^qs  ' 


dv\ 

mf  ’ 


(22) 

(23) 

(23) 
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and 


a 

(0 


a 


1  dU 
na^qs 


9e  qs  3i^ 


-L. 

n,2  \  e 


[2  ^ 
de 


+  2a 


3u\ 
aa  /  • 


(25) 

(26) 


(27) 


From  '27)  we  at  once  have  M  ,  in  view  of  (4). 

Before  we  proceed  further,  it  is  worth  malcing  explicit  the  conventional 
assumption,  which  holds  through  the  paper,  that  the  dot  notation  refers  to  total 
differentiation  with  respect  to  time,  d/dt  ,  rather  than  to  partial  differenti¬ 
ation,  ()/3t  .  hs  the  3/3t  notation  refers  to  differentiation  whilst  all  the 
elements  are  held  constant,  ie  in  relation  to  the  osculating  orbit,  3a/3t  etc 
are  trivially  zero,  and  a  misunderstanding  could  hardly  arise.  For  other  quanti¬ 
ties,  however,  the  distinction  can  be  important:  it  is  inherent  in  the  concept 
of  osculating  elements  that  r  -  3r/3t  ,  since  d/dt  =  3/3t  when  applied  to 
(position)  coord-nates  relative  to  fixed  axes,  but  in  general  the  two  derivatives 
are  not  equal;  thus,  v  *  3v/3t  ,  for  example,  as  appears  later  in  this  section. 

It  will  be  observed,  in  (22) -(27),  that  it  is  only  a  and  (l  that  have 
been  expressed  in  terms  of  individual  partial  derivatives,  but  the  other 
equations  can  be  combined  in  such  a  way  as  to  isolate  the  remaining  derivatives, 
and  these  combinations  can  be  very  useful.  Thus  from  (6)  we  have 


p  =  q^  a  -  2ae  e  , 


(28) 


after  which  (22)  and  (23)  yield 


2q  da 
na  3o) 


(29) 


Again,  combination  of  (24)  and  (29)  gives 

d(pc^)  _  2qc  da 
dt  na  3ft 


(30) 


and  this  is  an  important  result,  since,  for  a  disturbing  potential  that  is 
symmetric  with  respect  to  the  Earth's  polar  axis,  U  is  independent  of  ft  and 
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hence  pc^  is  constant  (conservation  of  yfip  c  ,  the  polar  con^onent  of  the 
angular  momentum) . 

Isolation  of  the  two  remaining  partial  derivatives  is  possible  if  we 
introduce,  as  in  previous  wor)c^'^,  quantities  and  p  defined  by 


and 


We  then  have 


and 


2  ^ 
na  da 


(33) 


(34) 


Though  y  and  p  are  not  the  rates  of  change  of  quantities,  y  and  p  ,  of  the 

familiar  typo,  since  c  and  q  in  (31)  and  (32)  are  not  constant,  first-order 
analysis  can  be  greatly  facilitated^' by  derivation  of  perturbations  6y  and  6p 
from  integration  of  (33)  and  (34),  the  assumption  being  that  the  variation  of  c 
and  q  can  be  neglected  over  the  period  involved.  (The  informal  use  of  the  8 
notation  here  is  without  prejudice  to  the  formal  meaning  to  be  attached  from 
section  3.1  onward.)  Clearly,  6y  then  remains  well  defined  as  i  approaches 
zero  or  k  ,  so  long  as  e  does  not  at  the  same  time  approach  zero,  and  this  is 
not  in  general  true  for  6<o  and  5(1  individually;  again,  6p  is  Invariably 
well  defined,  whereas  6c  is  in  general  not  so  as  e  aporoaches  zero. 

We  can  regard  p  as  the  non-singular  quantity  corresponding  to  d  ,  and 
Ref  3  also  introduces  L  ,  the  non-singular  quantity  corresponding  to  M  ;  thus 


L  M  +  qy  ,  (35)  j 

and  again  L  and  5li  are  useful,  though  there  is  no  familiar  integrated  'fast 
element'  t  .  It  would  be  wrong  to  give  the  impression  that  the  integrated 
quantities  y,  p  and  L  are  so  artificial  as  to  be  entirely  devoid  of  meaning, 
however;  it  is  not  hard  to  see,  in  particular,  that  y  is  the  intrinsic  rate  of 
rotation  of  the  perigee  direction  within  the  evolving  orbital  plane,  so  that  y 


i 
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exists  39  the  accumulated  angle  (with  an  arbitrary  origin)  swept  out  by  this 
direction.  It  is  to  be  noted  that,  from  (35),  using  (4)  and  (32),  we  also  have 


L  -  p  +  n  . 

In  regard  to  the  potential  associated  with  zonal  harmonics  only,  the  right- 
hand  sides  of  the  planetary  equations  can  be  expressed  as  closed  functions  of 
true  anomaly  but  not  mean  anomaly.  It  is  therefore  of  fundamental  iirportance  to 
be  able  to  transform  the  independent  variable  from  t  ,  which  is  effectively 
equivalent  to  M  ,  to  v  ,  where  v  is  a  function  just  of  e  and  M  .  We  have, 
in  fact. 


(37) 


where  the  partial  derivatives,  denoted  by  v^  and  Ve  ,  are  given  in  (42)  at  the 
end  of  this  section;  also  M  is  given  by  (4).  We  can  now  appreciate  the 
difference  between  the  total  time  derivative,  v  ,  and  the  partial  derivative, 
dv/3t  ,  referred  to  earlier;  for  the  latter  we  obtain,  on  replacing  M  and  e  , 
in  (37),  by  n  and  zero  respectively,  and  setting  vh  from  (42), 

3v 

-  nW  .  (38) 

at 

Transformation  by  (38)  permits  a  complete  first-order  analysis  for  any 
zonal  harmonic;  on  proceeding  to  second  order,  however,  we  must  allow  for  the 
e  term  in  (37),  and  also  for  the  d  term  that  is  implied  by  the  appearance  of 
M  ,  both  these  terms  being  of  first  order.  More  aw)iwardly,  the  first-order 
perturbation  in  (osculating)  v  ,  which  is  our  assumed  independent  variable, 
contains  short-period  terms  that  are  therefore  functionally  dependent  on  v 
Itself,  and  this  implies  something  of  a  vicious  circle;  further,  some  of  these 
terms  -  see  (182)  -  contain  the  e"^  singularity  factor.  In  second-order 
analysis  for  J2  ,  all  the  difficulties  are  overcome  if,  anticipating  the 
introduction  of  mean  elements  (and  semi-mean  elements)  in  section  3,  we  transform 
(from  t)  to  V  (or  rather  v)  instead  of  v  .  We  require  a  new  relation  to 
replace  (37)  and  (38),  and  this  turns  out  to  be 

V  -  nW  t  0{J^,  J3)  ,  (39) 
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since  n  can  be  defined  to  include  d  >  Mhilst  5  is  o(j^)  .  (The  second- 

order  term  in  (39)  can  be  ignored,  fox  change-of-variable  purposes,  since  the 
equation  is  operating  on  the  planetary  equations,  which  are  already  of  at  least 
first  order;  when  v  is  required  in  its  own  right,  however,  i  asust  be  allowed 
for  -  cf  (99) )  . 

Application  of  (38),  or  (39),  to  (22)-(27)  gives  (for  any  zonal  harmonic) 
first-order  perturbations  in  a,  e,  i,  fi,  O  and  O  ,  as  functions  of  v  .  For 
a  con^lete  first-order  solution,  we  still  need  the  perturbation  in  H  ,  however, 
and  this  is  coiqillcated  by  the  need,  through  (3),  for  the  perturbation  in  J  ,  to 
add  to  the  perturbation  in  a  .  Talcing  'perturbation*  and  'integral'  to  commute, 

t 

we  thus  require  /  Sn  dt  ,  where  Sn  is  available  at  once  from  Sa  (cf  (159)) 

0 

and  hence  is  a  function  of  v  .  On  applying  (39)  again,  it  is  immediate  that 

d(5/)  Sn 


where  the  right-hand  side  is  effectively  just  a  function  of  0  ,  since  the  barrad 
version  of  (5)  can  be  used  to  eliminate  the  in^licit  presence  of  f  in  W  . 

t 

Thus  the  derivation  of  |  8n  dx  involves  just  the  v-integration  of  Sn/fSW  . 

0 

Evaluation  of  all  the  partial  derivatives,  for  substituting  in  the  planet¬ 
ary  equations,  is  straightforward  if  appeal  is  made  to  the  following  particular 
derivatives  as  necessary; 


3  cos  V  , 


ae  sin  v 


sin  V  (2  +  e  cos  v) 


2 . 3  Assumed  potential 

As  we  are  concerned  only  with  the  Jj  and  J3  terms  of  the  geopotential, 
we  take  as  the  disturbing  function,  U  ,  the  combination  of  U2  and  U3  defined 


r  (7)  ^2 
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bnd 


(44) 


where  9  is  geocentric  latitude,  R  is  the  Earth's  equatorial  radius,  and  P2 
and  P3  refer  to  the  usual  Legendre  polyncmials. 

The  Legendre  argiinent,  sin  P  ,  can  be  eliminated  from  (43)  and  (44)  by  use 
of  the  relation  (which  underlies  the  name  'argument  of  latitude'  for  u  ) 

sin  P  ••  s  sin  u  ,  (45) 


and  r  can  be  eliminated  by  appeal  to  (5) .  The  analysis  can  be  significantly 
condensed  if  we  replace  Jj  and  <13  by  quantities  conveniently  denoted  by  K 
and  H  ,  respectively,  and  defined  by 


K 


A  J 
2  2 


(46) 


and 


(47) 


the  resultant  expressions  for  C'2  and  Uj  are  then 


and 


■  f 

^(t)  . 


(48) 


(49) 


the  notation  introduced  in  (18) -(20)  being  freely  employed. 

For  many  purposes,  K  and  H  can  be  thought  of  as  constants.  It  is  some- 
'  times  vital  to  recognize  that  they  are  functionally  dependent  on  p  ,  however,  so 

that,  in  particular,  they  have  non-zero  partial  derivatives  with  respect  to  a 
<  and  e  ;  thus  the  notation  X  and  H  is  appropriate  when  p  is  replaced  by 

its  mean  value,  p  . 

It  has  already  been  remar)(ed  that  for  an  axi-symmetric  disturbing  poten- 
,  tial,  le  for  all  the  Earth's  zonal  harmonics,  not  just  J2  and  J3  ,  the  polar 

component  of  the  satellite's  angular  momentum  furnishe.!  a  useful  constant  of  the 

i 

I 

i 
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motion,  but  a  etore  important  constant  is  provided  by  the  energy  integral.  The 
kinetic  energy  per  unit  mass  of  the  satellite  is  given  by 


(50) 


and  the  corresponding  potential  energy  by 


P.E.  -  "  (f'  ' 


(51) 


hence  H/2a  +  0  is  an  absolute  constant.  It  follows  that  if  a'  is  defined  by 


a’ 


(52) 


then  a'  is  a  totally  invariant  quantity  to  which  the  osculating  semi-major  axis 
approximates.  As  such,  it  is  the  obvious  candidate®  for  the  mean  a  that  we 
shall  need  to  adopt,  and  it  is  very  gratifying  that  this  choice  (a  •=  a')  is 
demanded  by  other  considerations. 

Since  a'  is  an  absolute  constant,  so  is  n'  ,  defined  by 


n'^a'^  «  , 


(53) 


and  we  will  find  this  of  great  importance®  in  the  analysis  of  M  ,  since  it 
provides  a  firm  point  of  reference  for  the  perturbation  in  mean  motion  -  the 
integration  of  6n  ha.s  already  been  discussed  in  the  derivation  of  (40)  . 

3  MCAN  ILBMSNTS  AND  COORDINATES 

3 . 1  Kean  elements 

Let  5  symbolize  the  generic  osculating  element,  standing  (in  particular) 
for  a,  e,  i,  (O,  <S  or  M  .  The  variation  of  M  is  given  by  (4),  but  the 
variation  of  the  other  d  is  (we  are  assuming)  entirely  due  to  and  J3  . 

Each  variation  is  made  up  of  a  series  of  short-period  terms  that  are  essentially 
trigonometric  functions  of  v  ,  together  with  a  long-term  effect  that  is 
independent  of  v  (mod  2jt)  and  the  orbital  period.  The  combination  of  all  the 
short-period  terms  in  ^  constitutes  the  short-period  perturbation,  6^  , 


i 
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removsl  of  which  from  C  gives  a  atean  element,  C  ,  which  for  most  purposes 
(because  its  variation  can  be  plotted  easily  and  accurately  over  long  periods  of 
time)  is  more  useful  than  the  osculating  element;  thus 

C  -  C  +  •  <54) 


The  variation  of  C  with  time  is  given,  of  course,  by 


t 


dt 


(55) 


and  C  (i  in  general  have  two  coaponents;  one,  a  purely  secular  component 

that  only  arises  with  some  of  the  elements;  the  other,  a  long-period  component 
induced  by  the  secular  perturbation  in  <o  ,  This  secular  perturbation  being  a 
first-order  effect  of  J2  (and  any  other  'even'  zonal  harmonic  that  might  be 
considered) ,  long-period  variation  does  not  arise  until  the  second  order  in  pure 
J2  analysis,  and  is  a  cross-coupling  effect  in  J2/J3  analysis;  we  shall  find 
that  there  are  second-order  long-period  perturbations  for  every  element  but  semi- 

A 

major  axis.  We  summarize  the  division  of  ([  into  its  components  by  writing 


c 


(56) 


where  Csec  can  only  have  an  unperturbed  component  wlien  C  is  M  • 

Now  (Jaec  i*  constant,  by  definition  of  'secular',  so  its  integral,  from 

epoch  to  current  time,  is  simply  Csec  t  ;  also,  the  integral  of  Cfp  rosy 
conveniently  be  denoted  by  AC  (Ref  4  used  SCjp)  .  Thus  the  variation  of  C 
from  epoch  is  given  by 

?  “  ^0  +  Lc  t 

There  is  an  inportant  distinction  in  form  between  short -period  perturba¬ 
tions  and  long-period  perturbations,  as  defined  here,  and  it  can  be  seen  at  once 
from  equations  (54)  and  (57) .  Equation  (54)  indicates  that  SC  is  defined  quite 
independently  of  epoch,  and  cannot  in  general  be  zero  at  epoch  unless  special 
epochs  (such  as  ascending  nodes)  are  chosen  to  suit  the  particular  definition  of 
SC  •  Equation  (57),  on  the  other  hand,  indicates  that  AC  is  automatically  zero 
at  epoch  (t  ^  0);  AC  is,  in  fact,  the  definite  integral  of  Cfp  f;:om  epoch  to 
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t  ,  wh«c««s  8C  is  esssntislly  an  indafinif  integral.  Many  writers  have  worked 
with  epoch-independent  long-period  perturbations,  so  that  dC  behaves  like  S^, 
but  as  the  present  author  has  reaiarked  before^,  a  heavy  price  has  to  be  paid,  in 
particular  in  the  introduction  of  an  intrinsically  spurious  singularity  at  the 
so-called  critical  inclinations  (where  tt,ec  vanishes) .  Nothing  is  lost  by 
having  epoch-dependent  dC  ,  since  this  simply  makes  the  long-period 

A 

perturbations  behave  like  the  secular  perturbations,  Cjec  t  being  necessarily 
epoch-dependent . 

laplicit  in  the  last  paragraph  is  the  understanding  that  a  short-period 
perturbation  la  only  defined  'up  to  an  arbitrary  constant',  and  this  is  inherent 
in  the  original  (54),  in  the  concept  of  an  indefinite  integral,  in  the  notion  of 
epoch  independence,  and  in  the  earlier  reference  to  a  'particular  definition  of 
55  ' .  The  choice  of  a  particular  constant  in  ,  or  the  choice  of  a  particular 
'reference  value'  for  (  (which  amounts  to  the  same  thing),  is  a  matter  of 
considerable  inportance  at  the  fir.st-order  level  In  J2  ,  and  it  is  unfortunately 
not  always  clear,  from  published  papers,  what  constants  the  authors  have  chosen 
or  the  motivation  for  their  choice.  As  was  remarked  before^  (and  more  recently 
in  Ref  6),  there  are  several  possible  rationales.  Thus,  the  constants  may  be 
chosen  to  auike  each  unbiased  in  time,  ie  such  that  its  value  is  zero  wlien 
averaged  with  respect  to  t  .  This  choice,  which  the  use  of  the  epithet  'mean' 
rather  naturally  suggests,  leads  to  awkward  constants  (see  equations  (194)  to 
(198)  of  Ref  3),  so  a  better  one  is  to  make  8C  unbiased  in  respect  of  the 
transformed  Integration  variable  v  .  The  best  strategy  of  all,  however,  is  to 
make  an  overall  choice  of  constants  so  as  to  get  the  sinplest  possible 
expressions  in  a  particular  set  of  quantities  derived  from  the  elements.  This 
was  the  strategy  of  Refs  2  and  3,  and  the  present  Report  continues  the  philo¬ 
sophy,  the  'particular  set  of  quantities'  t>elng  the  short-period  perturbations  in 
the  system  of  spherical-polar  coordinates  to  Ise  introduced  in  section  3.3. 

In  Ref  3,  explicit  arbitrary  constants  were  introduced  into  the  expressions 
for  5C  derived  from  first-order  J2  analysis,  so  that  the  chclces  adopted  by 
particular  authors  could  be  conveniently  quoted  and  conpared.  This  Involved  a 
generality  that  will  not  be  repeated  in  the  present  paper,  and  in  section  4  we 
use  the  'best'  constants  from  the  start.  In  the  present  section  we  must  briefly 
consider  just  one  of  these  constants,  the  one  associated  with  semi-major  axis, 
and  here  the  obvious  (and  best)  choice  recovers  the  energy  constant  a'  intro¬ 
duced  In  section  2.3,  this  being  sometimes  referred  to  as  the  semi -major  axis  of 
Brouwer^^.  It  was  used  at  RAE  until  the  advent  of  the  program  PROP^^"'®,  for 
which  the  decision  was  taken  -  mistakenly  with  hindsight  -  to  follow  the 
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Smithsonian  Astrophysical  Obsarvatory  by  switching  to  ths  semi-major  axis  of 
Koaail''. 

The  choice  of  i  also  bears  upon  the  choice  of  n  .  since  this  must  at 
first  sight  be  tied  to  i  by  a  barred  version  of  (1),  in  which  case  (53)  demands 
that  n  ••  n*  .  Kepler's  third  law  does  not  have  to  be  carried  over  in  this 
unmodified  way,  however^,  and  in  spite  of  the  vital  role  played  by  n'  ,  there  is 
a  simpler  and  more  natural  choice  of  n  .  To  arrive  at  this  choice,  we  set  C 
to  H  in  the  generic  (S')),  getting 

H-M.+M  t  +  AM,  (58) 

0  *ac  '  ' 


where 


"  o. 


(59) 


and 


AM  >  Ao  -I-  aJ  . 


(60) 


Here  /  is  the  shorthand  quantity  introduced  by  (3) ,  and  we  can  therefore  write 

t 

I  "  n't  +  J  (n  -  n')  dt  . 

0 

From  (1)  and  (53)  we  can  relate  n-n'  to  8a(“n-a')  and  then  integrate 
using  (40) .  (The  notation  fin  has  not  been  used  for  n  -  n'  ,  since  it  more 
appropriately  applies  to  n  -  n  .)  The  integral  in  (61)  has  a  secular  component, 
(ab[t')t  say,  such  that  n'  +  dn'  constitutes  .  Then  the  natural  definition 

of  n  is  given  by 


3  + 

sec 


n'  +  dn' 


(62) 


If  n  is  chosen  in  this  way,  it  becomes  part  of  the  task  of  orbit  analysis  to 
express  Kepler's  third  law  In  an  appropriately  modified  form^  such  that  n  can 
be  correctly  derived  from  a  . 

Me  shall  find,  in  section  4,  that  in  the  first-order  J;  analysis  n  is 
identical  with  n'  ,  so  that  the  Kepler  law  is  valid  for  mean  elements  without 
modification  -  in  fact  it  Is  (53).  in  section  5,  however,  we  shall  find  this  to 
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b«  no  longer  true  at  the  sc  i:on<l-order  level;  n  differs  from  n'  ,  and  the 
modified  Kepler  law  is  givon  by  (239) . 

It  is  natural,  and  sonvenient,  to  extend  the  concept  of  'mean*  from  the 
standard  orbital  element:^  to  other  quantities,  eg  v  and  r  .  The  basis  for 
definition  is  that  the  ‘Expression  for  a  mean  quantity  in  terms  of  mean  elements 
is  the  same  as  the  exp  'csslon  for  the  original  quantity  In  terms  of  osculating 
elements.  Thus  th«  e.^uations  of  section  2.1  hold  true  after  the  addition  of  bars 
to  all  quantities,  the  only  (possible)  exception  being  (1)  -  for  the  reason  that 
has  just  been  given .  As  an  exan^le,  (5)  yields  * 

p/r  -  1  +  e  cos  V  ,  (t>3) 

after  which  the  quantity  W,  anticipated  at  (39),  becomes  available.  (Obviously, 

V  and  r  are  nit  like  the  ^  in  being  free  of  short-periodic  variation!) 

3.2  SeaUL-atetn  elesMnts 

An  ia^ortant  ccnplicacion  arises  in  regard  to  the  derivation  and  use  of  the 
mean  elements,  (  ,  and  their  smooth  (free  of  short-period  variation)  rate  of 
cl'ange,  f  ;  it  relates  to  the  change  of  variable  from  c  to  v  when 
integrating  the  planetary  equations,  (22)-(27).  The  integration  of  s^Ecular  and 
)ong-period  terms  in  d^/dv  leads  to  long-term  variation  that  is  'smooth  in  v  ' 
rather  than  'smooth  in  t  '  ,  whilst  thr  integration  of  short-period  terms  leads 
to  a  pure  Poisson  series  (le  to  a  combination  of  trigonometric  terms  with 
arguments  linear  in  v  and  U  ) .  We  can  think  of  this  long-term  variation  as 
applying  to  a  9»mi-mean  element  ^  ,  rather  than  to  the  true  mean  element,  the 
relation  between  the  two  being  of  the  form 

C  -  C  +  CS  ;  (64) 

here  in  is  the  'mean  equation  of  the  centre*  defined  in  accord  with  (8),  wnilst 

A  » 

C  is  essentially  of  the  form  ^/n  as  wo  shall  see.  (The  qualification 
'essentially'  is  because  the  precise  def.inition  we  adopt,  as  most  useful  in 
practice,  Involves  only  the  secular  component  of  ,  being  given  by 
equation  (67),  but  this  is  irrelevant  for  the  moment.)  To  justify  (64),  suppose 

A 

that  C  has  been  introduced  to  denote  the  long-term  component  of  dC/dv  , 
defined  from  the  right-hand  side  of  the  planetary  equation  for  t  after  the 
change  of  variable.  This  integrates  (at  least  for  a  purely  secular  long-term 
variation)  to  C  -  Cg  r  given  by  C  (v  ”  vg)  .  But  this  may  be  rewritten  as 
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5(M  -  Mq)  +  C(m  -  Mo)  ,  whore  the  first  tenr  is  no*  smooth  in  t  and  is  to  be 
identified  with  Ct  •  This  term  constitutes  C  ”  Co  '  so  we  have  derived  (64) 
in  the  torn:  specifying  variation  from  to  to  t  ;  the  second  term  m^iy  be 
thought  of  as  an  induced  -carry-over*  term. 

It  follows  from  (54)  and  (64)  that 

c  -  c  -  8C  -  CS  ,  <65) 

A 

where  in  principle  (until  C  has  been  precisely  defined)  (65)  expresses  the 
pure-Poisson  component,  8pC  say,  of  the  short-period  perturbation,  to  which 
Cfn  must  be  added  to  give  the  complete  5C  .  Now  we  are  concerned  with  a 
relation  (operating  to  second  order)  of  the  form 

C  =  C  +  {^n)m  +  SpC  ,  (66) 

in  which  the  right-hand  side  reduces  to  a  two-term  sum  if  each  component  of  the 
middle  (induced)  term  is  combined  with  either  C  (contributing  to  C  )  5p^  . 

In  first-order  J2  analysis,  the  only  induced  terms  to  arise  are  those  associated 

A 

with  iljec  ^sec  (because  of  the  identity  between  n’  and  n) ,  and  it  is 

natural  to  Incorporate  them  with  O  and  a  ,  respectively,  leaving  just  the 
pure  6pC  to  be  amalgamated  in  the  form  of  perturbations  in  spherical 
coordinates.  On  proceeding  to  second  order  (as  in  section  5),  we  find  this 
first-order  policy  to  have  been  essential  (and  not  merely  'natural'),  at  any  rate 
for  0}  ,  since  aw)cward  carry-over  terms  can  only  be  avoided  if  &  ,  not  5), 
appears  in  the  arguments  of  the  first-order  terms.  It  is  then  obviously  sensible 
for  n  and  c5  to  include  the  second-order  terms  induced  by  ftsoc  sf'd  Ujoc  » 
as  well  as  the  first-order  ones,  and  it  will  be  found  in  section  5  that  there 
also  has  to  be  a  second-order  term  in  M  ,  Induced  by  Ojec  c(n'  from  (62);, 
there  is  no  effect  from  n*  (the  dominant  term  in  (62))  because  this  integrates 
directly  to  n't  ,  without  any  need  to  transform  the  variable  from  t  to  v  . 

The  Induced  effects  of  the  long-period  variation  could  be  treated  in  the  same  way 
(and  would  probably  have  to  be  if  the  analysis  were  being  ta)cen  to  third  order) , 

but  as  arises  for  every  element  (except  a  )  it  is  simplest  if  the  entire 

long-period  con^onent  of  the  middle  term  of  (66)  is  combined  with  6pC  rather 

than  d  ;  the  resulting  quantity  may  be  denoted  by  8^  ,  and  it  is  the  8?  that 

are  amalgamated  into  coordinate  perturbations. 
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Thus  th«  precise  definition  of 

A 

c  - 


t  is 


(67) 


except  that  Mtac  snist  be  replaced  by  n  -  n*  Mhen  C  is  M  .  With  this 
definition,  we  can  rewrite  (66)  as 


c  -  c  +  . 


(68) 


We  also  require  the  tine  derivative  of  (64).  With  m  defined  by  (8),  v 

*  A 

given  by  (39),  M  ta)cen  as  n  ,  and  C  ta)cen  as  constant,  this  derivative  is 
given  at  once  in  the  form 


C  -  C  +  Cn  (W  -  1) 


(69) 


A 

Now  (  is  given  by  (67) .  Also,  it  is  legitimate  and  convenient  to  replace  w 

by  W  ,  since  they  are  the  same  to  0(J2)  and  (  is  0(J2)  ,  so  we  can  rewrite 
(69)  as 


i 


(70) 


except  that  (70)  needs  the  extra  term  n' (1  ~  H)  when  (  is  M  . 

Explicit  results  for  the  perturb.ttlons  in  the  orbital  elements,  associated 
with  J2,  J2^  and  J3  respectively,  are  derived  in  sections  4,  5  and  7  of  this 
paper,  and  it  is  convenient  to  represent  the  three  sets  of  expressions  in  terms 
of  a  compact  notation  Involving  the  suffixes  1,  2  and  3.  Thus,  for  the  secular 
perturbations  we  write 


C  "  K  Cl  +  C2  (+  H  C3) 


(71) 


(all  the  C3  in  fact  being  zero) ;  for  the  long-period  perturbations  we  write 
(since  the  Ci.fp  are  all  zero) 
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(72) 

and  for  the  pure  short-period  perturbations  we  write 

8pC  “  K  Cl  +  ”  ^3  ■ 

It  then  follows  that 

5C  -  K  Cl  +  k'(C3  +  +  H(C3  +  C3,<pS)  ,  (74) 

it  being  iinmaterial  whether  we  write  m  or  in  in  (74)  . 

3 . 3  Satellit*  position  and  velocity 

Given  a  set  of  mean  orbital  elements  at  epoch,  the  unimaginative  and 
laborious  way  of  computing  the  satellite's  position  and  velocity  at  time  t  is 
to  incorporate  the  appropriate  perturbations,  long-term  and  short-period,  in  each 
element,  and  then  determine  position  and  velocity  from  the  resulting  osculating 
elements.  The  algorithm  for  determining  position  from  osculating  elements  was 
given  before^*^  (see  also  Ref  7),  but  for  completeness  it  is  given  again  here, 
with  the  necessary  extension  for  generation  of  velocity,  after  which  we  return  to 
the  use  of  mean  elements.  The  standard  system  of  geocentric  coordinates 
(x,  y,  z) ,  is  assumed,  with  x  measured  towards  the  vernal  equinox  and  z 
towards  the  north  pole.  Then  the  algorithm  for  position  and  velocity  from 
osculating  elements  (C)  is  as  follows: 

(i)  the  eccentric  anomaly,  £  ,  is  found  by  solving  Kepler's  equation^®- 

E  -  e  sin  E  -  M  ;  (75) 

(ii)  V  is  found  from  one  of  the  two  equivalent  formulae  (apart  from  the 
ambiguity  of  quadrant  in  the  first  formula,  >;hich  is  automatically  resolved 
if  the  Fortran  ATAN2  function  is  used) 

.  q  sin  E  ,T,. 

tan  V  -  '  7  "  (76) 

cos  E  -  e 

and 

1  .ft 

tan  )4v  -  tan  >4E  ;  (77) 
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(iii>  u  is  obtained  from  (7),  and  r  from  either  (5)  or  the  equivalent  formula 

r  a{l  -  e  cos  E)  ;  {78) 

(iv)  n  is  obtained  from  (i),  and  then  (38)  gives  v  ,  with  which  we  can 

identify  u  (the  dot  notation  is  legitimate,  because  we  are  referring  to  a 
frozen  osculating  orbit,  for  which  there  are  no  perturbations,  so  that 
dv/dt  -  dv/Bt  and  du/dt  “  du/3t) ; 

(V)  r  is  obtained  from  one  of  the  equivalent  formulae 


and 


r  =  naeq"^  sin  v 


f  -  ne  (3^/r)  sin  E  ; 


(79) 

(80) 


(vi)  X,  y  and  z  are  obtained  from  the  double  coordinate  transformation 
expressed  by  the  matrix  fontula 

(x  y  z)^  "  R3(-Q)  Ri(-i)  (t  cos  u  r  sin  u  0)^  ,  (81) 

where  Rj(  )  describes  rotation  about  the  axis,  so  that 


Rl<0)  -  /  1  0  0 

0  cos  0  sin  0 

\ 0  -sin  0  cos  0 


and  T  in  (81)  denotes  transposition; 


(vii)  X,  y  and  z  are  obtained  by  3-differentiation  of  (81),  which  gives 

.  ,  .  X  ,  .  T 

(x  y  z)  "  R3(-fl)  Ri(-i)  (r  cos  u  -  rii  sin  u  r  sin  u  +  rii  cos  u  0)  . 

. (83) 

If  Mo  ,  the  conceptual  value  of  the  mean  anomaly  at  epoch  (though  defined  at 
time  t  )  is  the  starting  point,  rather  than  the  mean  anomaly  at  t  itself,  then 
there  is  of  course  a  preliminary  step  given  by  M  -  Mq  +  nt  . 
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When  we  start  from  a  set  of  I^q  (as  opposed  to  Co  )  and  apply  (57), 
followed  by  (54),  to  obtain  the  C  from  which  the  foregoing  algorithm  can  be 
applied,  then  much  the  most  involved  part  of  the  procedure  arises  witn  the 
confutation  of  the  8C  to  substitute  in  (54)  (or  in  practice  the  to 

substitute  Jn  (68))  .  Hence  the  overall  procedure  will  l>e  greatly  sinflified  if 
we  can  operate  the  position-velocity  algorithm  in  a  modified  form,  starting  (at 
tine  t  )  from  t)ie  C  (^r  alternatively  the  C  )  instead  of  the  C  ■  '^^e 
modified  algorithm  would  then  lead  to  ’ (semi-) mean  position  and  velocity',  after 
which  it  would  (hopefully)  be  possible  to  add  greatly  simplified  combinations  of 
the  5C  directly  to  the  position  and  velocity  components. 

It  was  in  pursuance  of  this  philosophy  that  Kozai’^  obtained  first-order 
short-period  perturbations,  5r  and  5u  ,  associated  with  J2  ,  in  a  much  more 
compact  form  than  the  perturbations,  5a,  6e,  8(0  and  Sm  ,  that  they  displaced. 
These  perturbations  (5r  and  8u)  were,  in  effect,  the  perturbations  in  a  pair  ol 
two-dimensional  polar  coordinates  (r  and  u),  the  plane  of  the  coordinates  beirg 
simply  the  osculating  plane  of  the  orbit.  Much  more  recently,  it  has  been 
pointed  out  by  the  present  author^' ^  that  it  is  logical  to  complete  Kozii's 
approach  by  ma)cing  the  coordinate  system  three-dimensional,  rotating  but  with  a 
standard  geocentric  origin,  with  a  reference  plane  that  is  now  the  'mean  oroital 
plane'  instantaneously  defined  by  i  and  Q  rather  than  i  end  Q  .  (It  is 
the  'seirl-mean  orbital  plane’,  based  on  Q  rather  than  O  ,  that  we  strictly 
require,  but  this  may  be  regarded  as  an  unnecessary  complication  at  this  point.) 
On  this  basis  we  can  interpret  a  set  of  mean  orbital  eleme-uts  at  time  t  a.? 
defining  tno  position  of  a  'mean  .satellite',  having  mean  anomaly  M  in  a  'moan 
orbit'  defined  (instantaneously)  by  a,  e  and  ca  ;  the  mean  orbit  lies  in  (and 
moves  with)  the  meai’  orbital  plane.  This  is  a  helpful  point  of  view  for 
geometric  visualization,  but  care  is  needed.  Thus  tne  r.can  satellite  itself  ha.", 
a  set  of  evolving  osculating  elements,  and  these  cannot  be  identified  with  the 
mean  elements  of  the  actual  satellite,  in  particular  necau.se  the  osculating 
elements  of  the  mean  satellite  are  not  (in  genera) )  free  of  short-period 
variation  -  the  underlying  explanation  of  this  is  chat  the  Instantaneous  mean 
elements  only  define  the  position,  not  the  velocity,  of  the  utan  satellite*. 

*  For  satellite  motion  i:,  a  unitormly  rotating  Kepierian  e.il.ipsc:  in  a  uniformly 
rotating  orbital  plane  (non-conservative  example,  in  ge»ieral ) ,  mean  elements 
can  be  defined  such  that  the  'mean  satellite'  actually  coincides  with  th?  true 
satellite;  then  i  identifies  with  i  (osculating)  at  the  'apexes',  but.  not 
with  i  at  the  'nodes'.  (This  example  underlies  the  c.noice  of  I  in  fir.se- 
otder  J2  analysis  -  see  section  1.) 
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Refs  2  and  3  descriiaed  a  conflation  of  the  Kozai  approach  in  which  the  axis 
system  is  referenced  by  cylindrical  polar  coordinates,  the  coordinates  within  the 
mean  orbital  plane  (corresponding  to  r  and  u  )  being  designated  r'  and  u'  , 
and  the  out-cf-plane  coordinate  c  .  It  would  have  been  just  as  satisfactory  to 
have  used  spherical  polar  coordinates,  however,  and  in  the  complete  (untruncated) 
analysis  presented  here  it  turns  out  that  the  spherical  system  is  superior.  The 
radial  coordinate  of  the  (actual)  satellite  in  this  system  then  naturally  reverts 
to  r  ,  and  the  other  two  coordinates  will  be  denoted  by  b  and  w  ,  the 
(vector)  relationship  of  r,  b  and  w  to  rectangular  coordinates,  X,  Y  and 
Z  ,  based  on  the  mean  orbital  plane  as  XY-plane  (with  X  towards  the  node) , 
being 

0“  /r  cos  b  cos  w\  .  (84) 

I  r  cos  b  sin  w  j 
\  r  sin  b  / 

The  directions  of  the  axes  of  X  and  Y  are  not  well  defined  for  a  near- 
equatorial  mean  orbital  plane,  but  it  will  be  found  in  due  course  that  (so  far  as 
position  is  concerned)  this  causes  no  difficulty.  The  relationship  of  the  new 
spherical  system  (r,  b,  w)  to  the  (obsolete)  cylindrical  system  is  immediate, 
since  in  the  old  system  we  have  (X,  Y,  Z)  -  (r*  cos  u',  r'  sin  u',  c)  ,  but  in 
view  of  the  close  relationship  between  the  out-cf-plane  coordinates  b  and  c 
(since  sin  b  «  c/r  )  it  is  worth  an  explicit  observation  that,  whereas  c  was 
the  third  cylindrical  coordinate,  b  is  the  second  spherical  one.  (Clearly, 
also  w  "  u'  and  r^  "  r'^  +  c^  .) 

The  position  of  the  mean  satellite  is  conceptually  given  by  coordinates 
r,  b  and  w  ,  but  it  is  really  the  'semi -mean  satellite'  we  are  concerned  with, 
so  the  reference  point  for  (short-period)  perturbations  must  actually  be  expressed 
as  (r,  b,  w) .  Moreover,  since  Z  is  zero  for  the  semi -mean  satellite,  we  have 
from  (84)  that 

b  -  0  (85) 


and  it  is  also  clear  that 


(86) 


Since  the  whole  point  of  this  approach  is  that  we  apply  (68)  to  coordinates, 
rather  than  elements,  the  spherical  coordinates  of  the  actual  satellite  are  given 
by 
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(r,  b,  w)  -  (i,  b,  w)  +  (Sr,  Sh,  Sv)  ,  (87) 

with  a  natural  extension  of  the  notation  that  (68)  introduced.  From  (87)  we  also 
have 

(r,  b,  w)  -  (r,  b,  w)  +  (Sr,  Sb,  Sit)  ,  (88) 

on  the  understanding  that  Sr  (etc)  denotes  d(5r)/dt  ,  rather  than  5(dr/dt); 
the  distinction  arises  because  (in  particular)  r  is  not  identical  with  r  , 
assuming  the  latter  to  he  defined  by  converting  the  right-hand  side  of  (80)  to 
semi -mean  elements. 

We  can  now  focus  on  the  specific  objective  of  this  section,  namely,  the 
derivation  of  position  and  velocity  at  time  t  fro.n  mean*  elements  at  epoch 
without  any  reference  to  osculating  elements.  The  derivation  of  mean  elements  at 
t  is  necessary  as  a  preliminary  step,  carried  out  in  principle  by  application  of 
(57);  the  practical  difficulties  in  this  step  are  considered  in  section  3.5  (and 
further  in  Part  3),  so  we  may  assume  here  a  starting  point  consisting  of  the  six 
C  at  t  .  We  confine  ourselves  initially  to  the  derivation  of  position,  and 
start  by  applying  steps  (i)  and  (ii)  of  the  algorithm  that  was  given  for  oscula¬ 
ting  elements;  the  result  is  the  'mean  true  anomaly*  v.  Since  in  is  now  avail¬ 
able,  (64)  gives  the  appropriate  ^  »  ’'Ir  fl,  5  and  M  ;  then  a  repeat  of  steps 
(i)  and  (il)  (only  done  to  reflect  the  change  from  M  to  M  ),  followed  by  (iii), 
leads  to  v,  r  and  u  .  At  this  point,  in  view  of  (85)  and  (86),  we  have  the 
semi-mean  satellite  position  (r,  b,  w)  ;  on  the  assumption  that  formulae  for 
Sr,  Sb  and  5w  are  available,  it  follows  that  (87)  yields  (r,  b,  w),  and  then 
(84)  gives  (X,  Y,  Z) .  We  finally  require  the  formula  corresponding  to  (81)  in 
step  (vi)  of  the  algorithm,  and  thi  may  be  written 

(X  y  z)”^  -  R3(-ft)  Rj(-i)  (X  y  Z)’'  ,  (89) 

since  there  is  no  distinction  between  i  and  i  . 

It  is  a  straightforward  matter  to  express  the  generic  first-order  formulae 
for  Sr,  Sh  and  5w  in  terms  of  the  S^  ,  and  we  now  do  this  in  advance  of  the 

specific  first-order  J2  analysis  of  section  4  and  the  J3  analysis  of  section  7. 

*  To  avoid  confusion,  it  should  be  noted  that  the  basic  (reference)  elements  are 
mean,  not  semi -mean;  in  application  of  the  ^position-velocity  algorithm', 
however,  the  first  requirement  is  that  the  C  be  superseded  by  the  ^  i  and  it 
Is  the  latter  that  define  the  spherical-coordinate  system. 


TR  88068 


28 


(Second-order  formulae  are  much  harder  to  expresa/ 
ia  deferred  «intil  they  are  needed,  in  aection  6.) 
Identification  of  the  two  formulae  for  (x  y 
Uaing  alao  (84)  in  the  identification,  we  get 


however,  and  their  derivation 
Our  atarting  point  ia  the 
z)'^  ,  viz  (61)  and  (89). 


(coa  b  coa  w 
coa  b  ain  w 
ain  b 


) 


(I)  R,  (-i) 


(90) 


Now 


R3(n)  R3(-a)  -  R3(-5a)  , 


(91) 


exactly,  but  firat-order  approximation  ia  required  to  complete  the  reduction, 
which  leada  to 


and 


b  -  ain  u  5i  -  a  coa  u  SCI 

w  -  u  +  c  SCi 


(92) 


(93) 


In  view  of  (85)  we  can  write  ^b  in  place  of  b  in  (92),  and  in  view  of  (86), 
which  can  be  aubtracted  from  (93) ,  we  alao  get 

5w  ■  5u  +  c  50  .  (94) 


For  the  formulae  we  are  aeejcing,  it  only  remaina  to  expreaa  5r  and  5u 
in  terma  of  the  5C  /  and  the  req^uired  expreaaiona  are  immediate  from  (41)  and 
(42) .  We  thua  obtain  (genetically) 


- 

(r/a)  5a  - 

(a  coa  v)  5e 

+  (aeq'*  ain  v)  5m  , 

(95) 

Sb  - 

(ain  u)  5i 

-  (a  coa  u) 

50 

(96) 

and 


5w  -  5o)  -)■  c  50  -)■  q*2  sin  v  (2  +  e  coa  v)  Je  +  w5m  . 


(97) 


In  practice,  of  courae,  the  coefflcienta  in  (95) -(97)  are  to  be  interpreted  as 
aend-mean  quantitlea  and  not  aa  oaculating  onea. 
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It  was  remarked,  earlier  in  the  section,  that  no  difficulty  in  position 
computation  arises  from  the  indeterminacy  (for  near-equatorial  orbits)  of  the 
X-axis  (origin  for  w  )  in  the  rotating  system  of  spherical  coordinates,  and  this 
needs  to  be  justified  before  we  pass  on  to  velocity  (for  which  there  is  a 
difficulty  -  to  be  dealt  with  in  the  next  section) .  The  underlying  principle 
here  is  that  a  set  of  standard  elements  (osculating,  mean  or  semi -mean)  can 
always  be  safely  used  (for  position  computation) ,  so  long  as  their  values  are 
consistent.  All  this  means  is  that  no  matter  how  badly  defined  the  elements 
n,  (0  and  M  are  individually,  no  accuracy  is  lost  so  long  as  they  are  such 
that  their  non-singular  combinations  (cf  the  quantities  y  and  L  of  section  2) 
retain  full  accuracy. 

Given  a  consistent  set  of  elements  at  epoch) ,  the  first  step  in 

the  computation  of  satellite  position  at  time  t  is,  as  already  noted,  the 
application  of  long-term  perturbations,  to  produce  t  at  t  .  This  must  be  done 
with  full  allowance  for  the  possibility  of  singularity,  as  described  in 
section  3.5,  after  which  we  have  (at  time  t  )  a  consistent  set  of  mean  elements 
and  hence  a  consistent  set  of  semi-mean  elements.  If  Q  is  not  well  defined  in 
this  set,  then  the  X-axis  in  the  semi-mean  orbital  plane  must  also  be  ill 
defined;  but  this  will  not  matter,  assuming  that  the  value  of  (b  ,  to  be  added 
to  V  to  give  u  (  -  w)  ,  is  consistent  with  ft  .  In  other  words,  the 
coordinates  (X  and  Y)  given  by  (84)  will  only  be  arbitrary  to  an  extent  that  is 
precisely  compensated  for  in  the  application  of  the  corre.spondingly  arbitrary 
ft  in  (89) .  This  is  true  even  in  the  extreme  case  -  which  can  arise  in  practice 
as  will  be  seen  in  section  7.4  -  when  ft  becomes  infinite,  but  such  an  infinity 
inevitably  has  a  more  profound  effect  on  velocity,  which  is  why  the  velocity 
formulae  to  be  presented  at  the  end  of  this  section  are  only  provisional. 

Thus  a  nominal  indeterminacy  for  w  in  (r,  b,  w)  is  dealt  with  easily 
enough.  For  the  perturbations  (5r,  5b,  5w )  the  situation  is  even  better,  since 
these  are  well  determined  at  all  times.  This  is  essentially  because  of  the  5ft 
term  in  (94),  which  reflects  the  fact  that  5w  (unlike  5u  )  is  tied  to  the  same 
reference  direction  as  w  itself  (Kozai's  5u  is  free  of  'circular  singularity' 
but  not  of  'equatorial  singularity'). 

We  now  proceed  to  the  formulae  for  x,  y  and  z  ,  to  complete  the  algorithm 
for  position  and  velocity,  bearing  in  mind  that  the  formulae  are  only 
provisional.  The  adjustments  to  avoid  infinities  are  given  in  the  next  section. 

We  have  to  cover  the  three  remaining  steps  -  (iv) ,  (v)  and  (vii)  -  in  the 
adaptation  of  the  original  algorithm  to  the  use  of  (semi-) mean  elements.  It 
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£o11oms  from  the  way  in  which  r  is  defined  that  we  can  invoke  (41)  again,  as  in 
the  derivation  of  (95),  to  obtain  r  ;  since  a  ••  0  we  get 

•  -1  i 

r  -  (a  e  q  sin  v)  M  -  (a  cos  v)  e  .  (98) 

Clearly  S  is  zero.  Finally,  (42)  can  be  invoked  again,  with  (86),  to  give  w  , 
the  result  being 


_-2 


“  «a  +  W  M  +  q  sin  v  (2  t  e  cos  v)  e 


(99) 


The  absence  of  a  term  in  Q  (cf  (97))  is  due  to  the  point  just  made  in  relation 
to  Su  and  w  ;  is  allowed  for  separately  (as  we  are  about  to  see)  via  the 
motion  of  the  (aeml-)mean  orbital  plane,  whereas  SQ  is  resolved  directly  into 
contributions  5b  and  5w  . 

From  <r,  b,  w)  we  obtain  (r,  b,  w),  by  appeal  to  (88),  (.Si,  5b,  Si) 
being  given  by  differentiation  of  (Sc,  5b,  5w).  The  resulting  short -period 
expressions  Involve  combinations  of  v  and  &  ,  including  u  (  “  v  +  &)  ,  but 
first-order  representations  of  these  are  adequate  for  our  second-order  analysis, 
so  the  singularities  associated  with  long-period  perturbations  in  Q  and  (O 
cause  no  difficulties  in  these  expressions. 


It  remains  to  obtain  a  formula  having  the  same  relation  to  (83)  as  (89)  has 
to  (81) .  Differentiating  (89)  directly,  we  get 


y 

z 


{R3(-a)  Ri(-I)  I  +  Rj  (-n)  Bi(-I)n}  /  +  RjC-O)  RjCI)  j  x 


(100) 


where  Rj  (-0)  “  dRj(-6)/d0 


0  0  0  \  (etc) 

0  -  sin  0  ~  cos  0  j 

0  cos  0  -  sin  0  ^ 


(101) 
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Now  from  (84)  wo  hove 


(coa  b  coa  w  -  aln  b  coa  w 

coa  b  ain  w  -  ain  b  aln  w 

aln  b  coa  b 


Alao,  alnce  1  **  0  ,  we  have 

aoc 


whilat 


11.  +  w  n 

fp  soc 


(103) 


(104) 


by  (70) .  It  may  be  obaerved  that  the  coefficient  of  i  in  (100)  can  be 
expreaaed  *a  (z  ain  11  -  z  coa  11  cY  -  aZ)’’’  ,  whilat  the  coefficient  ot 

11  can  be  expreaaed  aa  (-  y  x  0)'*^  . 

3 . 4  Awoldanco  of  olngulazlty  In  voloclty  computation 

The  formulae  at  the  end  of  the  laat  aection  being  unauitablc  for  univocaal 
velocity  computation  aa  they  atand,  becauae  of  poaaible  singularity,  we  must 
start  by  identifying  the  potentially  infinite  quantities.  They  are  tb  and  M 
in  (99),  which  propagate  into  w  in  (102),  and  the  compensating  H  in  (100). 

The  quantity  M  la  potentially  singular  through  expressions  for 
that  contain  e  ^  aa  a  factor  -  see  (399)  in  particular.  This  causes  no  problem 
in  (98),  where  fi  occurs  with  a  multiplying  e  ,  but  in  (99)  the  potential 
infinity  is  real.  It  is  convenient  to  deal  with  the  two  formulae  in  the  same 
way,  however,  and  we  seelc  to  replace  M  by  as  close  a  non-singular  equivalent  as 
possible. 

From  the  special  case  of  (70)  when  la  M  , 


(106) 


32 


W*  now  recall  that  £  ,  given  by  (35),  cannot  b«  aingular,  so  we  introduce  a  new 
quantity,  N  (designated  nj,  in  Ref  4),  that  can  Ise  regarded  as  a  modified 
version  of  the  swan  motion  (though  it  is  not  free  of  short-period  variation,  let 
alone  constant);  it  is  defined  by 

N  -  n  +  (i  -  1)  (n  -  n')  +  L^p  .  <107) 

Then  from  (106)  and  (35)  we  have 

M  -  N  -  q^jp  .  (108) 

Hence,  since  i  0  ,  (38)  can  be  rewritten  as 

r  «  (a  e  q  ^  sin  v)  N  -  (i  cos  v)  e^p  -  (a  e  sin  v)  y^p  .  (109) 

All  terms  here  are  free  of  singularity,  the  last  term  (li)ce  the  original  M  term) 
being  so  because  of  the  e  factor. 

Our  real  concern  is  to  rewrite  (99) ,  and  this  introduces  the  complication 
that  (as  we  have  seen)  the  quantity  on  the  left-hand  side,  viz  w  ,  itself  has 
the  potential  singularity,  so  that  a  modified  quantity  must  be  expressed.  From 
(70),  with  C  ta)cen  as  <o  ,  and  (108),  wo  have 

fi  +  HM  -  w(n  +  -  qy^p)  +  .  (110) 

in  which  the  last  term  is  potentially  singular  for  (both)  circular  and  equatorial 
orbits,  but  even  the  term  in  y^p  is  no  longer  singularity-free.  However,  by 
addition  of  c  tiip  to  both  sides  of  (110) ,  we  may  write 

5  +  it  M  +  c  fi^p  "  W  (n  +  -  (q  W  -  1)  y^p  ,  (111) 

and  now  the  right-hand  side  is  well  behaved,  since 

5*"1  “  {2  cos  v  +  e(l  +  cos^  v)  }  .  (112) 
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It  -jmedlatsly  fol^.rjs  that  the  modified  (non-singular)  version  of  (99)  that  we 
require  is 


w*  “  5  1-  c  ft,  .  (114) 

«p 

(It  is  remarjced  that  ft*  is  the  appropriate  quantity  to  use  as  the  derivative  of 
a  ,  when  differentiating  the  trigonometric  arguments  of  the  individual  terms  of 

Si,  and  to  provide  Si,  Sh  and  Sit  ;  (86)  suggests  it  is  w  rather 

than  w*  that  is  required,  but  w  is  not  singularity-free;  also,  the  difference 
between  w  and  w'  can  be  neglected  in  Si  etc,  since  the  resulting  effects 
are  0(J2^,  J2J3)  . ) 

Since  w',  rather  than  w  ,  is  the  non-singular  quantity,  it  follows  chat 
X  and  Y  ,  as  given  by  (102),  should  be  replaced  by  quantities  computed  from  w' 

rather  than  w  ,  where  w'  is  obtained  from  addition  of  Sit  to  w'  instead  of 

to  w  .  This  will  ma)ce  the  last  term  of  (100)  non-singular  and,  as  we  shall  see, 
in  such  a  way  as  to  ma)ce  the  middle  term  non-singular  as  well .  What  we  are  doing 
to  (102)  amounts  to  the  creation  of  errors  -r  c  ft^p  cos  b  sin  w  and 

A 

cos  b  cos  w  ,  in  X  and  Y  respectively,  that  can  be  corrected  in  (100) 
by  the  additional  term 


-  -  .  .  .  T 

on  expanding  (-ft)  Rj  (-i)  ,  the  multiplier  of  (X  Y  Z)  .  But  the 
additional  term  can  be  rewritten  as 
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cos  b  cos  w 


c  sin  Q  cos  Q  0 
c  cos  Q  sin  Q  0 


cos  b  sin  « 


which  is  in  just  ths  right  foru  for  combination  with  the  coiqionent  of  the 

middle  term  of  (103),  which  can  itself  be  written  as 


sin  Q  -  c  cos  O  s  cos  Q 

cos  fi  -  c  sin  Q  s  sin 


cos  b  cos  w 
cos  b  sin  w 


On  combination,  and  in  view  of  (84),  we  get 


I  scosii  0  sinQ 


This  can  be  rewritten  as  s  Q.  R,  (-0)  R,  (-1)  (Z  0  -X)  ;  the  appearance 

«P  3  1  ^ 

of  8  ,  effectively  replacing  c  as  a  multiplier  of  ,  ma)ces  the  expression 

fully  non-singular,  and  its  form  is  such  that  it  can  be  gathered  back  into  the 
last  term  of  (100) . 

Only  the  secular  component  of  Q  in  (100)  now  remains  and  it  follows  at 
once,  in  view  of  the  simplificctlons  indicated  by  the  last  sentence  of 
section  3.3,  that  the  ecpiation  can  at  last  be  rewritten,  using  (104),  as 


I  /  X  sin  Q\  +  Rj(-I)/x'\  '  (115) 

j  -  X  cos  a  1*1  I  I 

\  cY  -  sZ  /  \  0  /  \  Z’/ 


where  x',  Y*  and  Z*  are  non-singular  versions  of  X,  Y  and  Z  ,  given  by 

x'  -  r  cos  b  cos  w  -  r  sin  b  (b  cos  w  -  sQ.-)  “  rw'  cos  b  sin  w  ,  (ii6) 


Y  "  r  008  b  sin  w  -  rb  sin  b  sin  w  +  rw*  cos  b  cos  w 
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I  and 

4 

I  i’  -  r  sin  b  +  t  cos  b  <b  -  50^^  cos  w)  ;  {118) 

t  ^  here  w'  ,  «s  indicated  after  (114),  is  given  by 

! 

w'  -  a  +  +  cQ^p  .  (119) 

} 

i 

) 

% 

3.5  Loag-tena  pestusbatlona  and  singularity 

( 

;  We  )iave  seen  (low  to  get  from  t)ie  mean  elements,  C  t  at  time  t  ,  to  t)ie 

con^onents  of  position  and  velocity  in  the  normal  (equatorial)  axis  system, 
without  any  real  problem  associated  with  singularity.  It  remains  to  consider 
such  difficulties  as  exist  in  the  propagation  of  the  mean  elements  themselves, 

I  from  their  values  ((o)  at  epoch,  this  having  been  described  in  section  3.3  as 

'a  preliminary  step*.  The  schematic  propagation  formula  is  (57),  and  the  only 
'  difficulties  turn  out  to  be  with  the  long-period  terms,  .  One  difficulty 

relates  to  the  general  method  of  evaluating  each  A(  as  a  well-behaved  integral 

A 

of  C|p  t  *nd  the  others  relate  to  the  way  in  which  the  circular  and  equatorial 
singularities  are  dealt  with.  Wa  give  a  general  consideration  here  of  how  the 
difficulties  can  overcome,  reserving  the  full  details,  for  <13)  pertur¬ 

bations;  to  Part  3. 

(a)  General  method  of  evaluating  the 

As  remar)ced  in  section  3.1,  the  AC  arise,  by  definition,  from  components 
of  C  1**  which  (0  ,  but  not  v  ,  is  an  explicit  argument,  the  point  here  being 
that  U  has  a  first-order  secular  variation  due  to  J2  •  E'or  the  general  zonal 
harmonic,  J|  ,  there  are  trigonometric  terms  in  kco  for  all  values  of  k  ,  up 
to  <  -  2  ,  having  the  same  parity  as  <  ,  so  that  for  J3  we  just  have  (as  we 
shall  see  in  section  7)  terms  in  cos  <0  and  sin  10  .  if  the  elements  are  to  be 
propagated  over  long  periods  of  time,  it  is  obviously  desirable  that  the  J3 
i  ‘  analysis  should  reflect  the  secular  variation  of  CO  in  the  integration  of  these 

'  terms,  even  when  the  analysis  is  only  taken  to  the  first  order  in  J3  .  (The 

i 

I  ,  question  of  the  precise  meaning  of  the  terms  'first-order',  'second-order',  etc, 

!  was  considered  in  section  9  of  Ref  3  and  that  discussion  will  not  be  repeated  in 

detail  )tere.)  On  taking  J2  analysis  to  second  order,  we  have  similarly  (as  we 
shall  see  in  section  5)  terms  in  cos  2a>  and  sin  2®  ,  integration  of  which 
should  again  allow  for  the  first-order  variation  in  ®  .  On  this  basis  (and 
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with  J3  regarded  as  second-order  in  J2)  we  derive  the  complete  formal*  second- 
oro  c  solution  that  we  require;  it  will  be  found  that  we  need  to  incorporate  some 
terms  that  (as  Ref  3  argues)  may  formally  be  regarded  as  part  of  the  third-order 
variation. 

To  specify  the  problem  precisely,  we  write 

J  cos  ko)  dx  and  Sj^  /  sin  ka  dx  ,  (120) 

0  0 

where  a  ,  written  at  to  avoid  ambiguity,  is  given  by 


a. 


+ 


a  X 

sec 


(121) 


observing  that  a  is  not  the  same  as  a  (except  at  X  -  0)  ,  because  it  does 
not  include  the  long-period  term  Aa  ,  Then  the  problem  is  to  have  formulae  for 
evaluation  of  and  that  are  'universally  computable'  in  giving  accurate 

values  for  all  values  of  >  this  includes  (in  particular)  zero, 

being  zero  (as  far  as  the  first-order  J2  formulae  are  concerned)  at  the  so- 
called  critical  inclinations  given  (as  we  shall  see  from  (152))  by  g  >■  0  . 

The  formal  evaluation  of  the  definite  integral  C);  is,  of  course,  given  by 
(sin  kWt  -  sin  k(ao)/kUsec  >  but  this  reduces  to  the  indeterminate  0/0  ,  rather 
than  the  determinate  t  co.s  kWo  t  as  Msec  tends  to  zero.  To  avoid  this 
situation,  we  follow,  with  modifications,  the  course  taken  by  Merson*'’,  who 
introduced  the  bounded  functions  Fj,  F2  and  F3  ,  given  by 

Fj  (0)  •>  0*^  sin  0  ,  (122) 


F^  (0)  »  0'^  (1  -  cos  0)  (123) 

and 

F^  (0)  -  0'^  (0  -  sin  0)  ,  (124) 


*  The  errors  in  the  formally  con^lete  second-order  solution  can  only  be  regarded 
as  truly  o(J2^,  J3) ,  ie  0(J2^,  J2J3)  ,  if  the  timescale  is  sufficiently  short, 
le  t  must  be  o(J2*')  ;  generally  speaking,  a  formal  solution  of  j'th  order 
in  J2  will  produce  errors  that  are  0(J2J*i  t)  in  the  long-term  variation 
over  time  t  .  This  point  was  not  made  in  Ref  3, 
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Foe  values  of  0  that  are  'far  enough  from  zero'  these  functions  can  be 
computed,  without  loss  of  accuracy,  from  their  defining  expressions,  whilst  for 
values  near  zero  they  can  be  congsuted  from  a  small  number  of  terms  of  their 
power-serles  expansions,  which  are  then  rapidly  convergent,  the  limiting  values 
of  Fi(0),  F2  (8)  and  F3(0)  being  1,  ^2  and  \  respectively.  (Merson's 

functions  are  equivalent  to  three  of  the  functions  introduced  by  Stunpff^®.)  We 
also  define,  for  convenience,  V<o  given  by 

Vm  =  m  t  .  (125) 

sec 

Merson's  formula  for  C|t  was  derived  on  the  basis  that 
sin  )c5^  -  sin  )cMg  =  sin  +  Vco)  -  sin 

sin  kVo)  cos  -  (1  ~  cos  kVm)  sin  ,  (126) 

so  that  he  obtained  (with  different  notation) 


Cjj  “  {Fj(kV(0)  cos  kSg  -  kVt!)  Fj(kV(B)  sin  kS^^ }  t  .  (127) 


However,  it  is  more  natural  to  write 


sin  kCOj,  -  sin  kM^  ••  2  sin  )ik((0j,  -  cos  )ik(Mg  +  co^) 


2  sin  )4kV{o  cos  kti)„  , 

/at 


(128) 


since  we  then  have  the  siit^ler  expression 


-  {FjOikVa)  cos  hSjjj,  }  t 


(129) 


the  corresponding  expression  for  is 


5^  »  {FjOikVo)  sin  kSjj^. }  t 


(130) 
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Use  of  (129)  and  (130)  removes  all  difficulty  associated  with  the  general 
evaluation  of  the  At[  ,  except  that  for  each  of  AH  and  Ao)  it  Is  necessary  to 
allow  for  a  further  tern  that  is  inportant  though  formally  only  of  third  order^. 

A  i 

These  terms  are  induced  in  the  first-order  effects  of  Clsec  and  a>sec  by  the 
second-order  Ae  and  Ai  ,  since  daec  and  Uisec  are  functions  of  e  and  i  . 
The  coirputation  of  these  terms  involves  t!ie  integration  of  Ae  and  Ai  ,  so  we 
require  universally  conputable  formulae  for  the  (definite)  integrals  of  Ci,  and 
Sk  -  specifically  of  S2  for  analysis  and  Ci  for  J3  analysis. 

The  formal  expression  for  the  integral,  from  x  •>  0  to  x  -■  t  ,  cf  Cj, 

(as  a  function  of  X  )  is  -  (cos  “  cos  kfflo)  /  (k&sec)  ^  "  (sin  kWo  /kC0sec)t  • 

He  can  write  this  as 


2  sin  )ikV(o  sin  kS 


(kCD  ) 
sec 


,,  2  sin  HkV(i)  cos  ka)„^  -  sin  ko)., 

nt  «t  /it 

—  +  - r - t  , 

ku 


and  then  rearrange  the  terms  to  give 


2  sin  )4kV£0  cos  k(o 


k(3 


sin  kS  . 

. . .  (kco  t  -  2  sin  KikVco) 

•  2 
{k(a  ) 

soc 


From  this  we  obtain  at  once  the  required  result/  viz 


t 

I  Cj^  dx  “  )4  [Fj  ( iikVo))  cos  km^,j,  -  )4kVci)  F^OikVo))  sin  t^  ;  (131) 

0 


similarly 


t 

J  Sj^  dx  «  )4  [Fj  ( iiikVo))  sin  ^<0^^  +  >4kVo)  F^O^kVo))  cos  t^  .  (1?2) 

0 

When  the  induced  contributions  to  Aft  and  Ato  are  allowed  for,  the 
propagation  of  the  ^  should  (apart  from  the  question  of  singularity,  still  to 
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be  discussed)  give  accurate  results  over  quite  long  periods  of  tiiive.  As 
indicated  by  the  last  footnote,  however,  there  will  be  an  increasing  need,  in 
modelling  the  long-term  variation  as  the  number  of  revolutions  from  epoch  builds 
up,  for  terms  that  are  formally  of  the  third  order.  Without  going  to  these 
terms,  accuracy  can  in  principle  be  in^roved  if  a  long  interval  is  split  into 
short  intervals,  as  a  form  of  'rectification'  involving  subsidiary  epochs.  This 
brings  in  the  hybrid  <or  semi -analytical)  component  of  the  orbit  generator, 
referred  to  in  section  1  and  to  be  described  in  Part  3.  Ephemerides,  generated 
with  and  without  rectification,  will  be  compared  in  Part  3,  and  assessed  against 
ephemerides  produced  by  pure  numarical  integration;  this  assessment  proviaas  the 
ultimate  justification  for  all  the  formulae  developed  in  the  present  section  of 
the  present  Report. 

(b)  Singularity  avoidance  for  near-circular  orbits 

It  will  be  found  from  equations  (389)  and  (399)  that  0)3, fp  and  M3,£p  ,  as 
defined  by  equation  (72) ,  are  infinite  when  e  is  zero,  this  being  a  manifesta¬ 
tion  of  the  circular-orbit  singularity  that  has  been  referred  to  several  times. 
However,  the  phenomenon  is  merely  a  consequence  of  the  use  of  elements  that 
purport  to  define  the  position  of  perigee  accurately  (relative  to  the  ascending 
node,  in  regard  to  CO  ,  and  relative  to  the  satellite's  (mean)  position,  in  regard 
to  M  )  in  circumstances  where  the  concept  of  'perigee'  ceases  to  be  meaningful. 

If  the  elements  e,  (O  and  M  are  replaced  by  the  quantities  4  and  t|  ,  defined 
by  (9)  and  (10),  and  U  (  ■  M  +  <B ) ,  then  the  difficulty  vanishes. 

It  is  not  a  satisfactory  procedure  to  \ioi)r  throughout  witn  \  ,  i)  and  U 
as  elements,  however.  For  one  thing  the  analysis  becomes  much  more  complicated, 
particularly  (as  we  shall  see)  in  respect  of  U  .  For  another,  4  ant*  T) 
themselves  become  undefined  when  the  orbit  is  equatorial  but  not  circular,  but  we 
defer  all  questions  of  the  equatorial  singularities  until  sub-section  (c) . 

Finally,  it  is  essential  that  the  dominant  (first-order)  component  of  the  secular 
perturbation,  ^sec  ^  applied  to  <0  directly,  since  the  perturbation  occurs 

precisely  and  naturally  (and  hence  necessarily  non-singularly)  in  this  form. 
Fortunately,  the  effect  of  woricing  throughout  in  terms  of  £  ,  t)  and  U  can  be 
obtained  by  appeal  to  their  differential  relationships  to  e  ,  a>  and  M  , 
without  introducing  them  explicitly  at  all.  Thus,  since  e  Am  is  finite  for 
ali  e  ,  and  may  be  thought  of  as  a  single  entity  to  be  operated  with  instead  of 
Am  ,  wo  compute  e  and  m  (the  mean  values  of  e  and  m  at  time  t  )  via  e 
and  m  ,  where  m  was  introduced,  at  (121),  to  cover  secular  perturbations 
alone,  and  e  (defined  on  the  same  basis)  is  eqtial  to  io  ;  the  formulae 
requited  are 
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8n<i 


e  C03  ib  "  (e  +  Ae)  cos  O)  -  <i  Ao)  sin  (O 


(133) 


e  sin  (O’*  (e  ')■  Ae)  sin  (i)  -t-  (e  Aco)  cos  to  . 


(134) 


(Recovery  of  e  trom  e  cos  a>  and  e  sin  W  is  trivial;  recovery  of  <b 
involves,  in  Fortran,  the  use  of  the  ATAN2  function,  unless  e  is  zero,  in  which 
case  the  arbitrary  value  of  zero  may  be  adopted  for  . ) 

The  value  of  M  still  has  to  be  found,  but  this  is  easy,  if  the  non¬ 
singular  AU  is  )cnown,  since  it  is  only  necessary  to  match  M  to  Au  .  Thus 
the  formula  in  terms  of  AU  is 


M  =  (M  +  (0)  +  AU  -  (b  ,  (135) 

in  which  any  indeterminacy  in  w  ,  as  derived  from  (133)  and  (134),  is  properly 
reflected  in  the  preservation  of  a  fully  accurate  U  . 

The  fact  that  e  and  u  are  legitimately  derived  in  such  a  simple  manner 
will  be  considered  further  in  Part  3.  It  is  assumed,  of  course,  that  Ae  etc 
are  evaluated,  following  sub-section  (a),  by  integration  of  e^p  (etc);  e^p  is 
given,  via  equation  (72),  by  ea.Fp  and  63, |p  (assuming  only  J2  and  J3  to 
be  operative),  and  these  are  obtained  in  sections  5  and  7,  The  essential  point 
about  the  first-order  Taylor  expansions  (133)  and  (134)  is  that  they  effectively 
allow  for  first-order  A^  and  At)  faithfully,  without  contamination  by  second- 
degree  combinations  of  Ae  and  A(0  that  would  introduce  intolerable  error  when 
Aci)  is  unbounded;  a  geometrical  way  of  looking  at  this  is  that  the  formulae 
'straighten  out'  the  effect  of  an  e  Aci)  which,  near  the  origin  of  the  (4,  ■>!)■ 
plane,  has  a  pronounced  curvature. 

It  was  indicated  in  section  2,1  that,  for  analysis  in  which  there  is  no 
e-truncation,  U  is  not,  a  natural  parameter  to  use,  and  this  is  because  (as  we 
shall  see)  the  direct  addition  of  Am  and  A(0  does  not  lead  to  simple 
expressions.  Simple  (and  hence  natural)  expressions  are  given  by  AM  t  q  Ao)  , 
or  AU'  say,  and  on  this  basis  (135)  should  be  replaced  by 
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M  -  (H  +  q  U)  +  AO’  -  q  a  ,  (136) 

but  an  aw)cward  point  now  arisea,  viz  that  when  Aa  Is  large  (because  3  is 
small)  an  ambiguity  of  Zk  in  the  recovery  of  a  from  (133)  and  (134)  can  lead 
to  an  error  of  2K(1  -  q)  in  the  derivation  of  M  from  (136) .  There  is  a 
similar  difficulty  in  relation  to  the  inclination  singularity,  and  a  double 
difficulty  when  the  singularities  occur  together,  as  we  shall  see  in  sub-sections 
(c)  and  (d) .  It  is  enough  to  remar)c,  in  the  present  Report,  that  the  procedure 
proposed  in  Ref  4  for  resolving  these  ambiguities  was  inadequate,  and  new 
resolution  criteria  will  be  given  in  Part  3. 

(c)  Singularity  avoidance  for  near-equatorial  orbits 

The  situation  for  near-equatorial  orbits  closely  parallels  that  for  near¬ 
circular  orbits,  and  being  Infinite  when  s  is  zero.  It  is  now 

the  concept  of  'node'  that  ceases  to  be  meaningful,  and  with  direct  orbits  the 
difficulty  would  vanish  if  we  replaced  the  elements  i,  il  and  o>  by 
(paralleling  5,  1)  and  U  )  the  quantities  JJ  and  5  ,  equal  to 
sin  i  sin  tl  ,  -  sin  i  cos  and  <1)  +  O  ,  respectively.  The  use  of  to  has 

been  trad.itional  in  celestial  mechanics,  where  orbits  are  almost  invariably 
direct,  but  it  is  inappropriate  for  retrograde  orbits,  for  which  it  would  be 
necessary  to  redefine  3  as  to  -  fl  ;  furthermore,  as  with  the  addition  of  AM 
and  Ao)  to  form  AU  ,  the  perturbations  in  to  and  Q  do  not  naturally  either 
add  or  subtract  directly.  Both  difficulties  (different  regimes  and  unnatural¬ 
ness)  disappear  if  we  use  the  generalized  quantity  y,  defined  at  the  differ¬ 
ential  level  by  (31),  as  the  replacement  for  (O  . 

As  with  5,  1)  and  U'  ,  it  is  only  necessary  to  use  JJ  and  \)/ 

implicitly  in  the  relations  corresponding  to  (133),  (134)  and  (136),  viz 

3  sin  il  “  sin  (i  Ai)  sin  il  +  (s  Aft)  cos  ft  , 

^  —  a  ¥t  m 

-  3  cos  ft  "  -  sin  (1  +  Ai)  cos  ft  +  (s  Aft)  sin  ft 

and 

3  -  (5  +  eft)  +  Ay  -  eft  . 


(137) 

(130) 

(139) 
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As  with  e  Ao  ,  i  AQ  is  to  be  regarded  as  a  single  entity  in  (137)  and  (138) . 
Clearly,  s  and  fi  can  be  recovered  at  once  from  these  equations,  but  a 
difficulty  arises  in  the  recovery  of  i  from  s  ,  since  near-polar  orbits  have 
to  be  covered  as  well  as  near-equatorial  ones.  The  inaccuracy  in  then  deriving 
i  from  s  has  already  been  referred  to  (in  section  2.1)  and  an  extreme 
situation  could  arise  in  which  the  confuted  value  of  s  became  greater  than 
unity.  The  problem  is  easily  dealt  with,  however,  at  the  expense  of  a  minuscule 

■i 

amount  of  extra  computing,  by  malcing  use  of  c  ,  given  by  cos  (i  +  Ai)  ,  as  well 
as  s  ;  this  amounts  to  operating  with  all  three  conponents  of  the  unit  vector 
(^,  Vi  C)  •  also  Ref  7.) 

It  is  worth  remarking  that  the  quantities  $  and  rj  turn  out  to  be  'less 
natural'  than  the  corresponding  4  ^  »  associated  with  e  and  (0  ,  in  that 

the  perturbation  expressions  lack  the  homogeneity  (in  regard  to  occurrences  of 
the  factor  c  )  that  might  have  been  expected.  This  can  be  seen,  for  example, 
when  (149)  and  (151)  are  combined,  to  give  and  rji  and  (as  will  be 

remarked  again  in  section  4)  the  inhomogeneity  is  apparent.  Two  pairs  of 
quantities  exist  that  are  homogeneous,  however,  viz  tan  hi  cos  Ci  and 
tan  )il  sin  il  ,  the  pair  appropriate  for  direct  orbits,  and  cot  hi  cos  Cl  and 
cot  hi  sin  ,  the  pair  appropriate  for  retrograde  orbits.  These  quantities 

belong  to  the  two  sets  of  'equinoctial  elements'  introduced  by  Broucke  and 
Cefola^^,  one  sot  being  fully  non-singular  for  direct  orbits  and  the  other  set 
fully  non-singular  for  retrograde  orbits.  (The  remaining  equinoctial  elements 
are  a,  e  cos  U,  e  sin  W  and  M  -f  S  ,  m  being,  as  defined  earlier,  a  ±  Cl 
as  appropriate.)  An  obvious  advantage  in  redefining  {  and  rj  to  be  a  pair  of 
equinoctial  elements  would  be  that  no  (  would  then  be  required,  but  the  need 
for  two  different  pairs  is  not  very  satisfactory.  Further,  there  is  really  no 
disadvantage  in  the  'unnaturalness'  of  our  {  and  rj  ,  since  for  explicit 
formulae  (associated  with  J2  and  J3  )  we  are  sticking  to  Ai  and  5  ACl  ,  { 

and  T)  only  being  (implicitly)  introduced  in  the  general  transformations  (137) 
and  (138) .  He  do  require  explicit  formulae  for  Ay  ,  however,  to  avoid  singular 
expressions,  so  it  is  only  here  that  homogeneity  is  desirable;  it  seems  ironic, 
therefore,  that  equinoctial  elements  are  'natural'  where  it  does  not  matter,  but 
'unnatural'  (  (5  and  M  +  (i)  )  where  it  does  matter. 

There  is  a  more  basic  respect  in  which  the  parallel  between  the  eccen¬ 
tricity  and  inclination  singularities  breaks  down,  since  the  variation  of  O) 
constitutes  the  very  essence  of  the  A(  ,  by  definition,  the  variation  of  Q 
being  irrelevant.  Thus  the  study  of  i  and  u  (or  ^  and  ii  )  is  essen¬ 
tially  self-contained,  whereas  the  study  of  1  and  n  cannot  be  divorced  from 
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a  .  (Xn  another  sense,  of  course,  the  reverse  of  this  is  true,  as  1  and  Q 
define  the  orbital  clane,  and  its  node,  independently  of  O  ,  whereas  A)  only 
defines  the  perigee  relative  to  the  node.)  The  in^rtance  of  the  distinction 
should  be  clearer  when  these  studies  are  pursued  in  Part  3. 

(d)  Avoidance  of  all  singularity 

Putting  together  the  results  of  the  previous  sections,  we  can  now  summarize 
an  approach  that  in  principle  avoids  all  singularity.  The  conceptual  U'  ,  of 
sub-section  (b) ,  must  be  replaced  by  L  ,  where  £  is  given  by  (35),  and  our 
first  step  is  the  evaluation  of  the  five  specific  long-period  perturbations  de  , 
Ai  ,  0  AQ ,  e  Ay  and  Al<  ,  using  the  expressions  (of  type  42,fp  ) 

to  be  derived  in  sections  6  and  7;  the  evaluation  is  based  on  (129)  and  (130), 
with  (131)  and  (132)  also  invo)ced  for  the  induced  contributions  to  A(1  and 
A(4  ;  since  these  contributions  are  always  non-singular  (free  of  s  and  e 
divisors)  they  can  most  conveniently  be  incorporated  directly  with  the  purely 

*  B  at 

secular  generation  of  O  and  0)  (required  together  with  that  of  M  )  . 

The  values  of  i  and  fl  can  now  be  obtained,  in  principle  from  (137)  and 

(138)  but  with  c  used  as  well  as  s  ,  as  has  been  described.  We  next,  again 
in  principle,  derive  e  and  y  from  versions  of  (133)  and  (134)  in  which  O) 
has  been  replaced  by  y  ,  with  y  set  to  w  +  eft  ;  then  u  is  obtained  from 

(139)  ,  re-expressed  simply  as  u  >  y  -  eft  .  The  arbitrary  nature  of  y  imposes 
no  difficulty  in  the  use  of  cos  y  and  sin  y  ,  since  (133)  and  (134)  are 
preserved  under  the  addition  of  an  arbitrary  quantity  to  5  ,  and  in  fact  it  is 
convenient  to  use  versions  of  these  formulae  in  which  a>  is  replaced  by 

5  +  6  (ft  -  ft)  ,  instead  of  y  ,  since  w  i.<3  then  obtained  directly;  however, 
this  does  not  get  round  the  difficulty  associated  with  the  2x  ambiguity  in 
Cl  -  a  ,  and  we  return  to  this  in  Part  3. 

It  remains  to  derive  M  .  The  formula  for  this  is  a  modified  form  of 
(136),  viz 

M  -  (m  +  qy)  +  AL  -  qy  ,  (140) 

and  again  the  only  difficulty  is  the  one  associated  with  the  Zit  ambiguity  in 
V  -  V  • 
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4  WIMV-OKOOL  AKAXiTtZS  VOR  Jj  VSKXOABACZ01I8 

4 . 1  ia  tha  aaoaZatiag  alanaats 

Substitution  of  the  potential  given  by  (48)  into  the  planetary  equations, 
(22) -(27),  gives  exact  equations  for  the  variation  of  the  osculating  elements 
(excluding  M  );  expressed  with  the  help  of  tha  Cj  and  Sj  families  given  by 
(19),  these  equations  are 


Kna 

2q5 


{f  (SeSj  +  4S2  -  eSj)  +  4eh  sin  v}  , 


(141) 


•{f  (SeS^  +  I4S3  +  12eS2  +  2Sj^  -  eSg)  +  4h  (e  sin  2v  +  2  sin  v)  J  ^ 


(142) 


i  -  - 


7  If)  ' 


(143) 


n  -  #  (•?rc(c  -1)  , 


3  vr/  -'-2 


(144) 


(b  -  -S2-  /£\  {f  (5eC^  +  I4C3)  -  2e(4  -  7,;)  -  f  (2C^  -  eCp) 

8eq^  '  ^ 

+  4h(e  cos  2v  +  2  cos  v)  +  2e(6  -  7f)  } 


and 


(145) 


Kn  /£\^ 
8eq' 


+  14C3  -  18eCj  -  2Cj  +  eCjj) 

+  4h(e  cos  2v  +  2  cos  v  -  3e)  J  . 


(146) 


The  independent  variable  is  now  changed  from  t  to  v  (or  strictly  v  ) ,  via 
(39),  after  which  the  first-order  Integration  of  each  equation  is  straight¬ 
forward.  It  is  noted  that  v  (in  one  form  or  another)  acts  as  the  only  variable 
in  each  integrand,  being  Involved  in  three  different  ways:  (i)  explicitly; 
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iii)  via  the  Cj  and  Sj  (for  j  #  0  )  by  definition;  and  (iii)  via  p/r  and 
equation  (63) .  For  ease  of  expression,  ee  will  from  now  drop  bars  (and  tildes) 
in  V  ,  etc,  wherever  this  should  cause  no  confusion.  (To  first  order,  the 
results  of  this  section  are  correct  without  the  bars,  but  bars  must  be  addsd  to 
all  first-order  expressions  when  we  extend  to  second  order.) 

The  results  of  the  integration  are  as  follows,  when  expressed  in  terms  of 
the  notation  introduced  by  (71)  and  (73),  by  using  the  bracket  conventions  of 
section  2.1,  and  after  adoption  of  the  best  set  of  integration  constants  (as 
discussed  in  section  3.1  and  to  be  justified  in  section  4.3): 

Si  "  {3f(e^Cj  +  6e^C^  +  3eI4,l)Cj  +  4(2,3)02  *  3s[4,ljCj  +  +  e^C_j/ 

+  4h(e^  cos  3v  +  6o^  cos  2v  +  3e[4,l]  cos  v  +  2(2,3])  },  (147) 

‘'l  “  ’it  +  18eC^  +  (28,17)03  +  eOeO^  +  3(4,ll)Cj  +  1860^  +  3e^C_j) 


+  4h(e^  cos  3v  +  6e  cos  2v  +  3(4,1)  cos  v  +  lOe)  }  ,  (148) 

ij  -  -i-  sc(eC3  +  30^  +  3e0j  +  3)  ,  (3149) 

A 

Oj  -  "  O  ,  (150) 

flj  ”  ~  c  (®S3  +  383  +  3eS^  -  6e  sin  v)  ,  (151) 

-  2g  ,  (152) 


fflj  -  ~  e  "  {se^fSj  +  18efS^  +  [28f,  -  (8,-39)  ]  S3  -  126(2,-5)83 

-  3  [4f,  (8,  -15)]  Sj  -  18efS^j  -  3e^fS_j 
+  4e^h  sin  3v  +  24eh  sin  2v  +  6[8h,  (14,-17)]  sin  v  }  ,  (153) 
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Oj  -  qh 


(154) 


and 

Cj  -  «!  +  18*S^  +  (28,-13]S3  -  -  3[l,17)Sj 

-  ISeSjj  -  3e^S_j) 

+  4h(e‘  sin  3v  +  6«  sin  2v  +  3 [4, -5)  sin  v)  }  .  (155) 

A  certain  inconformity  will  be  observed  in  (153),  in  coiiq>arison  with  (148) 
and  (155)  in  particular,  in  that  we  do  not  have  the  same  neat  division  into 
'pure  f-tetms“,  involving  the  Cj  ,  and  'pure  h-terms',  involving  cos  jv  . 

This  inconfozmity,  propagated  from  (145),  follows  from  the  presence  of  the  term 
in  dc;/di  in  (26),  and  the  natural  partner  of  (146)  is  not  (153)  but  (166),  ie 
the  expression  for  yi  given  in  the  next  section,  since  (33),  the  equation  for 
y  ,  does  not  contain  du/ii  .  (For  further  comments  on  this,  see  Part  2.) 

It  only  remains  to  obtain  the  perturbation  in  M  (effectively  to  supersede 
the  perturbation  in  o  ) ,  the  basis  for  the  analysis  having  been  given  in 
section  2.2.  As  explained  in  section  3.1,  this  in  principle  involves  the 
derivation  of  a  modified  version  of  Kepler's  third  law,  but  it  turns  out  that  no 
modification  is  needed  in  practice;  this  is  fortuitous,  and  it  does  not  apply 
for  any  other  .  Thus  we  require  nothing  but  the  formula  for  itself. 

The  starting  point  is  (40),  and  to  use  this  we  require  a  suitably  expressed 
formula  for  nj  ,  the  derivation  of  which  calls  for  some  preliminary  discussion 
of  (147),  the  formula  for  ai  . 

As  explained  in  sections  2.3  and  3.1,  the  integration  constant  in  (147), 
part  of  which  appears  in  the  3f-component  of  the  equation  and  part  in  the 
4h-component,  is  such  as  to  make  the  mean  element  a  identical  with  the  exact 
constant  a'  given  by  the  conservation  of  energy.  But  a'  is  related  to 
(osculating)  a  by  the  exact  equation  (52)  which,  on  substitution  for  the 
potential  given  by  (43),  gives 


a 


1 


(3fC2  +  2h) 


(156) 
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This  Is  in  non-stsndsrd  for*,  bscsus*  of  th«  factor  ip/s)^  ,  but  expansion  of 
this,  hy  us*  of  (5),  at  one*  confirms  (147)  and  incidentally  explains  the 
manifest  symmetry  there  of  the  Cj  terms  -  the  3f-co*pon*nt  of  the  integration 
cosntant,  viz  6e^Co  ,  is  clearly  essential  for  this  symmetry,  to  match  the  term 
6e^C<  .  Of  the  two  (non-standard)  expressions,  intermediate  in  powers  of  p/r  , 
that  emerge  as  (156)  expands  to  (147),  the  first  is  the  one  corresponding  to  the 
required  formula  for  n^  ,  but  the  second  has  an  application  as  well,  to  be  mst 
in  section  5.9.  The  first  expression  is 


{3f  (eCj  +  2C2  +  eCj)  +  4h(l  +  e  cos 


and  the  second  is 


a^  -  —  aq‘2  {35(0^0^  +  +  2(2, HC^  +  4eCj  +  e^C^) 

+  4h  (e^  cos  2v  +  4e  cos  v  +  [2,1])  }  .  (158) 

To  each  of  the  four  expressions  for  ax  there  is  an  immediately 
corresponding  expression  for  nx  ,  given  by 

"1  -  -  '  <159) 


on  the  provisional  assun^tlon  that  the  n  we  require  is  indeed  such  that  the 
Kepler  law  is  satisfied  without  modification.  The  reason  why  (157)  leads  to  the 
most  apposite  expression  for  nx  is  that  the  factor  (p/r)^  then  cancels  out 
when  8n  (derived  from  this  expression)  is  substituted  in  (40),  W  in  (40) 
being  defined  by  (17) . 

Thus  we  have 

nj/nW  -  -  ■i-q{3f(eCj  +  20^  +  eCj)  +  4h(l  +  e  cos  v)  }  ,  <150) 


which,  in  view  of  (40),  is  in  the  right  form  for  further  integration.  The 
integration  leads  to  a  secular  term,  as  well  as  short-period  terms,  and  the 
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resulting  wcprcssicns,  with  th*  usual  notation  involving  the  synbol  introduced  in 
(3) ,  are  given  by 

A 

(161) 


and 

+  aSj  +  3eSj)  +  4h  e  sin  v}  . 


(162) 


The  secular  perturbation  in  M  is  now  given  by  the  combination  of  (154) 
with  (151),  and  the  short-period  perturbation  by  the  combinatioi.  of  (155)  with 
(162) .  The  former  combination  gives 

“  0  (163) 

as  anticipated,  validating  the  provisicnal  assumption  about  Kepler's  third  law 
for  mean  elements,  whilst  the  latter  combination  gives 

Mj  -  -  e‘^q{f  (le^Sj  +  18eS^  +  (28,-1)83  "  3t4,5]Sj  -  ISeS^^  -  3e^S_^) 

+  4h(e^  sin  3v  ¥  6e  sin  2v  +  3(4, -1]  sin  v)  ]•  .  (164) 


4 . 2  Perturbations  la  related  quantities 

For  reference,  we  give  formulae  for  perturbations  in  a  number  of  quantities 
related  to  the  standard  osculating  elements  covered  in  the  preceding  section;  we 
start  with  the  conceptual  yf,  p  and  L  ,  whose  rates  cf  change  are  given  by 
(33),  (34)  and  (35).  He  have 

¥1  -  h  ,  (165) 

Yj  -  e*^{f  (36^85  +  IQeS^  +  (28,11)83  *  ‘ 

-  18e6(,  -  3e^8_j) 

+  4h  (e^  sin  3v  +  6e  sin  2v  +  3(4,3]  sin  v)  }  ,  (166) 


% 
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(167) 


pj  *>  ■^q{£{«Sj  +  SSj  +  3«Sj)  +  4h«  sin  v}  ,  (168) 

“  hq  <169) 

and 

Lj  “  -J-  q{£(«Sj  +  3Sj  +  3eS^)  +  4he  sin  v}  .  (170) 

In  practice,  of  course,  the  yf  and  p  results  were  obtained  before  the  (D  and 
O  results,  (1S2)-(15S)  being  derived  from  (150)-(1S1)  and  (165)-(168).  Further, 

A 

it  will  be  noted  that  pi  ”  -  2ii  and  I>i  “  “  Ji  (and  li)cewise  Pi  “  -  zjj  and 
I'l  **  ~  Jl)>  these  are  not  coincidental  results,-  and  it  may  be  seen  from  Refs  1 
and  22  (see  else  Part  2)  that  there  is  a  general  flrst"Order  result,  for  the 
zonal  harmonic  ,  given  by  8p  --%(<+  1)  6j  and  5l  -  -  Vi(2l  -  1)  8j  . 

Next,  we  can  use  (147)  and  (146)  to  derive  pi  ,  since,  from  (28) 

Pi  “  ^  *1  ~  2ae  ej  .  (171) 

Thus 

Pj  -  ip{f{eC3  4  aCj  4  3eCj)  4  2h}  .  (172) 

We  can  erqiloy  the  last  two  equations  in  reverse,  to  obtain  a  more  corpact  (but 
non-standard)  expression  for  ex  based  on  the  most  contact  of  the  four  formulae 
for  a^  ,  viz  (156) .  This  expression,  which  can  be  confirmed  by  direct  reduction 
from  (148)  and  which  closely  corresponds  to  the  first-order  expression  for  5e 
given  by  Kozai^’,  is 

•i  ■  i  (l/  (3fC2  4  2h)  -  q^  [f  (eCj  4  3C2  4  SeCj)  4  2h]  }  . 

(Though  contact,  the  expression  is  unfortunate  in  having  the  factor  e"^  ;  as 
can  easily  be  seen,  however,  this  cancels  out  when  the  brac)cetod  terms  are 
expanded. ) 
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Again,  since 


(pc^)^  -  Pj  -  2p8C  i^  , 


(174) 


(149)  and  (172)  lead  to 


(pc^)j  -  ipc^  (1  -  3f)  . 


(175) 


This  is  a  (firso-order)  constant,  as  it  must  be,  since  (as  was  remarked  in 
section  2.2)  pc^  is  an  exact  constant  for  a  geopotential  that  involves  only  the 
zonal  harmonics;  it  represents  the  polar  component  of  the  satellite's  angular 
momentum.  The  reason  that  (pc^) i  does  not  actually  vanish  is  that  the 
constant.,  integration  Cor  a,  e  and  i  are  not  (optimally)  such  as  to  force 
this  to  happen. 

Perturbations  in  4  1)  ,  defined  by  (9)  and  (10) ,  can  be  obtained 

using  (148),  (132)  and  (153),  but  it  would  actually  be  unhelpful  to  use  (152)  to 
Introduce  quantities  f  *hd  ft  because  this  would  involve  the  conversion  of  a 
straightforward  (purely  secular)  &  into  an  apparently  long-periodic  variation 
(cf  remarks  in  section  3.5(b)).  However,  the  pure  short-period  perturbations  can 
usefully  be  expressed  with  the  help  of  the  Yj  >  dj  <  Cj  and  Sj  f  anti  lies, 
given  by  (18)  and  (20)  (and  so  far  unused);  in  terms  of  these,  (146)  and  (153) 
lead  to 

{  3e^fY5  +  18efY<  +  2  [l4f,  -{2,-9)]Y3  ••  12e(l,-5)Y2  -  3e^4,-13)Yi 
+  e^(9,-7)c3  +  36e(l,-l)C2  +  6  [2(4, -5),  (10,  -ID]  c^ 

+  2e(20, -21)C(,  -  9e^  (4,-5)  c_j  } 


and 


(176) 


Tlj  -  {ae^fOj  +  18efO^  +  2  [I4f,  -(2,-9)]<yj  -  i2e(l,-5)02  -  3e^(4,-13)0^  -  Se^fSj 
+  12e(l,-3)S2  +  6(2(4,-7),3(2,-3))Sj  +  2e(20,-39)Sp  -  3e^  (12, -13)  s.^  }  . 

. (177) 
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Foim-lae  for  and  JJi  ere  not  given,  because  (effectively  due  to  the 
c-factor  in  ii  )  ii  and  Qi  do  not  conbine  in  a  natural  (homogeneous)  way. 
This  was  remar)ced  in  section  3.5(c),  where  it  was  also  remar)ced  that,  for  the 
equinoctial  elements,  ij  and  do  combine  homogeneously.  The  resulting 

formulae,  for  both  pairs  of  equinoctial  elements,  are: 

(tan  )4i  cos  ft) ,  -  {3  cos  n  -  3e  cos  (v  +  ft)  +  3e  cos  (v  -  Q) 

1  6(1  +  cos  i) 

+  e  cos  (3v  +  O  +  20))  +  3  cos  (2v  +  £1  t-  20)) 

+  3e  cos  (V  +  ft  +  20))  }  ,  (178) 


(tan  )4i  sin  ft) ,  -  {3  sin  ft  -  3e  sin  (v  +  ft)  -  3e  sin  (v  -  ft) 

^  6(1  +  cos  i) 

+  e  sin  (3v  +  ft  +  2o))  +  3  sin  (2v  +  ft  +  20)) 

+  3e  sin  (v  +  ft  +  2o))  }  ,  (179) 


(cot  )4i  cos  ft)  ^ 


and 


. iM  {3  cos  ft  -  3e  cos  (v  -  ft)  +  3e  cos  (v  +  ft) 

6(1  -  cos  i) 

+  e  cos  (3v  -  ft  +  20))  +  3  cos  (2v  -  ft  +  2o)) 


+  3e  cos  (V  -  ft  +  2ci))  }  (180) 


(cot  )ii  sin  ft)  ^ 


^  ^ ■  {3  sin  ft  +  3e  sin  (v  -  ft)  +  3e  sin  (v  +  ft) 

6(1  -  cos  i) 

-  e  sin  (3v  -  ft  +  20))  -  3  sin  (2v  -  ft  +  2(0) 

-  3e  sin  (v  -  ft  +  2»)  }  .  (181) 


For  the  remaining  quantities  covered  in  this  section,  it  would  be  even  less 
helpful  to  consider  long-term  perturbatio.is  than  for  ^  ai)d  T)  ,  but  pure  short- 
period  perturbations  can  again  be  usefully  expressed.  Thus,  vi  can  be  derived 
from  (148)  and  (164),  using  (42);  the  result  is 
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Vj  -  -  •^e'Mf(3e^S5  +  18eS^  +  [28,11133  +  32eS2  -  [12,-5)S, 

-  iSeSg  -  3e^S.; ) 

+  4h{e^  sin  3v  +  6e  sin  2v  +  [12,1]  sin  v)  }  .  (182) 

From  vj  and  <0i  ,  che  lather  given  by  (153),  the  important  expression  for  uj 
can  be  obtained;  it  is 

Uj  -  -  -^{26(1,-1)33  +  (6, -7)82  +  2e(3,-5)S^  -  4e(5,-6)  sin  v}  .  (183) 

The  expressions  for  ^j,  Hi  and  ui  can  be  used  for  a  good  overall  check 
that  also  invoves  the  expression  for  ri  in  the  next  section.  Thus,  from  the 
formula 

e  cos  V  -  4  cos  u  +  T)  sin  u  (IS-ii 

it  follows  that 

(e  cos  V)  j  cos  u  •*  Hj  sin  u  -  0^(4  sin  11  -  T)  cos  u)  ,  (185) 

which  evaluates  to 

(e  cos  v).^  ’i  -^{£(36^0^  +  ISeCj  +  2[10,7)C2  -i  32eC^  + 

+  4h  (e^  cos  2v  +  8e  cos  v  +  [6,1))  }  .  0  86) 


The  factor  (1  +  e  cos  v)  can  be  extracted  from  the  right-hand  side  of  (186)  and 
the  resulting  expression  leads,  in  view  of  (5) ,  to 

(p/r)^  “  ^  (p/r)  {f  (3eC3  t  lOC^  +  lleC.  )  +  4h(e  cos  v  +  3)  }  .  (187) 

This  result  ties  up  Immediately  with  pi  ,  given  by  (172)  ,  and  ri  given  by 
(188)  . 

A  certain  property  may  be  observed  in  the  structure  of  the  various 
expressions  given  in  sections  4.1  and  4.2.  for  a  fully  non-singular  quantity, 
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C  t  the  property  for  is  that  the  coefficient  of  the  sine  or  cosine  of 
jv  +  ku  (  “  j'v  +  ka  ,  with  j’  “  j  +  k  )  contains  (at  least)  the  factors  e'i' 
and  s'*''  ;  if  it  were  not  so,  a  change  to  non-singular  variables  (involving 
e  cos  V,  e  sin  v,  s  cos  u  and  s  sin  u  )  would  result  in  uneliminated  singular 
terras,  via  a  negative  power  of  e  or  s  .  The  property  operates  in  the  first- 
order  perturbation  analysis  for  an  arbitrary  zonal  harmonic,  and  it  is  related  to 
the  so-called  d'Alembert  characteristic^^.  It  may  be  seen  to  apply  here  with 
»!»  Pi,  H  and  Pi  . 

The  property  does  not  apply  directly  to  ei  and  ii  ,  because  the  elements 
e  and  i  ,  though  non-singular  in  themselves,  do  not  have  well-behaved  rates  or 
change  at  their  associated  singularities.  (This  is  easily  seen  by  considering 
the  relation  ^  in  particular.)  Thus  it  is  eei  and  sii  ,  rather 

than  ei  and  ii  ,  that  have  the  property.  For  quantities  that  are  intrin¬ 
sically  singular,  the  effect  on  is  naturally  more  severe,  the  property  being 
satisfied  for  s^fl,  s^e^oji,  e^di,  e^Mi,  e^Vi,  s^^i,  s^Hi,  e^vi  and  s^ui  . 

in  the  next  section  it  will  be  found  that  ri,  bi  and  wi  also  have  the 
property,  this  being  an  essential  aspect  of  their  merit  of  'outse.  Finally,  for 
the  second-order  formulae  to  be  developed  in  sections  5  and  6,  where  muiif ication 
of  the  direct  property  relates  to  second-order  rates  of  change,  the  property  will 
be  found  to  hold,  in  particular,  for  a2,  e^e?,  3^i2,  P2,  s^n2»  s^e^052>  t?. 

b2  and  W2  . 

4 . 3  Perturbations  In  spherical  coordinates 

As  explained  in  section  3.3,  the  complicated  expressions  for  short-period 
perturbations  in  the  osculating  elements  (as  given  in  section  4.1)  can  be  com¬ 
bined  into  compact  and  convenient  expressions  for  perturbations  in  the  system  of 
spherical  coordinates  (r,  b,  w) .  Expressions  for  pure  short-period  perturbations 
can  be  obtained  from  the  six  basic  of  section  4.1  by  substitut.^ on  in 

(1)5) -(97).  The  results  are  worth  writing  with  unabbreviated  sines  and  cosines, 
in  view  of  their  importance,  and  in  terms  of  explicit  mean  and  semi-mean 
elements;  thus 

r.  "  4-  p(£  cos  2u  -  2h)  ,  (188) 

1  0 

Oj  “  esc  {sir.  (u  +  v)  -  3  sin  w} 

and 
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{f[3in  2u  +  4i  sin  (u  +  5»  ]  +  Seii  sin  v} 


(190) 


The  independent  chec)c  of  (188),  based  on  (187),  has  been  mentioned.  There 
is  no  obvious  chec)c  of  (189)  or  (190),  however  -  verification  that  (94)  is 
satisfied  proves  little,  since  (94)  was  used  directly  in  the  derivation  of  (97) . 

We  can  now  justify  the  claim  that  we  are  -jaing  the  best  set  of  integration 
constants*  for  the  .  With  any  other  choice  of  constant  for  aj  ,  (188)  would 
(via  (95))  contain  a  complicating  term  in  r/a  .  Kith  any  other  constants  for 
ei  and  Mi  ,  (188)  would  also  contain  a  term  combining  cos  v  and  sin  v  ,  and 
there  would  be  considerable  complication  in  (190) .  With  any  other  constants  for 
il  and  Hi  ,  (189)  would,  j.i)cewise,  contain  a  term  combining  cos  u  and 
sin  u  .  Finally,  any  other  constant  in  (o  would  be  transmitted  directly  into 
(190)  . 

On  differentiating  (188) -(190),  we  have  expressions  for  ri,  bi  and  wj 
immediately.  Thus 

tj,  =  -  p  f  u  sin  2u  ,  (191) 

bj  “  e  s  c  { (u  +  v)  cos  (u  +  v)  -  35  cos  5 }  (192) 

and 

w^  “  •|•{f[u  cos  2u  +  2e  (u  +  5)  cos  (u  +  5)  ]  +  4ehv  cos  v}  .  (193) 


For  a  purely  first-order  solution  w..  car.  set  5=0  and  v  =  u  =  Wn  in  these 
expressions.  For  incorporation  in  the  second-order  solution,  however,  we  must 
allow  for  the  first-order  variation  in  5  we  deal  with  u  just  by  identi¬ 
fying  it  with  w'  ,  as  given  by  (113)  -  cf  the  remac)i.  after  that  equation. 


*  Some  of  the  'constants'  are  represented  by  two  distinct  terms  in  the  Cl  in 
U47),  eg,  there  is  a  term  in  Cq  ,  wliich  contributes  a  first-order  constant  to 
Kai  ,  whilst  the  final  term  may  be  regarded  as  a  'pure  constant',  ma)cing  a  con¬ 
tribution  that  is  constant  to  second  order.  In  the  analysis,  similarly, 

Ca  may  contain  three  distinct  'constant'  terms,  of  which  the  first  two  are 
'less  constant'  than  the  third. 
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Equations  (188) -(190)  lead  at  once  to  Sc ,  Sb  and  5w  ,  since  there  are 
"o  Ci.tp  in  (74) .  If  we  want  the  con^jlete  short-period  perturbations  (5r,  Sb 

A  ^ 

and  Sw  ) ,  we  must  allow  for  the  carry-over  effects  of  O-i  and  (0i  .  The 
resulting  formulae  are  given  by 


5r  ” 

Kr^  , 

(194) 

5b  ” 

K  (b^  +  cam  cos  u  )  ^ 

(195) 

5w  “ 

K  (wj^  +  iim)  . 

(196) 

The  full  expressions  can  easily  be  shown  to  be  equivalent  to  the  first-order 
terms  of  (366) -(368)  of  Ref  3,  so  long  as  an  error  in  (366)  ibid  is  corrected  by 
changing  ...Ka{(fq^  cos  2u  . . .  into  ...Kp{(f  cos  2u  . . .  .  A  second  error  in 
(366)  ibid  is  worth  remarking,  viz  that  ...  +  4(4  -  12f  -  19f^)  ...  should 
read  . . .-  +  (16  -  48f  +  23f^)  ...  ;  this  pair  of  errors  (in  the  same  equation) 
are,  apart  from  rather  obvious  ones  in  equations  (35)  and  (81) ,  the  only  errors 
known  in  Ref  3.  (The  quoted  equations  from  Ref  3  are  from  the  'Note  added  in 
proof;  this  was  written  after  the  author  suddenly  registered  the  superiority  of 
a'  as  mean  semi -major  axis,  having  started  by  recommending  an  a  only  ©(Ke^) 
different  from  Kozai's  a  .)  From  (194)-(196)  we  can  get  the  corresponding 
expressions  for  8r,  5b  and  6w  .  For  the  reason  indicated  in  section  3.2, 
however,  the  first-order  solution  with  carry-over  effects  included  is  not  a  good 
basis  for  development  of  the  complete  second-order  solution;  so  it  is  best  to 
regard  (194) -(196)  as  of  only  academic  interest  -  the  important  equations  in  this 
section  are  (188)-(190)  and  (191)-(193). 

5  J2^  PERTURBATIONS  IN  OSCULATING  ELEMENTS 

Before  proceeding  to  the  second-order  solution  for  the  J2-only  field,  it 
is  worth  stressing  again  that  the  first-order  solution  given  in  section  4  should 
in  general  be  understood  to  be  oxpressed  in  terms  of  semi-mean  quantities  and  not 
osculating  ones.  Thus  the  arguments  of  the  trigonometric  functions  should  be 
interpreted  as  if  tildes  were  carried,  and  the  coefficients  of  these  functions  as 
if  bars  were  carried,  there  being  no  distinction,  for  a,  e  and  i  ,  between 
mean  and  semi -mean  elements.  The  omission  of  bars  and  tildes  was  a  matter  of 
vi  lal  convenience,  and  was  in  any  case  legitimate  when  second-order  terms  were 
not  under  consideration. 
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Me  can  susnarlze  the  first-order  solution  for  the  generic  osculating 
element  as  follows.  For  Q  and  to  only,  there  is  a  secular  perturbation  given 
by  K^^nt  (cf(71)),  that  propagates  the  mean  element  C  >  together  with  a  carry¬ 
over  effect,  in  ,  that  (added  to  C  )  yields  the  semi-mean  element  C  • 

There  are  no  long-period  perturbations,  but  for  each  element  there  is  a  pure 
short-period  perturbation  given  by  (cf  (73)). 

The  second-order  solution  for  con^rises  the  first-order  terms  just 
summarized  and  the  second-order  terms  now  to  be  developed.  The  understanding 
about  the  omitted  bars  and  tildes  in  the  first-order  terms  now  becomes  vital,  but 
(because  we  are  not  concerned  with  third-order  terms)  it  will  be  legitimate  (and 
certainly  easier  on  the  eye)  if  bars  and  tildes  are  omitted  in  the  new  terms 
wherever  possible.  In  addition  to  the  second-order  secular  and  short-period 

A  A 

perturbations  specified  by  and  respectively,  long-period  perturbations 

will  now  be  majclng  their  appearance.  Once  the  expres.sions  for  the  {p 
available  (cf  (72))  -  and  they  are  summarized  in  section  6.1  -  the  J2^-component 
of  each  is  effectively  )cnown,  since  the  method  of  generating  it  was  given  in 
section  3.5(a).  Further,  it  will  be  found  that  these  components  of  the 
(unlike  those  associated  with  J3  )  are  all  non-singular.  The  carry-over  effects 
of  the  roust  not  be  forgotten;  they  are  incorporated  in  the  5^  (cf  (74)) 

and  the  amalgamation  of  the  into  Sr,  Sb  and  Su  is  the  subject  matter  of 

sections  6.2  to  6.4. 

A 

The  basic  idea  in  deriving  expressions  for  the  and  Cj  is  to 

•bootstrap'  on  the  first-order  solution  by  substituting  it  in  the  right-hand 
sides  of  the  planetary  equations  and  then  re-integrating.  A  special  procedure, 
free  of  integration,  is  possible  for  C  -  a  ,  and  this  is  developed  in 
section  5.1,  paralleling  the  derivation  of  (147)  via  (156).  We  then  proceed,  in 
section  5.2,  to  an  effective  rederivation  of  a2  (33  and  32, 4p  do  net  arise) 
by  the  general  procedure,  and  extend  it  to  the  other  elements  in  sections  5.3  to 
5.10.  It  will  be  found,  in  sections  5.3  and  5.10,  that  n  can  no  longer  most 
conveniently  be  identified  with  n*  . 

5 . 1  Perturbation  in  a  (spacial  method) 

We  want  82  such  that 

a  -  a  +  -i-  Ka  q'^(p/E)^  (3f  62  +  2h)  +  K^a^  ,  (197) 


the  starting  point  for  which  is  the  exact  equation  (52),  which  leads  (still 
exactly)  to  (cf  (156)) 
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a  -  a*  +  K  a'q"^  (p/r)^  (3f  Cj  +  2h)  .  (198) 

On  conparing  (197)  and  (198)  it  follows  that  there  will  be  five  sources  of  terms 
contributing  to  a2  ;  four  of  these  sources  are  essentially  the  same  as  four  of 
the  five  listed  in  Refs  2  and  3,  but  source  (ii)  of  these  no  longer  arises  now 
that  we  have  unreservedly  identified  a  with  a'  ,  and  it  is  replaced  by  a 
source  that  we  could  previously  neglect  (due  to  e-truncation) .  The  sources  are 
as  follows:  (i)  the  variation  of  p  in  K  ,  which  could  easily  be  overlooked; 

(ii)  the  variation  of  e  in  q"^  (this  is  the  new  source) ;  (iii)  the  variation 
of  (p/r)^  ;  (iv)  the  variation  of  i  ,  which  affects  f  and  hence  h  ;  and 
(v)  the  variation  of  u  in  C2  .  The  contributions  of  these  sources  are  as 
follows,  with  bars  and  tildes  suppressed  in  line  with  the  general  principle  laid 
down: 


(i)  -7  q"*  (p/r)^  OfCj  +  2h)pj  ; 

(ii)  laeq'"  (p/r)^  i3£C^  +  2h)e^  ; 

(iii)  aq‘^  (p/r>^  (3fCj  +  2h)  (p/r)^  ; 

(iv)  2aq'^  (p/r)^sc(C2  -  Di^  ; 

(V)  -2aq‘'^  (p/r)^  fS^Uj  . 


Here  pi,  ej  (p/r)i,  ij  and  uj  are  given  by  (172),  (148),  (187),  (149)  and 
(183)  respectively.  In  using  these  quantities,  we  are  effectively  identifying 
5^  with  6p^  (as  appropriate),  this  being  legitimate  as  82  is  to  be 
multiplied  by  ;  the  appearance  of  £2  ,  as  opposed  to  £2  ,  in  (197) 
validates  the  use  of  ui  (associated  with  fii  and  ^u  rather  than  8u  )  . 

It  is  very  convenient  that  (187)  expresses  (i/r)i  in  a  form  containing 
(p/r)  as  a  factor,  since  this  ensures  that  each  source  of  a2  contains  (p/r)^ 
as  a  factor.  On  multiplying  out  the  remaining  factors  in  each  source,  and 
combining  the  resulting  expressions,  we  eventually  derive  the  required  result, 
viz  that  (with  bars  and  tildes  omitted) 
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(1+0  co.-i  v)^  [  3f^  {  3e¥^  +  ISe^g  +  15et2,lirg  +  4[4,ll)r^ 

•  e[46,-l]r3  +  ISe^r^  +  3eYJ 

+  12f  {ae^hCg  +  12e^hC^  +  e  [2  (22,-31),  -(14,-17)]C3 
+  8h(3,2)C2  -  e[2{2,3),  -(34,-39)]Cj 

+  12e^hCg  +  2e^hC_j  } 

{22  2 

e  (8h  +  9f  )  (e  cos  3v  +  6  cos  2v) 

+  3e[2(24,-56,49),  -{8, 8, -37)]  cos  v 

+  2  [2  (20h^  +  3f^  +  39e^f^  ]  }  ]  ; 

. (199) 


the  brac)cet  conventions  of  section  2.1  ate  used  in  this  result. 

5 . 2  Perturbation  In  a  (general  method) 

This  section  serves  as  a  prototype  for  the  application  of  the  general 
'bootstrapping'  procedure  to  one  orbital  element  after  another.  The  fact  that  it 
permits  a  rederivation  of  (199)  acts  as  a  check  on  the  correctness  of  the  general 
method,  and  also  on  the  algebraic  reductions  involved  (for  a  )  in  both  methods. 

The  starting  point  is  the  exact  equation  (141),  converted  into  an  equation 
for  da/dv  by  use  of  (39) .  To  have  a  valid  basis  for  second-order  analysis  it 
is  essential  that  we  distinguish  carefully  between  barred  and  unbarred  quantities 
in  this  conversion,  the  required  form  of  the  equation  being 


~  -  )4Ka  (n/n)  (q^/q^)  {(p/r)'*/ (p/r)^}  Q  ,  (200) 

dv  ® 


where 


Qa  “  f(5eS3  +  4S2  -  eSi)  +  4eh  sin  v 


(201) 


To  express  da/dv  as  a  function  of  mean  c[uantities  only,  we  must  eliminate 
the  osculating  quantities  a  (and  the  basically  equivalent  n  ),  i,  e,  v  and 
(0  ,  and  we  cio  this  by  substituting  the  first-order  solutions  for  these  five 
quantities  into  the  right-hand  side  of  (200),  via  a  first-order  Taylor  expansion 
relative  to  the  corresponding  (semi-) mean  quantities.  A  special  situation  arises 
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in  regard  to  a  (iopliclt  in  via  S3/  S2  and  Sx  ,  and  explicit  when  each 

Sj  is  replaced  by  sin<jv  2(D)),  and  this  requires  prelindnary  consideration. 

First-order  analysis  effectively  involved  the  identification  of  barred  and 
unbarred  quantities  in  (200)  and  (if  we  ignore  for  the  moment  the  coiqilication 
caused  by  the  resulting  factor  (p/r)^  )  the  interpretation  of  each  Sj  in  (201) 
S3  (zero-order)  sin(jv  +  2(6)  ,  with  a  taken  as  constant  in  the  integration 
with  respect  to  v  ;  the  result  of  the  integration  was  then  naturally  taken  as 
-{l/j>  cos(jv  +  2u)  ,  with  CD  evaluated  (after  the  integration)  to  allow  for 
.  In  the  second-order  analysis  we  require  a  replacement  for  Sj  that  is 
first-order  correct,  ie  given  by  sin{j(v  +  Kvi)  +  2  (a»  +  Kcoi)  }  ,  but  in  the 
first-order  con^onent  of  the  extended  solution  we  will  still  not  be  making  proper 
allowance  for  the  variation  of  a>  ,  since  we  want  to  go  on  taking  the  integral 
in  this  con^onent  to  be  -  Cj/j  ;  this  means  that  there  is  a  source  of  second- 
order  perturbation  (for  each  Sj  in  Qa  )  additional  to  the  direct  Taylor 
effects  associated  with  vi  and  (Di  (viz  given  by  jic  Cj  vx  and  2K  Cj  (Di 
respectively) .  It  is  this  additional  source  of  perturbation  that  constitutes  the 
'special  situation*  we  are  considering,  being  associated  with  the  error  in  taking 
the  integral  of  Sj  to  l>e  -  Cj/j  when  (to  second  order)  it  is  really 
-Cj/(j  +  2KMi  )  .  But 


1 

j  +  2KU. 
*'  1 


1 

T 


0{K^ 


(202) 


so  the  error  can  be  rectified  by  introducing  the  'carry-over'  term  (2K£)i/j^)  Cj  . 
Since  this  term  is  the  integral  (with  respect  of  v  )  of  (-2K(0i/j)  Sj 
(neglecting  terms  that  would  carry  over  to  third  order!),  we  can  cope  with  the 
'special  situation'  by  supposing  that  there  is  an  'd)  contribution'  to  da/dv  , 
as  well  as  the  obvious  (O  contribution.  This  d)  contribution  is  such  that  each 
Sj  in  Qa  would  (apart  from  the  complication  of  the  (p/r)^  factor  that  has 
been  Ignored  so  far)  make  a  contribution  to  da/dv  proportional  to  (-2Kmi/j)  Sj 
as  well  as  the  contribution  proportional  to  2K  Cj  (0i  . 

At  first  sight  we  can  only  deal  with  the  complication  of  the  (p/r)^  factor 
by  multiplying  each  Sj  by  this  factor,  interpreted  as  (1  +  e  cos  v)^  ,  to 
yield  a  set  of  terms  in  Sj-2,  Sj-i,  Sj,  Sj,i  and  Sj+2  ,  and  applying  the  general 
procedure  to  the  individual  terms.  This  is  exactly  what  we  have  to  do  with  the 
V  and  (0  contributions  to  da/dv  ,  but  for  the  d)  contribution,  which  we  are 
currently  concerned  with,  we  can  escape  the  complication  by  making  use  of  a  known 
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expression  for  Si  directly.  The  argtment  of  the  last  paragraph  effectively 
shows  (on  nailtiplying  through  by  -j  )  that  to  each  term  of  aj  in  the  form  of  a 
multiple  of  a  particular  Cj  there  corresponds  a  carry-over  contribution  to 
da/dv  consisting  of  the  same  multiple  of  2K^  (bi  Sj  .  (The  same  principle  applies 
to  the  analysis  for  every  element,  ,  of  course;  when  a  multiple  of  Sj  ,  as 
opposed  to  Cj  ,  is  involved  in  ,  arising  from  the  appearance  of  S^/j  in  the 
first-order  Integration  of  ,  then  the  second-order  carry-over  to  d^/dv  is 

that  same  multiple  of  -2K^  tbi  Cj  .)  The  important  point  is  that  the  factor 
2K^  (b,  is  independent  of  j  ,  and  this  means  that  the  argument  is  not  affected  by 
an  unexpanded  power  of  (p/r)  in  aj  ;  indeed,  we  can  take  the  d>  contrioution 
(to  da/dv  )  from  the  sinplest  possible  expression  for  ax  ,  viz  (156)  .  (The 
general  principle,  being  independent  of  j  ,  applies  even  when  j  =  0  ' ) 

Mith  a  self-explanatory  notation,  we  may  now  list  the  six  contributions  to 
the  second-order  component  of  da/dv  as  follows,  dropping  all  bars  and  tildes, 
even  from  v  (it  is  remarked  again  that  the  Da  contribution  covers  the 
variation  in  n  as  well  as  a  ) ; 

DJda/dv)  -  ,  (203) 

2  2 

Dj  (da/dv)  -  -  K  aq  (p/r)  sc(5eSj  +  -  eSj  -  6e  sin  v)  i^  ,  (204) 

D^  (da/dv)  -  -  )iK^aq'^  (p/r)  { [  9eq'^  (p/r)  +  4  cos  v]  +  (p/r)Q^^}ej  ,  (205) 

where  =  5fSj  -  fS^  +  4h  sin  v  ,  (206) 

D^(da/dv)  -  )4K^aq'^  (p/r)  {  4eQ^  sin  v  -  (p/r)Q^^}vj  ,  (207) 

where  =  f (15eCj  +  SC^  -  eCj)  +  4eh  cos  v  ,  (208) 

D^(da/dv)  -  -  K^aq'^  (p/t)^  f  (5eCj  +  4C2  -  eCj)a)j  (209) 

and 

D^ida/dv)  -  4K^aq’^  (p/r)^  fg  .  (2io) 
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Th«  six  contributions  are  most  conveniently  dealt  with  by  grouping  into  the 
four  terms  that  contain  Qa  *  common  factor  <the  Da-term,  two  De-terms  and 
one  Dy-tezm)  and  the  five  other  terms  (one  associated  with  each  of  D^,  De,  Dy, 

D|a  snd  Dm  ).  Seven  of  the  nine  terms  contain  (p/r)^  (at  least)  as  an 
explicit  factor,  the  reason  that  only  the  first  power  of  (p/r)  appears  in  the 
other  two  terms  being  that  these  are  the  ei  and  vi  conponents  of  what  is 
really  a  term  in  <p/r)  i  ;  but  (p/r)  i  itself  contains  p/r  as  a  factor,  joeing 
given  by  (187);  thus  (p/r)^  is  an  overall  factor. 

The  combination  of  the  two  grcups  of  terms  may  L  .*  denoted  by  D2  (da/dv)  , 
and  we  find  after  some  tedious  algebra  that 

Dj  (da/dv)  -  -  .^K^aq'^  (p/r)^  x 

X  [f^  {iSe^^g  +  102e^£.,  +  66^(38,11)1:^  +  2a(103, 128)  +  4  (16,  86, 3) 

+  2e(73,32)S3  +  ee^te,!)^^  -  6e\  -  3e\  } 

+  8f{4e‘*hS,  +  26e^hS^  +  e^ (3 (30, -43) ,  4  (5, -6) IS. 

0  3  4 

+  e[6(16, -23),  (2, -9)) S3 

+  (24h,  6(2, -5),  (34,-39)1  52  -e(2(4,-J),  -(22,-27))Sj 

+  e^((2,3),  -2(8,-9))So  “  -  2e'’hS.2  } 

+  2  {  e^  (8h^  +  9f^)  (e  sin  4v  +  6  sin  3v; 

+  2e^ (8 (8, -20, 19),  -(8,8,37)  1  sin  2v 

+  2e((64,-176,14M,  -4  (2, -2, -11) )  sin  v  }  ]  . 

. (211) 


It  is  not  obvious  how  we  can  integrate  (211)  without  first  expanding  the 
factor  (p/r) 2  ,  but  we  are  in  the  happy  position  of  knowing  what  we  expect  the 
answer  to  be.  Rather  than  attempt  to  integrate  (211),  therefore,  we  differenti¬ 
ate  (199)  with  respect  to  v  ,  using  the  fact  that,  eg. 
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<!  +  •  cos  v)^Cj}  -  -  54(1  +e  cos  v)^  {«(j  +  3)S^^j  +  2j  +  e(j  -  3)S^_j  }  . 

. . <212) 

We  find,  after  soine  further  tedious  algebra,  that  da2/dv  ,  derived  in  this 
way,  is  identical  with  D2(da/dv)  ,  as  given  by  (211).  This  completes  the  chec): 
that  the  special  and  general  methods  yield  the  same  result. 

In  sections  5.9  and  6.2  we  shall  need  the  expression  for  a2  with  (p/r)^, 
rather  than  <p/r)^  ,  as  a  factor.  This  can  be  obtained  at  once  from  (199),  of 
course,  and  will  be  given  here  for  con^leteness .  In  section  5 . 5  we  also  need  the 
expre.ssion  with  no  p/r  factor  at  all,  but  this  is  so  lengthy  that  it  seems  best 
to  suppress  it.  (He  also  suppress  the  version  of  a2  that  has  a  single  factor 
p/r  ;  finally,  it  is  easy  to  show  that  a  further  such  factor  cannot  be  removed 
from  (199),  which  is  on  this  basis  the  simplest  possible  expression  for  .) 

The  expression  required  in  sections  5.9  and  6.2  is: 

"  28?  (1  +  e  v)^  X 

X  [sf^fse^Fg  +  24e^r.j  +  6e^(ll,3jrg  +  4e(19,23)rj  +  2(16,82,7)r^ 

+  12et9,  5)  +  2e^(41,l)r^  +  24e^r,  +  3e  } 

+  24f{e^hCg  +  8e^hC5  +  e^t(34,-49),-(6,-7))C^  +  4e[2  (7,  -10) ,  -  (0, 1 ) ) 

+  C24h,  2(18,-29),  (10,-11)JC2  *  4e  [2  (1, -3 ) ,  3  (4,  -  5) ) 

+  e^((10,-21),3(6,-7)  JCg  f  8e\c_j  +  e‘’hC_2  } 

+  2  {  e^  (8h^  +  9f^)  (e  cos  4v  +  8  cos  3v) 

+  2e^[3(40,-104,103),-(8,24,-69))  cos  2v 
+  4e[(112, -288,243), -3(0, 16, -45))  cos  / 

+  (8(20h^  +  3fS, 6(24, -56, 75), -3(8, 8, -37)}]  _ 

. (213) 
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5 . 3  fUEturbstion  1a  t 

Ths  stacting  point  is  given  by  converting  the  exact  ecpiation  <142)  into  an 
equation  £or  de/dv  ,  making  use  ot'  (33) .  The  result  is 

-  -  -i-KCn/n)  (q^/q^)  {  (p/r)  (p/r)  ^}q  ,  (214) 

dv  ”  ® 

where  Q-  “  f  (SeS  +  14S  +  12eS  +  2S,  -  cS  )  ■••  4h(e  sin  2v  +  2  sin  v)  .  (215) 

We  proceed  as  in  section  5.2,  listing  six  contributions  to  the  second-order 
corrponent  of  de/dv  .  With  self-evident  notation  they  are: 

D^(de/dv)  -  xV^  (p/r)Qgaj  ,  j2i(;) 

2 

(de/dv)  -  --^K  (p/r)sc(5eS^  +  IsS^  +  12eS2  +  28^  -  eS^^ 

-  5e  sin  2v  -  12  sin  v)i,  ,  (217) 

(de/dv)  -  -•|•K^{Qg[7eq■^(p/r)  +  3  cos  v]  +  (p/r)Qce}e.  ,  (213) 

-  f(5S^  +  12S2  -  Sg)  +  4h  sin  2v  ,  (219) 

D^(de/dv)  -  sin  v  -  (p/r)Oev)v,  ,  <220) 

where  Qev  "  2f(10eC^  +  21Cj  +  12eC2  +  Cj)  8h(e  cos  2v  +  cos  v)  ,  (221) 

2 

D^(de/dv)  "  -•^K  (p/r>f(5eC^  +  Isc^  +  1260^  +  2C^  -  oCjj)(0j  (222) 

and 

D^(de/dv)  -  -^K^fg  {30^85  * 

+  3(4  +  lle^)Sj  +  10e8j,  +  3e^S_^  }  .  <223) 
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As  in  section  5.2,  the  six  contributions  may  be  grouped  into  four  terms  with 
Qe  as  a  conmon  factor  and  five  other  terms.  The  first  two  Qg  terms  (and  only 
these)  are  found  (after  substitution  for  aj  and  ei  )  to  involve  the  factor  , 

but  their  combination  yields  a  cancelling  factor  q^  ,  so  q"^  does  not  appear  in  the 
overall  combination.  Six  of  the  nine  terms  contain  the  factor  p/r  and  it  can  be 
extracted,  as  in  section  5.2,  from  a  combination  of  two  of  the  others,  but  it  does 
not  occur  in  the  06  term,  so  p/r  cannot  !>e  factored  out  of  the  overall 
combination.  Some  lengthy  algebra  eventually  leads  to: 

1  2-1 

D,  (de/dv)  «  — —  K  e  x 
2  1536 

X  [f^{l5e‘'Ejg  +  162e^£g  *  8e^(82,13)£g  +  14e[84,53)£^  +  [78  4,  1712,3911^ 

+  20e(54,l]E  -  8(56,58,53]£,  -  12e[66, 103)£, 

D  H  3 

+  (48,-624,-371)£^  +  2e[36,-91l£j  -  18e^£.i  -  3e'’‘£.2  } 

+  8f  {4eSSg  +  42e^hS,  +  e^ [172h, -(110, -14 1))  Sg  +  5e  [64h,  -  3  (38, -49) )  Sj 

+  [224h  -8(94, -123), -(338, -459)  IS, 

-  2eC4(14,-23), 5  (78, -113)153  -  2[80h,-2{26,  9),  3  (10,  -23))S^ 

-  6eri6{3,-4),-(98,-103)lSj  -  8e^q^ [14, -151 

-  e[48h,  (158,-165)  lS_j  -  [52h,  3  (18, -19)  1 

-  ISe^hS.,  -  2e'’hS.4  } 


+  2  {  e^(8h^  4  9f^)  (e  sin  6v  +  i0  sin  5v) 

+  4e^[2{40h^  +  37f‘) ,  -  (24, -8, -31)  )  sin  4v 
+  2er36(8h^  +  5f^),  -(296,-312,-165)  1  sin  3v 
+  [16(24h^  +  7f^),  -16  (72,-136,33), -(216,120,-551)  1  sin  2v 
-  48[2(88h^  -  21f^),  (168,-08.-207)  1  sin  v  }  ]  .  (224) 
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The  term  in  Sq  is  independent  of 
long-period  variation;  the  absence  of  a 
noted,  and  this  applies  to  every  element 
write,  in  the  notation  of  section  3.2, 


V  and  represents  our  first  source  of 
corresponding  term  in  £q  will  be 
except  a  (see  section  5.8).  We  may 


e 


2.fp 


--^eq^f(14  -  15f)SQ  , 


(225) 


the  contribution  to  e^p  being  given  (c£  (72))  by  n  62. Bp  •  This  leaves  the 
short-period  carry-over  term  given  (cf  (74))  by  e2,fp  •  The  method  of 
integrating  62.tp  with  respect  to  time,  to  give  the  J2^  component  of  Ae  ,  was 
discussed  in  section  3.5;  see  also  Part  3. 


The  integration  of  the  rest  of  (224)  is  immediate,  except  that  we  have  tc 
incorporate  the  right  'arbitrary  constant'.  This  is  naturally  expressed  via 
terms  in  Fq  and  Cq  ,  together  with  a  'pure  constant'  (see  the  footnote  in 
section  4.3)  that  is  a  function  of  e  and  i  alone,  the  'coefficients'  being 
chosen  to  suit  the  derivation  of  r^  in  section  6.2.  Anticipating  the  rj 
analysis,  we  derive  (from  the  integration)  the  required  result  for  <32  ,  viz 


1 

9216 


e 


-1 


X 


108e^r„ 


+  6e^t82,13)rg  +  12e[84,53)r^  +  [784,  1712, 39irg 


+  24e[54,l]rj  -  12t56,58,53)r^  -  24e [66, 103) 

+  3[48,-624,-3Tljrj,  +  12e [36, -91 J +  6e^  [66, -35] 

+  108e^F_j  +  9e‘*F_2} 


+  4f{6e\Co  +  72e\c,  +  2e^ (172h,  -  (110, -141) ) C, 
a  /  $ 

+  12e[64h, -3 (38,-49)  IC5  +  3[224h, -8  (94,-123), 

-(338,-459)]C, 

-  88[4(14,-23),5(78,-113))C3  -  12  [80h, -2  (26, 9) , 

3(10,-23)]C2 

-  728(16(3, -4), -(98. -103))Cj  + 
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^(di/dv)  -  -  { 7eq'^(p/r)  +  3  cos  , 


(231) 


^(di/dv)  -  K^{36Q^  sin  v  -  2(p/r)  scC2}vj  , 


(232) 


and 


D  (di/dv) 

(D 


Djj  (di/dv) 


2K  (p/r)  scC^fDj 


scg  (eSj  +  3Sj  +  3eSj) 


(233) 


(234) 


There  are  four  terms  with  Qi  as  a  common  factor,  and  four  other  terms. 
As  in  section  5.3,  the  q"^  factor  in  the  first  ej  term  disappears  when  it  is 
combined  with  the  aj  term.  The  overall  combination  is  given  by 


D,(di/dv)  -  — K^sc{e^(2,5)i;,  +  lOe  (1,2)Z,  +  2  [2  (3,  5) ,  (4,  3)  ]  +  6e(3,2)Z^ 

2  48  6  5  4  3 


+  e'(6,l)i:. 


20“^  (26,-33)8^  -  20e(5,-7)"^ 


+  4(2(5, -3),  -  3(2,-3))S2  +  4e(31,-33)Sj  +  2eM14,-15)S^ 
-  4e^(3,-7)  sin  2v  -  8e(3,-7)  sin  v}. 


(235) 


The  term  in  So  is,  as  in  section  5.3,  a  source  of  long-period  variation, 
given  by 

“  7T  -  15£)S„  ,  ,236) 

with  the  inevitable  short-period  carry-over  term.  The  remaining  terms  of  (235) 
lead  to  the  following  formula,  in  which  the  term  in  Co  and  the  'pure  constant' 
are  chosen  to  suit  the  derivation  of  b2  in  section  6.4: 


- s  c  X 

288 

X  {e^2,5)rg  +  12e(l,2)r^  +  3(2  (3, 5) ,  (4, 3)  )r^  +  12e(3,2)r^  +  3e%b,l)r2 
-  3e^(26,-33)C,  -  40e(5,-7)C3  ^  12  [2  (5, -3)  , -3  (2, -3) ) 


+  24e(31,-33)Cj  +  4e  (9,1>C^  -  12e  (3,-7)  cos  2v 

•  48e(3,-7)  cos  v  +  2  [3  (17, -25) ,- (90, -91)  ]  }  . 


(237) 
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5 . 5  Chack  thus  far  by  parturbationa  in  p  and  p  eoa^  1 

The  absolute  constancy  of  pc^  ,  for  a  disturbing  function  associated  with 
only  the  zonal  harmonics,  has  been  referred  to  in  sections  2.2  and  4.2.  Now  that 
we  have  second-order  solutions  for  a,  e  and  i  .  we  can  derive  the 
corresponding  solutions  for  p  and  c^  ,  and  use  their  combination  to  check  that 
pc^  is  indeed  constant  to  second  order. 

Since  there  is  no  long-period  variation  in  a  ,  whilst  the  long-period 
variation  in  e  is  given  by  (225),  it  follows  at  once  that  the  long-period 
variation  in  p  is  given  by 

P2.I,f  "  -12  .  (223) 

We  obtain  pj  from  the  second-order  identity 

p  +  Kp^  +  K^pj  -  (a  +  ica^  +  K^a^)  {  (1  -  e^)  -  2Kee^  -  K  (e^  +  2ee^)  }  ,  (239) 

which  gives  (on  dropping  bars  as  usual) 

Fj  -  •  (240) 

From  (147)  and  (148)  we  get 

aej  +  2eaj  -  •—  aq"^  [f  {3e(l,3]  (eC^  +  6C^)  +  (28,133,  19)03  +  12e[13,7]C3 

+  3(4,55,l)Cj  +  3e[l,3)  (60^  +  eC_j)  } 

+  4h{  (1,3)  (e^  cos  3v  +  6e  cos  2v  +  3(4,1]  cos  v) 

+  2e(13,71  }]  .  (241) 

We  use  (148)  again,  to  obtain  the  product  of  (241)  with  ej  ,  then  subtract 
the  result  from  the  combination  of  82  and  62  specified  by  (240);  the  required 
form  for  32  was  suppressed  in  section  5.2,  it  will  be  recalled,  whilst  62  is 
given  by  (226) .  We  find,  after  the  subtraction,  that  most  of  the  terms  have 
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cancelled  out,  and  a  factor  ,  needed  to  convert  aq"^  in  (241)  into  p  ,  can 
be  extracted  from  all  that  remain.  The  resulting  expression  for  P2  is; 

Pj  -  -  p  [f^{e^rg  -  3[2,5ir^  -  48er3  -  SSeV^  } 

-  f  {3e^(26,-33)C^  +  8e<32, -41)03  24h[2,3)C3  -  72e(8,-9)Cj 

-  36e^{8,-9)Cg} 

-  2{6ef  (4,-5)  (e  cos  2v  +  4  cos  V)  +  [(8,-2‘j,l),3(8,0,-15))  }  ]  . 

. (242) 

From  (236)  and  (238)  it  is  immediately  confirmed  that  there  is  no  long- 

>  —2 

period  variation  in  pc^  .  Also,  expansion  of  cos^(i  +  hi]  +  K  ij)  leads  to 

(0^)3  =  (2f  -  Dij^-  2sci3  ,  (243) 

from  which  it  follows  that 

(pc^)^  =  c^p^  -  2psci3  +  ij{p(2f  -  Dij  -  2scp,  }  .  (244) 

From  (149)  and  (172)  we  get 

p(2f  -  l)ij  -  2scpj  -  -  ■i-psc{(X,2)  (eCj  +  30^  +  3eCj)  +  (11,-18)}  .  (245) 

Using  (149)  again,  for  the  product  of  (245)  with  i]  ,  and  adding  the  result  to 
the  required  combination  of  p?  and  i2  ,  given  by  (242)  and  (237)  respectively, 
we  eventually  obtain  the  final  result  we  see)t,  viz 

(PC^),  ”  -  -3^^-  pc2{e2f(63,-82)Cj,  -  [4(1, -6, 2),  (12,-50,13))  }.  (246) 
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The  absence  of  v-dependent  terms  in  (246)  coiqpletes  our  chec)c  of  the 
constancy  of  pc^  ;  on  inverting  our  stand^int,  we  have  a  powerful  check  on  the 
expressions  for  82  r  e2  and  i2  •  (In  practice,  a  number  of  errors  were  made  in 
the  e2  analysis  and  were  located  via  this  check.)  Expressions  (175)  and  (246) 
provide  the  second-order  relation  between  the  mean  and  osculating  values  of 
pc^  ,  connected  (to  second  order)  by 

pc^  -  pc^  +  X(pc^)j  +  K^(pc^)^  .  (247) 

In  the  light  of  (247)  and  the  absolute  constancy  of  pc^  ,  (246)  may  appear 

as  a  paradoxical  result,  since  the  first-order  variation  cf  Co  (interpreted  as 

~  -2  , 

COS  2(0  ,  with  (0  varying  secularly)  induces  third-order  variation  in  K  (pc')j. 

But  the  second-order  variation  of  p  and  f  in  (pc^)i  ,  given  by  (175)  with 

bars  added,  induces  third-order  variation  in  K(pc^)l  ,  and  this  may  be  regarded 

as  cancelling  the  variation  in  K  {pc^)2  -  However,  the  real  resolution  of  the 

paradox  stems  from  the  observation  that,  in  the  present  analysis,  it  is 

meaningless  anyway  to  regard  pc  ,  in  (247) ,  as  contact  to  third  order,  since  we 

have  not  carried  the  definitions  of  (semi-) mean  elements  beyond  the  second  order. 

5 . 6  Perturbation  in  Q 

The  starting  point,  from  the  exact  equation  (144),  is 

4“  -  -  K(n/n)  (qVq^)  {  (p/r)  V(p/r)2}Qo  ,  (248) 

dv  “ 

where  *  c (1  -  C^)  .  (249) 

The  six  contributions  to  the  second-order  component  of  dil/dv  may  be 
listed  as  usual,  viz 


D  (dll/dv) 

d 

-  ~  K^a'^p/r)Q^aj  , 

(250) 

(dXl/dv) 

”  K‘(p/r)s{l  -  Cj)ij  , 

(251) 

D  (dO/dv) 
e 

“  -  K^{7eq'’^(p/r)  +  3  cos  v}Q^^e^  , 

(252) 
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D^(dn/dv)  -  K^{3eC^  ain  v  -  2(p/r)cS2}Vj 


(253) 


D^(dfl/dv)  -  -  2K^(p/r)cS,Ci),  (254) 

CP  £  I 

and 

2 

D^(dfll/dv)  -  -  cg(eC3  +  +  3eCj)  .  (255) 

There  are,  as  in  section  5.4,  four  terms  with  as  a  common  factor,  and 

four  other  terms,  the  combination  (free  of  the  necessity  for  a  q'^  factor,  as 
usual)  being  expressible  by: 

D3(dn/dv)  »  -  k'^c  X 

X  {e^(4,3)r^  +  10e(2,l)rj,  +  214(3,1),  (8,-l))r^  +  6e(6,-l)rj 
+  e^(12,-5)r3  -  52e^(l,-l)C^  -  8c(ll,-8)C^ 

+  8((8,-15),  -(3,1)]C3  ^  Se(17,-30)C, 

+  4e^(7,-15)C^j  +  2e^(4,5)  cos  2v 

+  '2e(2,l)  cos  V  +  2(4  (1, -1) ,  (4,  b)  ]  }  , 

.  (256) 

The  final  term  in  (256)  is  responsible  for  the  first  appearance  of  second- 
order  secular  variation;  thus,  it  is  given  oy 

ftj  •  -  c(4(l  -  f)  +  e^(4  +  5f))  .  (?57) 

The  term  in  Cq  in  (256)  leads,  as  in  earlier  sections,  to  the  long-period 
variation  given  by 

^2  «p  “  -  —  c{l  -  \bt)C^  .  (258) 
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The  remaining  terms  of  (256)  lead  to  the  following  formula  for  Q2  ,  the 
•constant’  term  in  So  being  chosen  to  suit  the  derivation  of  b2  in 
section  6.4: 


X  {e^(4,3)i:g  +  12e(2,l)Z5  +  3f4  (3, 1) ,  (8,-1)  +  120(6,-1)1:3 

+  3e^(12,-5)r2  -  78e^(l,-l)S^  -  166(11,-8)83 

+  24((8,-15),-(3,l))S2  +  48e(17,-30)Sj  +  2e^  (18, -19) 

+  6e^(4,5/  sin  2v  +  72e(2,l)  sin  v} 

. (259) 


5 . 7  Perturbation  in  O) 

The  starting  point,  from  the  exact  (145),  is 

~  -  •^Ke'Nn/n)  (qVq^)  {  (p/r)^/ (p/ f)^  }  Q„  ,  (260) 

dv  >’ 


where  0.^  -  f(5eC^  +  I4C3  +  2eC,  -  2Cj  +  eC^  +  4e) 

-  4h(2eC2  -  e  cos  2v  -  2  cos  v  -  3e)  .  (261) 


Equat.ion  (261)  has  a  certain  inconformity,  in  comparison  with  (215)  ,  similar  to 
that  remarited  upon  in  the  first-order  analysis  in  relation  to  (153)  (propagated 
from  (145))  since  we  do  not  have  the  neat  division  in  which  the  f-terms  involve 
precisely  the  Cj  and  the  h-terms  involve  precisely  the  cos  jv  . 

The  six  contributions  to  the  second-order  component  of  dCi)/dv  may  be 
listed  as  usual;  thus 
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D  <d<o/dv) 


(262) 


-  -  K^a'^e'^  (p/r)Q  a  , 

16  CD  1 

Dj,{d(0/<lv)  =  •l-KV^(p/r)  sc  (5eC^  +  14C3  +  1460^  -  2C^  + 

-  2  (3e  cos  2v  +  6  cos  v  +  7e)j  i.  , 

(263) 

D  (du/dv)  “  [7eq"^(p/r)  +  3  cos  v]  -  ie'^(p/r)Q  }e  ,  (264) 

o  °  0)  (i>e  ^ 

where  Q  “  f (7C  -  C  )  +  4h  cos  v  ,  (265) 

e)e  3  1 

D  (du/dv)  -  -  — K^e“^{3eQ  sin  v  +  2(p/r)Q  }v  ,  (266) 

V  8  to  tov  1 

where  Q  “  fdOeS,  +  21S,  +  2eS,  -  S.)  +  4h(e  sin  2v  +  sin  v  -  2eS,)  ,  (267) 

(Ov  4  3  2  1  2 

D  (dco/dv)  -  -i-K^e'^p/r)  {f (5eS.  +  14S,  +  2eS,  -  2S,  +  eSJ  -  8ehS,}(0,  (268) 

4  43^iJ  /  i 

and 

D.  (dco/dv)  -  -  {  3e^fC,  +  18efC,  +  [28f ,  -  (8,  - 19)  )  C, 

<0  12  5  4  J 

-  12e(2,-5)C2  -  3[4f,  (8,-15)  ICj  -  18efCQ  -  Se^fC,^  }  (269) 

There  are  four  terms  with  Q<o  as  a  common  factor,  ard  five  other  terms, 
the  combination  (free  of  as  usual)  being  expressible  by: 

D2(dt0/dv)  -  -  ^  K^e"'  [ef‘{l5eV,g  +  162eV^  +  8e[82,15)rg  +  14[84,67)r,  } 

+  (784£^,2<4$f’,  -  (64,16,-277)  ]  +  20e  [10  5  f  -  (16,  0, -69)  ) 

-  8e^(2(24, -4, -1C3),  (32,-20, -23>ir^  -  12e[18f\  (48, -32, -35) 

-  [48f^208f^  (192,-176.-9)  3r^  -  ef*^  {  2  [3b,  19]  -  ISe^P  ^  -  3eY_,  }  + 
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+  8  {  4e^hfCg  +  42e^hfC^  +  e^f  [172h, -3  (14, -17)  ]Cg 

+  5ef[64h,-(34,-39))C5  +  [224hf, -12f (2, 3), 

(104,-242,159)]C^ 

+  8et2f  (21,-31),  (22,-17, -5)]Cj  -  4e^  [2  (16, -138, 141) , 

-  (12,  22,  -33)1C2 

+  8e(f  (4,-3),-(34,-173,150)]C^  +  2e^  [2f  (14, -15) , 

-  (28,  -158,  135)]  Cg 

+  ef[48h, (110,-129) ]C_j  +  e^£ [52h,  (34, -39) )C  ^ 

+  18e^hfC  ,  +  2e'’hfC  ,  } 

“3  “4 

+  2  { 6^(8)/  +  9f^)  (e  cos  6v  +  10  cos  5v) 

+  4e^[2(40h^  +  37f^), -(8,  8,-37) )  cos  4v 
+  6e(12(8h^  +  5f^),-{24,24,~137)  1  cos  3v 
+  tl6(24h^  +  7f^, -32(4, 8, -41), -(264, -136,  -301))  cos  2v 
+  4e(2(8h^  +  43f^), -(256, -400, 9)  1  cos  v 

-  4e^[2(64,-180,95),  (56,  -36,-45))  }]  . 

. (270) 

The  final  term  in  (270)  represents,  as  in  section  5.6,  the  secular 
variation,  which  is  thus  given  by 

[2(64  -  180f  +  95f^)  +  e^(56  -  36f  -  45f^  ]  .  (271) 

96 

The  term  in  Cq  leads,  as  usual,  to  the  long-period  variation  given  by 

“7  p  “  "TT  t2f(14  -  15f)  -  e^(28  -  158f  +  135f^)]c.  .  (272) 

w  f  C  p  4  8  V 
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The  retnalning  terms  of  (270)  lead  to  the  following  formula  for  (O2  >  the 
•constant*  terms  in  Xq  So  being  chosen  to  suit  the  derivation  of  W2  in 

section  6.4: 

(0  - - L-e‘2  X 

2  4608 

X  [SefMse^Zjg  +  36e^E^  +  2et82,151Eg  +  4[84,67)£^} 

+  [784f^,2448f^,-(64,16,-277)  )E,  +  24e(106f^  -  (16, 0, -69) ) E, 

o  b 

-  12e^[2(24, -4, -103),  (32,-20, -23))E^  -  24e  [18f\ (48, -32, -35) )  E^ 

-  3(48f^  208f^  (192,-176,-9)1  E^ 

-  3ef^{  4  [36,191Ej  +  2e[66,-l)EQ  +  36e^E_j  +  3e^E_^ } 

+  8ef  {3e^hS„  +  36e^hS,  +  e{172h, -3 (14, -17))  S,  +  6 [64h,  -  (34, -39)  1  S,  } 

01  o  5 

+  4{3(224hf,-12f  (2,3),  (104, -242, 159)  ) 

+  32e[2f (21,-31), (22, -17,-5) ISj  -  24e^ (2 (16, -138, 141) , 

-(12,22,-33)  152 

+  96e[f  (4,-3), -(34,-173,150))Sj  -  [288hf ,  16f  (116, -147) , 

(144,-418,311)  )Sp} 

-  24ef  {2[48h,  (110,-329)lS,j  +  e(52h,  (34, -39)  ) 

+  12e^hS_3  +  e\s_^} 

+  2{e^(8h  +  9f  )  (e  sin  6v  +  12  sin  5v) 

+  6e^(2(40h^  +  37f^) , -  (8,  8, -37) )  sin  4v 

+  12e(12(8h^  +  5f^, -(24, 24, -137)1  sin  3v 

+  3(16(24h^  +  7f^, -32(4, 8, -41), -(264, -136, -301)  j  sin  2v 

+  24e(2(8h’  +  43f^  -  (256, -400,  9)  1  sin  v}]  . 

. (273) 
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5 .  t  V*rtuxb«tioa  la  a 

The  starting  point,  froa  the  exact  equation  (146),  is 

~  «  --^Ke'^Cn/n)  (qVq^)  {  (p/r)^/ (p/r)^}Q  ,  (274) 

dv  -  o 

where  »  £(5eC4  +  14C3  "  18eC2  -  2Ci  +  eCo)  +  4h(e  cos  2v  +  2  cos  v  -  3e)  . 

. (275) 

The  six  contributions  to  the  second-order  component  of  do/dv  may  be 
listed  as  usual;  thus 


C^(dff/dv)  =  ^K^a'^e"'q(p/r)Q^a^  ,  (276) 

(do/dv)  =«  -■iK^e"^q(p/r)  3  c  (5eC^  +  140^  -  18eC^  -  2C^  +  eC^ 

-  6  (e  cos  2v  +  2  cos  v  -  3e))  ,  (277) 

b  (dO/dv)  »  -iK^e''q{Q  [6eq'^(p/r)  )  3  cos  v'  -  2e'^{p/r)Q  }e  ,  (278' 

e  8  <j  •  (Je  ^ 

where  -  f(7C3  -  Ci)  +  4h  cos  v  (  ~  Q^^)  ,  (279) 

D^, do/dv)  -  •i'K^e‘^q{3eQg  sin  V  +  2(p/r)Q^^}vj  ,  (280) 

where  fdOeS^  +  2IS3  -  18eS2  -  Si)  +  4h(e  sin  2v  +  sin  v)  ,  (281) 

n^(dO/dv)  -  ■i.KV^q(p/r)  f{5eS^  +  I4S3  '  '  ^®)  ^  ^  “1  <282) 

and 

L^(dO/dy)  -  K^e'’q£g  {  3e^C3  +  18eC^  +  (28,-13103  -  36eC3 

-  3[4,17)Cj  -  18eCg  -  3e^C_^  }  .  (283) 


There  are  four  tenas  with  Q<j  as  a  common  factor,  and  ive  other  terms, 
but  (because  the  first  0^  term  has  coefficient  6  instead  of  the  usual  7)  it  is 
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not  possible  this  time  to  extract  the  factor  when  this  term  is  combined  with 
the  Da  term.  The  overall  combination  is  expressible  by: 

D  <d0/dv)  -  X 

^  1536 

X  [f^{l5e^(2,-ljr^Q  +  162e^t2,-l]rg  +  2e^[65S,-208,-97]rg 

+  2etll76,350,-737ir  +  [1568, 4112, -3814, -21) r 
»  6 

+  4e[1060,-704,-203)r  +  60e^  [16, -32, -5) T, 

5  H 

-  36e[12,8.35)r.  -  [96, 368,  350, 473) P,  -  2e [72, 2, 121)  T, 

^  i  L 

-  18e®rg  +  18e^[2,-l)r_^  +  {2,-l]r_^} 

+  8f  {4e^h[2,-l)Cg  +  42e\[2,-l)C, 

+  e^[344h.-8(32,-45),  (86,-105)  )C^ 

6 

+  e[640h,-30(22,-29),  (382, -4 53)) 

+  (448h, -8 (34, -33), 2(298, -327), -(418, -531) JC, 

+  2e [16 (21, -31), 4 (60, -61), -(532, -749) 

+  2e^  [16  (17, -21), -4  (118, -147), -(10, -63) 

+  2e(8  (4,-3), -8(95,-111),  (506, -531)  )Cj 

+  4e^[2(14,-15),-4(44,-51),  (106, -111) ) 

+  e[96h, 2(86, -93), -9(26, -31))C_j 
+  e^[104h,lS,-3(22,-25)  ]C_2 

+  18e^h[2,-l)C  ,+  2e'’h[2,-l)C  ,}  + 

-J  -  4  ' 
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+  2{e^(8h^  +  9f^)[2,-l)(e  cos  6v  +  10  cos  5v) 

+  2e*[8(40h^  +  37f^)  ,-8  (24, -56, 45) ,  (24, 56, -143))  cos  4v 

+  6e(24(8h^+  5f^), -2(72, -120,1),  (40,72,-253))  cos  3v 

+  [32  (24h^  +  7f^), -16  (40, -40, -103), 

2(312, 88, -1331), -(504, -744, 101)}  cos  2v 

+  4et4(8h^  +  43f^),4(76,-156,-17), 

-(488,-1176,689))  cos  v 

+  2e^(16a6, -56,31), -16(16, -88, 95), -(280, -328, -79)  )}  ]  . 

.  (284) 

The  final  term  at  once  yields  (with  the  usual  notation)  CT2  >  and  the  long- 
period  variation  is  given  by  the  terms  in  Pq  and  Co  .  .  However,  there  is 
little  point  in  giving  the  explicit  expressions,  or  in  giving  the  formula  for  the 
short-pericd  Oj  /  since  o  (from  its  definition  in  section  2.1)  is  merely  an 
intermediate  element  that  has  to  be  analysed  as  part  of  the  analysis  for  M  .  The 
existence  of  a  non-zero  long-period  term  in  Po  ,  in  addition  to  the  usual  term 
in  Co  ,  is  another  aspect  of  the  ephemerality  of  O  ,  and  it  will  be  found  in 
sections  5.9  and  5.10  that  in  the  M-analysis  the  Pq  term  is  cancelled  by  an 
equal  and  opposite  term  from  the  quantity  1  defined  by  (3) . 

5 . 9  Perturbations  In  n  and  J 

The  second-order  perturbation  in  n  ,  lilce  the  first-order  perturbation,  is 
purely  short -periodic  and  may  be  obtained  from  the  perturbation  in  a  by  use  of 
the  appropriate  version  of  Kepler's  third  law  for  mean  elements.  The  generalized 
second-order  form  of  this  law,  with  (non-dimensional)  off-sets  (fj  and  42  to 
suit  the  definitions  of  n  and  a  ,  may  be  written 


-2_3 
n  a 


(285) 


but  we  have  seen  in  section  4.1  that  the  natural  first-order  definition  of  a 
and  fi  leads  to  the  tajeing  of  to  be  zero.  There  is  no  reason  to  expect  42 

to  be  zero  as  well,  however,  and  it  will  appear  in  section  5.10  that,  with  a 
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defined  to  be  Identical  with  the  exactly  conatant  a*  ,  and  n  chosen  to  second 
order  such  that  the  secular  variation  in  M  vanishes, 

-  8f  -  5f^)  .  (286) 

<  24 

If  we  write,  for  convenience,  aj  and  nj  for  aj/i  and  n^/n  (j  “  1 
and  2) ,  it  follows  that 

n^a^  -  n^a’’  (l  +  Ka^  +  ic^a^)  (l  +  Kn^  +  K^n^)  .  (287) 

From  (1),  (285)  and  (287),  we  have  at  once  that 

-  -  (3aj  +  2nj)  (288) 

and 

“  SUj  +  ea^n^  +  Sn^  -  Sa^  -  2n^  .  (289) 


Equation  (288)  just  leads  to  (159),  on  the  basis  (from  section  4.1)  that 
)il  ”  0  ;  the  new  material  is  in  (289),  which  with  ni  set  to  "7  3]  yields 


A 

n 


2 


2  ^2 


A  2 


(290) 


The  appropriate  expression  for  aj  is  given  by  (158)  .  when  this  is 
squared  and  combined  with  the  expression  for  $2  given  by  (213)  ,  then  foe  )l2 
given  by  (286)  it  is  eventually  found  that 

"2  ”  v)^  X 

X  [3f^{3e^rg  +  24e¥.,  +  12e^[8,-l)rg  +  8e[22,-l]rj  +  2  [56,  32,  17) 

+  24e(2,5]r3  *  4e^(8,13ir3  +  240^3  +  3e''rp }  + 
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+  24f{eSCg  +  Be^hCj  -  2«^t8(l,-l)  ,-(22,-29)  ]C^ 

-  4e[4(4,-5),-(30,-41))C3  -  [16h, -32  (3, -4) ,  (IC .  1)  1C^ 

+  4e(4(8,-9),-3<6,-5))Cj  +  2e^  [8  (5, -6) ,  -  (26, -27) ) 

+  8e^hC_j  +  e^hC_2} 


+  2{e^(8h^  +  9f^)  (e  cos  4v  +  8  cos  3v) 

-  4e^  (96f  (1,-2), -(56, -72,-3)  ]  cos  2v 

-  8e(2(32,-48,-27),-3(40,-88,45))  cos  v 

-  (8(40h^  -  33f^), -96(4, -16, 15), -3(72, -88, -13)!  }  ] 

+  7rnq^8  -  8f  -  Sf"’)  . 

'IS 

. (291) 

The  second-order  perturbation  in  M  ,  to  be  derived  in  section  5.10,  is 
given  by  the  combination  of  the  perturbation  in  a  with  the  perturbation  in  /  , 
the  integral  of  n  ;  since  (i.i  section  5.8)  we  only  took  the  analysis  as  far  as 
D2(da./dv),  however,  what  we  need  here  is  D2(dJ/dv),  to  combine  with  D2(dCT/dv)  to 
give  ^2  (d'rt/dv)  . 

Now  (291)  displays  n2  as  a  sum  of  many  terms,  all  with  (1  +  e  cos  v)^ 
as  a  factor,  together  with  a  final  term  that  emanates  from  )l2  .  Apart  from  a 
complication  to  be  dealt  with  in  the  next  paragraph,  the  'sum  of  many  terms' 
gives  an  immediate  contribution  to  D2  (dj/  Iv) ,  since  (with  bars  omitted  as  usual) 


n(l  +  e  cos  v) 

by  (39)  and  (17).  Thus  the  factor  (1  +  e  cos  v)^  in  (291)  is  conveniently 
canoelled  and  the  required  contribution  from  the  sum  of  many  terms  is  given  on 
replacing  nq"^  (1  +  e  cos  v)^  by  K^q"^  .  This  just  leaves  the  'final  term', 
which  contributes  -i- K^q^  W  '(8  -  8f  -  5f^)  ,  since  dv/dM  =  W  . 


dv  nW 
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The  con^lication  in  the  ’sum  of  liiny  terms',  referred  to  in  the  last 
paragraph,  arises  because  to  cover  some  carry-over  terms  from  ni 

as  well  as  the  direct  terms  associated  with  n?  .  The  carry-over  terms  arise 
from  the  first-order  secular  variation  in  (S  ,  and  constitute  a  final 
manifestations  of  the  ’Do)  source*  of  terms  introduced  in  section  5.2  and 
encountered  in  all  the  subsequent  sections.  Here,  the  origin  of  the  terms  is  the 
first-order  evaluation  of  by  v-integration  of  <160).  As  in  section  5.2,  to 
con^nsate  for  the  neglect  of  the  term  (-2K5i/j2)Sj  ,  when  replacing  each 
integral  of  Cj  by  Sj/j  ,  we  need  a  contribution  to  D^idJ/dv)  given  by 
^qf(4  -  5f)  (eC  +  3C  +  3eC, )  ,  and  the  effect  of  this  on  the  conversion  of  (291) 
into  the  expression  required  for  D2  (dj/dv)  is  to  replace  the  terms  in  Cj,  C; 
and  Cl  (within  the  '24f  brac)cet')  by 

-4e(2(4,-5),-(22,-31)  ICi  +  [16(5,-6)  ,-8*^, -(10, 1)  )C2 
+  4e[2(28, -33)  ,-3(14,-15)  )Ci  . 

5.10  Perturbation  in  K 

The  expression  for  D2(dM/dv)  is  given  at  once  by  the  combination  of  (284) 
with  the  expression  for  D2(dl/dv)  obtained  from  (251)  as  indicated  in  the  last 
section.  The  result  is: 

D,(clM/dv)  -  X 

2  1536  ^ 

X  [f^{l5e‘'(2,-lirjp  +  162e^[2,-l)rg  +  16e"  [82, -26, -11) 

■*•  14e(168,50,-95]r^  [1563,  4112, -3238, -93) )'g 

+  20et212,-8B,-43]rj  +  96e^  [17, -16, -1)  r  ^ 

-  lO80[4,O,5]r^  -  [96,368, 158, leijT^  -  2e (72, 2,  49) 

+  18e"[2,-l)r_j  f  Se"  (2.-l)r_2} 

+  8f  {4e^h[2,-i:Cg  V  42e^h[2, -1)C^ 

+  2e^[172h,-4  ;32,  -45),  (66,-57))C, 

+  I448h, -8(3-1, -33), 2  (250, -279), -{154, -183)]C^  + 
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+  2e[16(21,-31),4(36,-31),-(318,-377)lC3 
+  8e^[4  (32,-39), -(118.  -141), -5 
+  2e(8(4,-3),-8(ll,-12),  (2,9)]C^ 

+  8e^q^n4,-lS)CQ  +  e[96h,2  (86, -93' , -3  (62, -69)) 

+  2e^[52h,8,-3(10,-ll)  JC_2  +  18e^h  (2, -1) 

+  2e'*h[2,-l]C  ^ 

+  2{e^(8h^  +  9f^)[2,-l)(e  cos  6v  +  10  cos  5v) 

+  8e^(2(40h^  +  37f^) , -2 (24, -56, 45) , (8, 8, -29) )  cos  4v 

+  2e(72(8h^  +  5f^), -6(72, -120,  1),  (184,24,-543))  cos  3v 

+  (32(24h^  +  7f^), -16(40, -40, -103), 2(312, -296, -563), 

-(56,-168,125))  cos  2v 

+  4e(4(8h^  +  43f^),4(12,-60,37),-  (8,-120,149))  cos  v 

-  16eV(8,-8,-5)  }] 
+  -i-K^qV'  (8  -  8f  -  5f^)  . 

40 

. (293) 

The  basis  for  the  choice  of  )l2  defined  by  (286),  which  detern.ined  the 
final  terms  of  both  (291)  and  (293),  should  now  be  apparent,  since  the 
combination  of  the  final  term  in  (293)  with  the  last  of  the  preceding  collection 

of  terms  is  ■^K^q^(8  -  8f  -  5f^)  (w’^  -  1)  ,  the  mean  value  of  which  (with 

respect  to  v  )  is  zero  oy  (39) .  The  penultimate  term  of  (293)  in  fact  defines 
M2  ,  viz 

-  8f  -  5f^)  ,  (?94) 

*  48 
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^  A 

and  ]l2  was  set  to  2M2  to  allow  for  the  resulting  secular  effect,  as 
anticipated  in  the  last  section  and  (before  that)  in  section  3.1.  Though  our 
choice  of  H2  i''  principle  allows  for  M2  ,  making  n  and  Mjec  identical  as 
required,  we  must  not  overlook  the  short-period  carry-over  effect  of  K2  .  This 
can  be  allowed  for  by  applying  (64)  with  ^  =  M  . 

As  forecast  in  section  5.8,  (293)  contains  no  term  in  Pq  .  The  term  in 
Co  leads,  as  usual,  to  the  long-period  variation  given  by 

.  (295) 

The  remaining  terms  of  (293)  lead  to  the  following  formula  for  K2  ,  the 
'constant'  terms  in  £0  and  So  being  chosen  to  suit  the  derivation  of  in 

section  6.2: 


”2  ”  "5^  ‘  [f^{9e''(2,-H£jg  +  108e^t2,-l]£g  +  12e^[82,-26,-U)£g 

+  12e[168,50,-95)£,  +  (1568,  4112, -3238, -93) £g 

+  24e(212, -88, -43)15  +  144e^ [17, -16, -1) -  216e [4, 0, 5) 

-  3(96, .368, 158, ISDEj  ’  12e  [72, 2,  49)  £j  -  12e^  [66, -34 , 13)  £, 

-  108e^[2,-l)£_,  -  96"  [2, -!)£_,} 

+  4f  {6e^h[2,-l)S„  +  72e^h[2,-l)S., 

+  4e^il72h,-4  (32,-45) ,  (46,-57)  )  S, 

+  12e[128h,-6(22,-29),  (86,-105)  jS^ 

+  3  [448h, -8  (34, -33),  2  (250, -279), -(154, -183)  )S^ 

+  8e [16 (21, -31), 4 (36, -31), -(318, -377) ) S3 
+  48e^[4(32,-39),-(118,-141),-5(2,-3))S2 
+  24e[8(4,-3),-8(ll,-12),  (2,9)  )Sj  - 
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-  (576h, 16 <214,-267) , -2(802, -963) , - (830, -861) J 

-  12e[96h,2(86,-93),-3(62,-69)lS_j  -  12e^[52h,8,-3(10,-ll) 

-  72e^h[2,-l)S_3  -  6e‘*hl2,-l]S_^} 

+  2{e^(8h^  +  9f^)[2,-l](e  sin  6v  +  12  sin  5v) 

+  12e^[2(40h^  +  37f^) , -2  (24, -56,  45) ,  (8,  8, -29) )  sin  4v 

+  4e(72(8h^  +  5f^),  -6(72,-120,1),  (184,24,-543)  1  sin  3v 

+  3(32(24h^  +  7£^), -16(40, -40, -103), 2(312, -296, -563), 

-(56,-168,125)  1  sin  2v 

+  24e[4(8h^  +  43£^), 4(12, -60, 37)  ,-(8,-120, 149)  1  sin  v}] 

. (296) 


6  ADOPTED  SOLUTION  FOR  PERTURBATIONS 

6.1  Secular  and  long-period  terms,  plus  s)tort-period  carry-over 

Solution  of  the  Jj  problem  to  second  order  amounts  to  the  extension  of  the 
first-order  solution  of  section  4  by  additional  terms.  The  terms  expressing  the 
long-term  variation  in  the  mean  elements  ate  compact  and  small  in  number  and, 
since  they  are  not  easily  to  be  found  among  the  much  longer  expressions  of 
section  5,  it  is  convenient  to  repeat  them  here.  Though  the  Cj.fp  ate  actually 
all  non-singular,  propagation  of  the  •>  will  nevertheless  be  via  3"^ 

L2,fp  (rather  than  C52,4p  3"^^  **  explained  in  section  3.5  and  Part  3, 

so  we  inc.lude  here  also  the  expressions  for  these. 

Each  C2,io  induces  a  short-period  carry-over  (cf  (74)),  and  these  are 
dealt  with  by  combining  them  into  components  of  the  perturbations  in  the 
spherical  coordinates,  as  indicated  in  section  3.2;,  the  resulting  expressions 
(components  of  5r ,  5b  and  5w  )  are  included  in  the  present  section.  The 
pure  short-period  components,  representd  by  rp,  bp  and  wp  ,  are  derived  in 
sections  6.2  to  6.4. 

'ne  only  secular  terms  are  in  fl  and  TO  .  They  are  given  by  K^nftpt  and 
K  nTO2t  where,  to  repeat  (257)  and  (271)  with  bars  added. 
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FrcHB  (302) -(304)  we  get  the  expressions  to  replace  (303)  and  (304)  in  practice, 
viz 

^2  tp  ”  (1^  -  15f)  cos  2a  (30b) 

and 

L,  ,  -  --2-i^qf(14  -  15f)  cos  2b.  (306) 

2««p  48 

The  short-period  effects  induced  by  (300)-{304)  are  combined  into  spherical 
coordinate  perturbations  by  use  of  (95) -(97).  The  results  may  be  expressed,  in  a 
new  but  obvious  notation  (in  which  the  carry-over  r^^o  etc  will  be  added 
directly  to  the  rj  etc  to  be  developed  in  sections  6.2  to  6.4),  as 

'2co  “  -  15f)  ms,  ,  (307) 

b^co  -  -  e^sem  [ISfy,  -  (28  -  45f)c_,]  (308) 

and 

“200  "  ■  15f)m{eC2  +  4C,  -  4eCp)  .  (309) 

(In  equation  (43)  of  Ret  4,  the  factor  e^  s  c  is  missing  in  the  terms  of  b./f,o  .) 

6 . 2  Pure  short-period  perturbation  in  r 

Since  r  is  a  function  of  the  three  elements  a,  e  and  M  ,  with  first- 
order  partial  derivatives  given  by  (41),  evaluation  of  rj  in  principle  entails 
the  combining  of  nine  terms  from  a  Taylor  expansion,  three  of  them  involving 
e2  and  M2  through  first-order  derivatives,  whilst  the  other  six  involve  the 
second-degree  products  of  ai,  ej  and  Mj  through  second-order  derivatives. 

Thus  the  analysis  is  potentially  much  more  complicated  than  for  v?  and  b2  (to 
follow  in  sections  6.3  and  6.4),  each  of  which  involves  only  a  pair  of  orliAtal 
elements  (not  three)  and  hence  only  five  Taylor  terms  (two  plus  three) .  A 
somewhat  different  approach  was  followed,  therefore,  ’tailored'  rather  than 
'Taylored';  it  is  based  on  the  use  of  eccentric  anomaly,  r  being  t  hs  dmpxO 
function  of  a,  e  and  E  given  by  equation  (78) . 

We  need  expressions  for  Ej  and  £2  ,  defined  to  represent  pure  short- 
period  perturbations  ir  the  usual  way,  and  these  can  be  obtained  from  (75)  .  We 
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r«plac*  *  by  (a  -f  Rai  -f  K  a2)  in  this  equation,  and  likewise  for  M  and  E  , 

_? 

then  expand  sin  E  (and  thence  e  sin  E)  to  0(K  )  ,  and  can  thus  identify 
expressions  equivalent  to  Mx  and  Mj  .  Rearrangement  of  these  expressions, 
with  algebraic  sisplification  through  the  introduction  of  v  where  possible, 
gives  us 


(310) 


where  e,  q  and  v  are  of  course  to  be  read  as  barred,  and 

Ej  -  (a/r)M2  +  q’^e^  sin  v  +  )4(a/r)Ej  [e^  (cos  v  +  cos  E)  -  q*^  eM^  sin  v]  . 

..,....(311) 


Equations  (310)  and  (311)  can  now  t>e  used,  as  required,  in  the  expressions 
obtained  by  expansion  of  (78) .  At  the  first-order  level  this  expansion  gives 

Tj  “  (r/a)aj  -  ae,  cos  T  -  aeEj  sin  E  ,  (312) 

from  which  the  use  of  (310)  leads  at  once  to  a  formula  for  ri  that  is 
equivalent  to  (95) .  At  the  second-order  level  tlie  expansion  yields 

-  (r/a)a2  -  ae^  cos  E  +  aeE^  sin  E  -  aj(e^  cos  E  -  eE^  sin  E) 

+  )4aEj(2ej  sir  E  +  eE^  cos  E)  ,  (313) 

from  whi.;h  rhs  ess  ot  (310)  and  (311)  leads  to 
-  (r./a)3.  -  =.  (a,  ccs  v  -  ncr'^H,  sin  v) 

*  ^  i  _  -  * 

+  «6q*'' [Sj  Sin  V  T  q"' M,  (1  +  e  COS  v)  ]  X 

X  [e.  ;4  sin  v  +  e  sin  2v)  +  eq”“Mj  (e  +  2  cos  v  +  e  oos  2v)  ] 

■t  aeq  '  sin  v  -  ae^  cos  v  . 


(314) 
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The  five  terns,  of  varying  coii()lexity,  in  (314)  have  been  presented  in  an 
order  that  reflects  the  way  in  which  they  are  most  naturally  combined.  To  start 
with,  we  combine  the  first  two  terms.  Since  <r/a)  is  sin^Jly  q2(l  +  e  cos  v)-i, 
the  lirst  term  is  given  by  (213)  with  the  overall  coefficient  replaced  by 
(1  +  e  cos  v)  .  Also,  (148)  and  (164)  give 

288 

ej  CdS  V  -  eq'*Mj  sin  ^  (3e^C^  +  12eCj  +  2[5,41C2  +  12eCj  +3e^Cg) 

+  4h(e^  cos  2v  +  4e  cos  v  +  3)}, 
. (315) 


the  product  of  which  with  (158)  yields  an  expression  for  the  second  teim  of 
(31-»),  with  (1  +  e  cos  v)  as  a  factor.  To  avoid  another  long  expression  here 
(and  likewise  at  some  other  points  of  section  6)  we  give  the  result  in  a  skeleton 
form  in  which  only  the  'end  terms'  are  quoted;  it  is 


-•j^aq'^/l  +  cos  V)  [3f^{3e^rg  +  ...  +  30^^^}  +  8fh{3e''Cg  +  ...  +  3e‘'c_2  } 

r  2{e^(8h^  +  9f^)  cos  4v  +  ...  +  8h^[12,22,l)  +  3f'’ (20, 74, 1 1 )  }  ] 


The  first  two  terms  of  (314)  can  now  be  combined,  whereupon  a  considerable 
simplification  occurs  since  ten  pairs  of  terms  cancel  out  and  a  factor  4q^  can 
be  extracted  from  the  remaining  terms.  The  result  is  worth  expressing  in  full; 
thus 

~  ajfCj  cos  V  -  eq”*Mj  sin  v) 

a(l  +  e  cos  v)  r3f^{e^r  +  3er  +  2q^r  -  5er  -  } 

-  8f  {e^(7,-9)C,  +  e(16,-21)C,  +  [2  (2 . -3) ,  -  (4, -3)  ] C, 

n  j  C 

-  e(20,-21)Cj  -  e^U,-12)Cg} 

-  2{e^(3, 0,-15)  cos  2v  +  e(32,-48,-3)  cos  v 

+  [(16,-48,27),  (3,0, -15))}]  . 

. (315) 
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Turning  to  the  third  term  of  (314),  we  find  that 
Oj  sin  V  +  q"^  Mj  (1  +  e  cos  v) 

-  -  ~  [f  {ae^Sj  +  6e(3,-l)S^  +  [28,-13183  +  ISeq^S^  -  3[«,1]S, 

-  6e[3,-l]S-  -  3e^S  .  } 

O  “A 

+  4h{e^  sin  3v  +  2e(3,“ll  sin  2v  +  [12,-7]  sin  v}] 

. (317) 

and 

e,  (4  sin  V  +  6  3in  2v)  +  eq”^  M,  (e  +  2  cos  v  +  e  cos  2v) 

[£{3e^Sg  +  66(3, 1]  $5+  2(14,19)8^  +  26(25,11183 

-  48q^8j  -  26(41, -5]8j  -  6(2,918(j  -  6e(3,l]S.j  -  36^5.^} 

+  4h{e^  sin  4v  +  2e(3,ll  sin  3v  +  2(6,5)  sin  2v  +  2e[ll,l)  sin  v  }  ]  . 

. (318) 


Multiplying  the  product  of  these  two  expressions  by  a  further  factor  14aq'^  ,  we 
get  the  required  'third  term';  in  skeleton  form  the  resulting  expression  is 


ae'V^  +  ...  +  Se^r,}  +  Shf{3e’C.  4  ...  +  3e'c.J 

+  2  {  e^  (8h"  T  5f  )  cos  7v  +  ...  -  6ef2  il68h^  +  73f^),  -  (fi8h^  -  197f^ ) , 


-4  {6h^  t  25f  )  1  }j  . 


This  expre.ssion  combines  well  with  the  fourth  term  of  (314),  given  by  multiplying 
the  SKpressicn  tax  Ms  ,  viz  (296),  by  aeq’’  sin  v  ;  the  result  of  that 
multlpXlcaticn,  in  skeleton  form,  is 
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-  j~  ae'V^  [fM9e"l2,-l]r^,  +  ...  +  9e^(2,-l]r_3} 

+  4f  {6he*[2,-llCg  +  ...  +  6he^ t2,-l]C_5 } 

+  2{e^[2,-l]  <8h^  +  9f^)  cos  7v  +  ...  -  24e[4(8h^+  43f^) ,  4  (12, -60,  37), 

-  (8,  -  120,  149) )  }]  .  * 


The  virtue  of  combining  the  third  and  fourth  terns  of  (314)  is  that  a  common 
factor  emerges,  to  cancel  the  existing  q"^  in  both  terms;  the  resulting 

expression  (in  full)  is 


-  ae'^  [f^{9eVj,  +  108e¥jg  +  3e‘^[164,29]r^  +  24e(42,31)r3 

+  [784, 2204, -267)r^  +  36o(64, -SSJPg [112, -1944, -325)r3 
-  16e(93,25ir^  -  [528, -328, -347) -  12e [96, -55) T, 

+  9[16,-164,-i9irj  +  24e[18,-4l)r3  +  3e^ [132, -67] r_j 

+  108eV_3  +  9e'’i\J 

1-  4£{6e'’hCj  +  72e\Cg  +  e^[344h,  -  (214, -273)  jC.^ 

f  24e[32h,-3(i3,-23)C,  +  [67  2h, - 4  (478, -609) ,  (558, -54  5) C,. 

+  4e(8(10,-13),  (798,-925)]C^  -  [288h, -16  (435, -519) , 

-(2,-147)  ]C, 

+  16er2(38,-4l),-(133,-139)  JCj  -  (672h,  8  (914, -1047) , 

(1195,-1413)  :C, 

-  36e[8(10, -13), (50,-59) JO^  +  [288h, 4 (398, -381) , 

V  « 

5 (iC2,  -105) ! C_j 

+  24e[24h,  (82,-37)c  ,  f  5e^ (1045.,  (11 0,  - 1’ 7)  ] C  , 

f  72e^hC_,  S=''hC.5}  +  I 

i 

-  ! 

?  I 

-  } 

J  ( 

I  ! 
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+  2{e^I8h^  +  9£^I  cos  7v  +  12e^[8h^  +  9f^]  cos  6v 
+  e^tl2(40h^  +  37f^), -(136, -24,-213)]  cos  5v 
+  8e[18<8h^  +  5f^>, -(136, -120,-123)1  cos  4v 
+  C48(24h^  +  7f^), -12 (248, -424, 5),  (232,-312,-81)1  cos  3v 

-  4e[48(16h^  -  9f^) ,  -  (232, -696, 351)  1  cos  2v 

-  (16(72)1^  +  37£^), -24(104. -152, -61),  (104,-312,-97)  1  cos  v 

+  3efl0(24h^  -  37fS  ,  (8, 264, -275)  1  }]  . 

The  terms  in  Fi,  F-i,  Ci  and  C-i  here  ate,  of  course,  determined  by  the 
apparently  arbitrary  choice  of  'constant*  terms  in  Sq  at'd  Sq  that  was  made  in 
(296),  but  we  are  now  close  to  seeing  why  the  particular  choice  of  coefficients 
was  made. 

We  are  left  with  only  the  last  term  of  (314).  MultipJying  the  expansion 
for  62  ,  viz  (226) ,  by  -a  cos  v  ,  we  obtain  for  this  last  terra,  quoting  in 
s]celeton  form  only, 

-ae,  cos  V  -  — i —  ae"^  [  f  ^  {  9e^r  +  ...+  9e^r  ,} 

2  18432  '  n  -3' 

+  4f{6he^Cg  +  ...  +  6he''c_^} 

+  2  { e^  (8h^  +  9f^)  cos  7 V  +  ...  -  24e[2  (88h^  -  21f ^)  , 

(168,-88,-207)1  }]  . 

. (319) 

When  the  full  version  of  (319)  is  combined  with  the  long  preceding  expression,  a 
tremendous  sin^liflcation  occurs,  since  no  less  than  20  pairs  of  terms  cancel  out 
and  a  factor  16e  can  be  extracted  from  the  remaining  terms.  The  result,  which 
represents  the  last  three  terms  of  (314),  is  worth  giving  in  full,  being 
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— a  [f^{24e¥,  +  120e¥g  +  e[185,-17]  +  16[5,-8ir^  -  eti81,131]r3 

-  264e^r2  -  72e¥j} 

-  4f{l6e^  (7,-9)C5  +  48e^(10,-13>C^  +  et  (538,-717) ,  (86,-147)  IC3 

+  32 [(10, -13), -2 (6, -5)1  C3  -48e[12(l,-l),5(l,-l)]C3 

-  48e^(14,-15)Cp  -  16e^  (11, -12)C_j  } 

-  16  {e^  (8, 0,-15)  cos  3v  +  3e^  (16, -16, -11)  cos  2v 

+  36(16(2,-4,1), (8, 0,-15))  cos  v  +  [ (48, -144, 77) ,  8 (4, 0, -7) ) }  ]  . 


We  can  now  make  the  grand  combination  of  this  last  expression,  which 
represents  all  but  the  first  two  terms  of  (314),  with  (316),  which  (on 
multiplying  out  the  factor  1  +  e  cos  v  )  compatibly  represents  the  first  two 
terms.  In  the  process,  another  12  pairs  of  terms  cancel  out  and  a  common  factor 
can  be  extracted,  which  has  the  effect  of  changing  an  external  factor  from 
'a'  into  *p'  ;  the  final  result  is 


-  4f{e(38  -  51f)C3  -  32(6  -  7f)C^} 

+  16(16h^  -  13f^)]  .  (320) 


As  one  of  the  definitive  formulae  of  the  paper,  this  has  been  written  without  use 
of  the  bracket  conventions;  its  brevity,  in  comparison  with  the  individual 
components  of  (314),  is  striking,  and  it  is  likely  that  there  is  some  physical 
significance  in  this  compression. 

It  should  now  be  clear  how  the  coefficients  of  Fq  and  Cq  ,  together  with 
the  'pure  constant',  were  chosen  in  (226),  and  likewise  for  the  coefficients  of 
Eq  and  So  in  (296)  already  referred  to.  They  were  triginally  set  as  algebraic 
unknowns,  and  carried  through  the  complete  analysis  for  r2  .  The  five  unknowns 
were  then  derived,  by  solution  of  linear  algebraic  equations,  so  that  the 
coefficients  for  Fi,  F-i,  Ci,  C-j  and  cos  v  in  (320)  should  all  be  zero. 


TR  88068 


93 


C.S  9«s*  Mlwct  la  w  mA  « 

Aa  a|ipt««eh  siailar  to  tkot  a«o4  tec  tj  tioald  «l*o  bo  posoiblo  for  v;  , 
but  ia  this  COM  t)M  tflMct  ow>ceoeli«  eoly  sliybtJty  addd^CiM.  is  siaplor  tbo 
stoctiag  poiat  is  (42>.  oOieh  yivos  tOo  poxtiol  docivotivos  t,  osd  that 
occur  in  the  (icst~oxdsc  focaoio 

».  •  »,o.  •  *,>t.  ■  021) 

The  ccrcospoodia^  fonsolo  for  vj  r=«oiros  frvo  toracs.  dcrivM  ires 
Taylor  oxpoasion  as  toaarkod  la  soctlcc  €.2.  Out  .'irfforoeciatios  of  v,  and  vv 
with  rospoct  to  M  is  sot  iaoadtlato.  ufiosoas  «rta  rospo^rt  t;  t  it  is.  so  it 
is  profocablo  to  start  with  tho  follcwizsq  srs-tors  foraaila! 


-  " 


%*‘2 


34  (*, 


W,  M.w.  ) 

-Hv  ;  . 


here 


q*^  sin  V  (q*  cos  v  +  2o<2  ♦  a  cos  v)  ]  , 

-  q'* («  sin^v  -  cos  v  (2  +  a  cos  v> ]  , 
q  *  (1  +  a  cos  V)  [2q^  cos  v  +  3e(l  +  a  cos  v)  ) 


and 


Mv 


-  2eq  sin  v  (1  +  e  cos  v) 


(323) 

(324) 

(325) 


(326) 


(The  hybrid  nature  of  these  partial  derivatives  must  be  noted;  thus  v^  is  the 
e-derivative  of  v  with  M  held  constant,  but  v^g  is  the  e-derivative  of  Vg 
with  V  held  constant . ) 

It  is  convenient,  in  substituting  these  second-order  partial. -derivative 
expressions  into  (322) ,  to  split  their  components  into  three  combinations  of 
terms.  The  first  combination  is  given  by  the  first  half  of  each  of  Vgo,  Vg^ 
and  vne  ,  together  with  the  whole  of  v^y  ;  the  terms  involved  are  precisely 
the  ones  that  arise  when  the  differentiation  of  Vg  and  v^  ,  given  by  (42),  is 
restricted  to  the  variation  of  e  cos  v  ,  and  this  first  combination  of  termr 
factorizes  in  a  convenient  way.  The  second  combination  involves  the  second  half 
of  Vgg  and  Vgy  ,  and  the  third  combination  involves  the  second  half  of  v^g  . 
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N«  now  express  the  working  formula  for  /  corresponding  to  (314)  for  r2  ,  as 
a  sum  of  five  terms,  of  which  the  first,  third  and  fourth  represent  the  three 
combinations  of  partial -derivative  terms  just  described;  the  second  term  of  the 
'working  formula'  covers  the  ej  term  from  (322),  whilst  its  final  i jrm  covers 
the  M2  term.  Thus  we  get 

Vj  -  14?”^  cos  V  -  eVj  sin  v)  {e^  sin  v  +  2q’^Mj(l  +  e  cos  v)  } 

+  q  sin  v  (2  +  e  cos  v) 

+  )4q  ^6^  (2eej^  sin  v  +  q^v^  cos  v)  (2  +  e  cos  v) 

+  Vs  eq  +  e  cos  v)^  +  q'^M2(l  +  e  cos  v)^  .  (327) 

For  the  first  term  in  (327),  we  have,  from  (148)  and  (182), 

Sj  cos  V  -  eVj  sin  v  -  •^[f{3e^C^  +  leeC^  +  2[10,7)Cj  +  32eCj  +  lle^C^^} 

+  4h{e^  cos  2v  +  86  cos  V  +  (6,1)}]  , 

(328) 


whilst  from  (148)  and  (164)  we  have 

Oj  sift  V  +  2q''^Mj  (1  +  e  cos  v)  - 
-  e'^  [f{3e^Sg  +  30e^S.  +  10e(10,-l)S^  +  2(56,-5)83  ^  ‘’8eq^Sj 

-  6(8,9)Sj  -  €6(14,1)83  -  30e^8_j  -  36^8,3} 

+  4h{e^  sin  4v  +  lOe^  sin  3v  +  6e(6,-l)  sir.  2v  +  2(24, -7)  sin  v}  ]  . 


. (329) 

Multiplying  the  product  of  these  expressions  by  a  further  factor  )4q"^  ,  we  get 
the  required  'first  term';  the  result,  in  3)teleton  form,  is 
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-  +  8£h{3e^Sj  -  ...  -  7o®S_J 

+  2{e®<8h^  +  21f^)  »in  6v  +  ... 

+  4[16(72h^  +  25f^),2{168h^  +  457f^)  ,-(l04h^  +  33f^)J  sin  v}]  . 


For  the  second  term  of  (327) ,  we  multiply  out  the  ex.sansion  for  e;  ,  viz 
(226) ,  by  q“2  sin  v  (2  +  e  cos  v)  .  The  result,  in  skeleton  form  is 


e'^q’^  [f^{9a^Zj2  +  ...  -  9e^Z_^}  +  8f{3e^hS^g  +  ...  -  Se^hS^^} 


1 

36864 


+  3{c^(ah  +  9f^)  sin  8v  +  ... 


-  16[4(216, -648, 56j», 24(44, -124, 77),  (648,-312,-907))  sin  v}] 


Combination  of  the  full  versions  of  the  first  two  terms  of  (327)  now  q.ves  the 
following  complete  expression: 

■  slfe"  ■*  144e^E^j  +  6e^  [154, 19)2^^  +  24e^  (124, 61)1, 

+  2e(2408,3592,39)Zg  +  8  [392, 2040, 133) 

+  2e[7824, 2470, -859)2,  +  32(98,267,-314)2, 

0  S 

+  64e[46, -251, -68)2^  -  32[J8,2i5,41b)23  -  6e [688,  882,  343)23 
-  24(24, 192, -49]2j  -  6e [312, 184, -115) 2^  -  24e^[84,-ll)2_, 

-  6e^[138,-13)2  ,  -  144e’2  ,  -  9e^2  ,  } 

+  8f{3e^hSjg  +  48e'*hSj  +  e^[316h,-(98,-123>lSg 
+  4e^[268h,-(234,-291)  [S, 
t  e[1872h,-4{670,-807>,-(526,-717))S, 

n 

+  4l336h.-2(130,~91*,-(782,-llt5!)S,  + 
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+  e [16  <266, -391) , -2 (914, -1787),  -  (34 6,  -687)) 

+  48(28)1,3(62, -73), 12  (7,-6)383 

-  4e (8 (78,  -77) , - (5074,-5785) , - (422, -459) ) 

-  16(36)1, 2(28, -23), -(518, -625) 

+  3e[16  (26,  -15) , -2  (2066,-2235) ,  3  (42,-71) ) 

-  4  (144)1, 2  (766,  -969) ,  (2210,-2367)  ]  S_3 

-  6(1296)1,1080(6,-7;,  (1418,-1503)18.3 

-  46^(228)1,  (494,-561)18.3  -  [300)i,  (190, -213)  1 8 

-  486^)18.5  -  3e^)iS.g} 

+  2{aN8)i^  +  9f^)  (6  sin  8v  +  16  sin  7v) 

+  2e^(6(72, -216, 235), -<56,24, -171)  1  sin  6v 

+  86^(24  (16,-48, 49), -(152, 24,-465)  5  sin  5v 

+  2e [24 (120, -360,337) , -8 (280, -24,-925) , - (392, -24,-661) ) 

sin  4v 

+  8(24(24, -72, 61), -4(168, -24, -817), -{904, -504, -?211) 1  sin  3v 

+  2e  [16  (72, -216, 671), -2  (4808, -5592, -4277), -0  368, -840, -1215)1 

sin  2v 

+  32(6(24, -72, 61), -(360, -984, 89), -8(47, -39, -40)  1  sin  v}]  . 

For  t)>e  tliird  term  of  (327)  we  require  (to  be  multiplied  out  by  )4q"^  ei  ) 
t)ie  following  result,  derived  using  (148)  and  (182)  again; 
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(2  +  •  COS  V)  (2e  Sj  sin  v  +  cos  v) 

-  -  -1-  [£{3e^Cl,-3]S^  +  30e^[l,-3)Sj  +  a[10C,-283,-51JS5 

X92 

+  4  [28, -53, -8918^  +  2etl22,-287,-45]S3  +  64(1,3,-4]S2 
+  1001*6,35,1)83  -  4 (12, -37, -89) 8q  -  a[84,-251, -67) 

-  30e^[l,-3)S_2  -  36^(1,-3)8_3} 

+  4h{e^tl,-3)  sin  5v  +  10e^(l,-3)  sin  4v  +  3e(12,-35,-3)  sin  3v 

+  4(12,-31,-19)  sin  2v  +  2e(18, -85, -3)  sin  v}]  . 

. <330) 

Again,  the  fourth  torm  requires  3eq”^  Mj  <1  +  e  cos  v)^  to  be  multiplied  by  this 
same  factor,  Mq*^  ex  .  It  is  natural  to  combine  the  two  terms,  therefore,  the 
result  of  the  combination  being  expressible,  with  the  factor  )iq*^  included,  as: 


-  «‘V\  [f{  30^83  +  300^85  +  edOO,  53,  -36)85  +  4(28,85,-56)8, 

+  26(290,-119,-66)83  +  64(1,6,-7)82  -  2e[42,17,  4  618j 

-  4(12,53,-8)8g  -  6(84,37,-4)8.3  -  300^2  -  36^5.3} 

+  4h{e^  sin  5v  +  lOe^  sin  4v  +  3e(12,5,-4]  sin  3v 

+  4(12,23,-16)  sin  2v  +  2e[90,-49,-6)  sin  v}  ]  . 

On  substituting  for  ex  from  (148)  and  multiplying  out,  we  get,  in  skeleton 
form, 

•  "  8fh{3e^S3„  +  -  3e-^3.5} 

+  2{e®(8h^  +  9f*)  sin  8v  +  ...  +  16(12  (24h^  +  7f^), 

2(984h^  +  643f^),-8(116h^  -  43f^),-(184h^  +  427fS )  sin  v}]  . 
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This  leaves  just  the  last  term  of  (327) .  on  nultiplv  inr,  the  expansion  for 
M;  ,  vix  (296),  by  q~^  (1  e  cos  v)^  ,  we  obtain  for  this  last  term,  in 
skeleton  form. 


(1  +  8  cos  v)^ 

-  — i—  e"V^  [3f^{3e*(2,-l)Z,,  +  ...  -  36*(2,-l]r  J 

3SS64  ’ 

+  8f{3e®h(2,-l]SjQ  +  ...  -  3e*h[2,-lJS_g} 

+  2{e®(8h^  +  9f^)t2,-ll  sin  8v  +  ... 

+  16e (48 (16, -48, 61), -24 (0,64, -121), 

6(64, 16, -253), -(8, -312, 453))  sin  v}  ]  .  (331) 


On  combining  the  full  versions  of  (331)  and  the  preceding  expression  (that 
represents  the  third  and  fourth  terms  of  (327)  together) ,  a  common  factor 
can  be  extracted,  for  cancelling,  and  the  resulting  expression,  in  skeleton  form, 
is 


— e'V^  +  ...  -  9e®Z  .}  +  8f{3eSs,„  + 

36864  ^  ‘  12  I  10 


+  2{e*(8h^  4-  9f^)  sin  8v  +  ... 


3e®hS_g } 


+  32e(6(40, -120, 183), -(744, -1464, 307), -2(44, -60, 97))  sin  v}]  . 


The  grand  combination  of  the  full  version  of  this  expression,  which 
represents  the  last  three  terms  of  (327),  with  the  earlier  expression  (given  in 
full)  for  the  first  two  terms,  is  now  possible.  A  further  (common)  factor  of 
8q^  emerges  in  the  combination,  and  ten  terms  cancel  out  coitpletely.  The  final 
result  is 
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X  [£*{9«*Ijj  +  108«^I,  +  6#*t82,15)Ig  +  12et84,67]I^  +  [784,2448,19111^ 

+  16«tlS9,771Sg  +  4«*[482,25jr^  -  48e[9,-81£3  -  (144,624, -SlZ^ 

-  12e(36,19]rj  -  6**[66,-l)rQ  -  108*^£_3  - 
4f{6«SSg  +  72e^hS3  2e^Il‘^2h,-3(14,-17)  )Sg 

+  12«[64h,- (34, -39)185  +  (672h,-36(2,3),-(94,-lS3)lS^ 

+  8e[8(21, -31),  (130,-139)183  +  2e^  [32  (77, -92),  (230,-243)  1  S3 
+  32«[3(4,-3),  (79,-91)lSj  -  [288h,  16 (116, -147) , -(266, -303)  1 
-  12e[48h, (110,-129)  lS_j  -  6e^[52h,  (34, -39) 1 

-  72e^hS.3  " 

+  2{«^(8h^  +  9f^)  (e  sin  6v  +  12  sin  5v) 

+  6s^[2(40,-120,127),-(d,8,-37)l  sin  4v 
+  126(12(8,-24, 23), -(24, 24, -137)1  sin  3v 
+  [48(24, -72, 61), -96(4, 8, -41), -(312, -168, -515)1  sin  2v 

+  86(6(8,-24, 61), -(144, -96, -257)1  sinv}]. 

. (332) 


With  th6  formula  for  vj  to  hand,  th6  formula  for  uj  may  ba  obtainad  at 
onca,  simply  by  Incorporation  of  the  formula  for  «2  •  (2V3) .  Th6  sinyillfi- 

cation  is  now  dramatic,  sines  20  pairs  of  tarms  cancsl  in  ths  combination,  and  a 
factor  26^  can  bs  sxtractsd  from  ths  rsmaining  tarms.  On  tha  dsbit  sida,  it  is 
no  longar  ths  cass  that  all  ths  £-tsm>s  contain  as  a  factor,  whilst  tha 
S-tsrms  contain  f  ,  but  this  property  is  rsstorsd  whan  ws  finally  procaad  to  w^ 
in  section  6.4.  Ths  result  for  112  is  given  by 
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u,  -  -i—  r{«*(32,«,-43)r.  +  4* (48,0, -S3) X.  +  8 [2 (18, -3, -17),  (24,-15,-11)  ] I, 

2  2304  L  *  «  5  4 

*  12*(48,-32,-l»)X,  *  *^(288, -264, -11)  Xj} 

-  4 {2«* (78,-158,81)8^  ♦  4«(88, -198,119)83 

-  [8(48, -106,55), -(144,34, -153)]8j  -  16«(102,-440,359)Sj 

-  4«^(18,-19,l)Sj} 

+  4  {•’(120,-60,-97)  sin  2v  +  8s (156, -276, 71)  sin  v}]  .  (333) 

Th*  co«££ici«nts  of  Xj,  Sj  and  sin  jv  ,  for  |j|  82  ,  in  both  (332)  and 
(333)  ar«  colourad  by  the  choica  of  th«  'constant*  tarns  (for  ej  and  M^)  in 
(226)  and  (296),  that  was  made  to  suit  r2  ,  so  it  is  noteworthy  (and  gratifying) 
that  (333)  contains  no  terns  in  X^,  X.i,  X.2,  S-i  or  S-2  .  No  tsm  in  Xq 
occurs,  because  it  has  been  avoided  by  the  choice  of  £0  coefficient  (for  (i);)  in 
(273);  there  is  still  a  tern  in  80  ,  however,  l>ecause  the  coefficient,  in 
(273),  was  chosen  to  avoid  a  term  in  W2  ,  not  U2  . 

6 . 4  fare  sliost-peslod  perturbetleae  la  b  end  « 

Ne  now  have  expressions  for  the  pure  second-order  perturbations  in  r,  u,  i 
and  Q  ,  given  (in  the  form  of  r2,  U2,  i2  and  Qj  respectively)  by  (320), 
(333),  (237)  and  (259).  These  are  the  four  fundamental  quantities  of  Kozai^''  and 
it  remains,  in  pursuance  of  the  philosophy  of  section  3.3,  to  compress  U2,  i2 
and  Q2  formulae  for  b2  and  W2  ;  in  combination  with  r2cor  ^200 

**700  as  given  by  (307) -(309),  the  spherical  polar  representation  of  second-order 
short-period  perturbations  will  then  be  complete. 

Our  starting  point  is  (90),  which  resulted  from  the  identification  of  two 
formulae  for  (x  y  r)‘  ,  one  in  terms  of  osculating  elements  and  the  other  in 
terms  of  semi-mean  elements.  Ma)(lng  use  of  (87)  and  (91),  we  have,  for  the  three 
cosponents  of  (90), 


cos(b  +  5b)  coa(w  +  ^w)  -  cos  8Sl  cos  u  -  c  sin  5(i  sin  u  ,  (334) 

cos(B  +  5b)  sin  (w  +  5w)  -  c  sin  50  cos  u  +  {cos  5i  -  ccd  -  cos  50)}  sin  u 

. (335) 
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and 


8in(b  -f  jb)  "  -  s  sin  SQ  cos  u  +  (sin  jl  4-  cs(l  -  cos  iO)  ]  sin  u  .  (3^6) 


As  thosa  fonmilaa  include  the  long-period  carry-over  termr  of  the  short- 
period  perturbations/  we  have  (ignoring  Js  ,  of  course/  in  the  definition  of  the 


Si 


Ki.  + 


-2 

X 


(i,  + 


(337) 


(and  similarly  for  jfll ) ;  also 

5b  -  Kbj  +  (bj  +  bj^^)  (338) 

(and  similarly  for  5w  ) .  He  have  dealt  with  carry-over  terms  separately/ 
however/  so  we  can  interpret  (334)-(336)  as  applying  to  pure  short-period 
perturbations/  ignoring  ia.ip  in  (337)  and  b2<;o  in  (338).  We  also  replace  u 
in  the  three  formulae  by  u  4  Xui  4  . 

It  is  now  straightforward/  in  principle/  to  expand  each  side  of  the  three 
_2 

formulae  to  0(K  )  and  then  to  Identify  the  expressions  that  emerge  as 
coefficients  of  X^  .  (Identification  of  the  coefficients  of  K  merely  leads 
again  to  the  formulae  (96)  and  (94)  that  were  used  in  the  generation  of  bi  end 
**1  f  9lven  by  (189)  and  (190).)  He  naturally  malce  use  of  the  generic  second- 
order  formula 


8ln(8g  4  xOj  4  X*9j)  -  sin  9^  (l  -  HiX^ej^  )  4  cos  (kOj  4  x^e^)  (339) 


and  the  corresponding  cosine  formula. 

On  this  basis  we  derive  from  (336) ,  since  it  is  now  permissible  to  omit 
Isars  from  all  the  terms/ 

bj  ■  sin  u  ij  -  s  cos  u  Qj  +  ^  "  c  sin  u  4  Uj  (cos  u  ij  4  a  sin  uflj)  .  (340) 
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The  first  two  teims  of  this  fonmila  are  given  by  (237)  and  (259),  ana  their 
combination  is  given  by 

sin  u  ij  -  s  cos  u  Oj 

-  -i-  3  c{2(l  -  f)  (e*  sin(Su  +  2v)  +  6e  sin(5u  +  v)  +  3  [3,2]  sin  5u 
576  ‘ 

*  18e  sin(4u  -f  (0)  -f  9e^  sin(3u  +  2(0)) 

+  (a^(6,-13)d5  +  4e(15,-29)<j^  +  6[  (21, -41) ,  3  (2, -7)] 

+  I2e(15, -53)02  +  54e^{l,-l)Oj) 

-  (36^(32,-41)33  +  8e(ll, -18)32  "  4(3(9, -11)  ,2  (27, -35)  Is^ 

-  120e(13,-22)s^  -  4e^  (21, -37)  s_3 )  }  ; 

. (341) 


the  coefficients  of  Oi,  si  and  s-j  here  are  determined  by  the  choice  of 
'constants'  made  in  (237)  and  (259),  these  choices  being  now  close  to 
explanation. 

The  first  five  terms  of  (341)  are  of  the  form  sin(jv  +  5(0)  ,  but  there  is 
no  point  in  having  a  special  notation  for  these  as  they  all  disappear  when  the 
third  term  of  (340)  is  combined  with  (341).  This  combination  reduces  to 

-i-  sc  {e^(32,-39)o.  +  8e(18,-25)0.  +  6[4  (6,-il) ,  (16, -31)  JO, 

576  04  J 

+  24e(6, -25)02  -  9e^  (20, -23) Sj  -  56e(2,-3)S2 
+  8(3(6, -7), 4(11, -13)  JSj  +  120e(14,-23)Sp  +  32e^  (3, -5)  s_j  }  . 


It  remains  to  evaluate  t  ^e  last  term  of  (340).  From  (149)  and  (151)  we  get 
cosuij  +  s  sinufl^  -  Ki  s  c  (2ec2  +  30^)  ,  (342) 

and  then  multiplication  by  (183)  yields  the  expression  required  O'-,  combining 
this  expression  with  the  preceding  expression  for  the  first  three  terms  of  (340), 
we  finally  derive  our  result  for  ba  .  As  with  rj  ,  given  by  (320),  we  express 
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this  without  usa  of  tha  brsckat  convantions.  Its  bravity  being  almost  as 
tanarkabla  as  for  X2  ;  thus. 


bj  -  -  a  c  [f{7a^Oj  +  40aO^  +  2(48  +  130^)03  +  360003} 

-  {e^(140  -  177f)s3  +  40e(8  -  9f>S3  +  8e(120  -  229f)SQ}]  , 


(343) 


(In  aquation  (43)  of  Raf  4,  the  factor  s  c  was  unfor'wunately  not  included. ) 

This  result  is  free  of  terns  in  Oi,  si  and  s-i  ,  which  would  arise  if 
the  right  choice  had  not  been  made  for  the  coefficient  of  Cq  and  the  'pure 
constant'  in  (237),  and  the  coefficient  of  So  in  (259).  In  practice,  of 
course-,  the  'right'  choice  was  made  by  setting  algebraic  unknowns  in  (237)  and 
(259);  these  ware  carried  through  the  analysis  for  bj  and  then  determined  such 
that  the  coefficients  of  Oi,  si  and  s.i  would  be  zero. 

To  derive  wj  ,  we  must  return  to  (334)  and  (335),  prepared  to  use  either 
or  both  of  these  equations  after  elimination  of  cos (5  +  5b)  or  substitution  for 
its  expansion  as  (1  -  >4K  bi^)  .  Again  recalling  that  the  carry-over  terms  have 
already  bean  dealt  with,  we  can  express  these  equations,  after  the  aforesaid 
substitution  and  an  expansion  of  the  right-hand  side,  as 

JOS  w  -  cos  u  -  K(c  sin  u  fl^) 

+  K  {sin  u  (si^flj  “  cQj)  -  )4  cos  u  ~  )  }  (344) 

and 


sin  w  -  sin  u  +  K(c  cos  u  flj) 

2 

+  K  {c  cos  u  fij  -  )4  sin  u  (i^*  +  -  b^)}  .  (345) 


It  is  now  convenient  to  write,  as  in  Refs  2  and  3, 

U  -  w  -  u  -  0(K)  ,  (346) 

since  (339)  and  the  corresponding  cosine  formula  then  immediately  give 
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and 


cos  w  “  cos  u  -  ICttj  sin  u  -  K  (tt^  sin  u  +  )4\)j  cos  u) 


*  2 
sin  w  »  sin  u  +  10)^  cos  u  +  K  (t)j  cos  u  -  sin  u)  , 


<347) 


(348) 


aquations  that  nay  ba  idantified  with  (344)  and  (345)  respactively.  In  both 
identifications  the  terms  in  K  siaply  lead  to  (94)  in  the  form 

t)j  -  cQj  ,  (349) 

^2 

but  the  coefficients  of  K  le».d  to  two  alternative  expressions  for  y>2  ,  given 
by 

(bj  -  cflj)  sin  u  -  -  sin  u  -  )4  (bj^  "  +  V^)  cos  u  (350) 

and 

(Vj  -  cOj)  cos  u  -  )S  -  c^Qj*  “  sin  «  •  (351) 

In  Refs  2  and  3  it  was  thought  best  to  develop  (350)  and  (351)  separately,  show¬ 
ing  how  the  same  formula  for  b2  -  cSi2  could  be  derived  from  either,  but  the 
most  natural  procedure  is  to  multiply  the  ecjuations  by  sin  u  and  cos  u 
respectively,  and  then  add  them  together.  When  this  is  done,  we  get,  with  the 
terms  in  bi^  (and  also  Vi^  )  cancelling  as  was  inevitable, 

\)j  -  cflj  “  -  )4{23<1  -  cos  2u)ijn^  +  sin  2u  (i^^  -  ftl^^)  }  ,  (352) 

since  it  is  now  legitimate  to  drop  all  bars. 

The  terms  on  the  right-hand  side  of  (352)  can  be  obtained  at  once  from 
(149)  and  (151) .  Thus,  direct  expansion  of  products  gives 
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ijQj  -  ~8(1  -  f)  +  SeZj  +  3t3,2ir^  +  ISeE^  + 

-  6e^S^  “  laeSj  +  6[3,-2]S2  +  36eSj  +  ISe^S^  -  36e  sin  v}  ,  (353) 


ij2  -  i  f(l  -  «){«^rj  +  esFj  +  3[3,2ir^  +  ISerj  +  +  12eC3  + 


36C, 


+  SeeCj  +  6*^  cos  2v  +  24e  cos  v  +  [27,10]  }  (354) 


and 


~  +  esfj  +  3t3,2)r^  +  iser^  +  oe^r^ 


-  12e^C^  -  36eC3  -  24e*C2  +  36eCj  +  36e^C^ 


+  30e^  cos  2v  -  24a  cos  v  -  [9,46] }  , 


(355) 


and  the  last  two  of  t)iasa  glva 


i.^  -  fQ,* 


—  £(1  -  ft  {e^r.  +  6er.  +  3[3,2]r.  +  IBeP,  +  9e^r, 

36  o  I)  %  J  K 


6a^C^  -  12aC3  +  6[3,-2]C3  +  36eCj  +  ISe^Cp 


+  18e^  cos  2v  +  9[1,-21  }  .  (356) 


It  Is  now  natural  to  coniblne  the  cos  2u  consonant  of  the  first  term  of  (352) 
with  the  whole  of  the  second  term,  the  result  )Hilng  (with  the  factor  -V* 
included) 

-  -i-  £(1  -  f)  (8e^S.  +  12eS-  -  24eV  -  36eS,  +  24e^  sin  2v  +  48e  sin  v)  , 

144  s  j  i  i 

after  which  the  full  evaluation  of  (352)  yields 

\)  -  O,  cos  i  -  — i-f  (1  -  f)  {eV  +  6er,  +  3[3,2]E.  +  18eS,  +  9e^Z,  +  2o's^ 

2  *  144  ft  J)  ^  J  2 

+  18[l,-2]Sj  +  ISe^Sg  +  24e^  sin  2v  +  12e  sin  v}  .  (357) 
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It  now  follows  from  (357)  and  (259)  that. 


»  -  -  -i-  <1  -  f)  {e^(4, 5)1^  +  246(1, 1>S.  +  3t2(6,5),  (8,3)]!^  +  24M3,1)Z, 

2  288  b  o  H  j 

+  36^(12,1)12  -  26^(39,-41)8^  -  166(11,-8)83 
+  12[(16, -27), -2(3,4)]S2  +  486(17, -30)8^  +  2e^(18,-l)Sg 

+  6e^(4,13)  sin  2v  +  48e(3,2)  sin  v}  . 


Finally,  from  (358),  (346)  and  (333)  it  follows  t)iat 


(358) 


i-  [f^{3«^E,  +  20er,  +  16(2  +  e^)£.  +  36er.  -  13e%} 

304  '■  6  5  4  ^  2  ^ 

+  4f{2e^(2  -  f)S^  -  4e(46  -  566)83 

-  £8(23  -  26f)  +  6^(14  -  39f))S2  +  166 (158  -  179f) Sj  } 

-  4{e^(72  -  168f +  59f^)  sin  2v  +  8e  (120  -  264f  +  95f'^)  sin  v}] 


. (359) 

As  with  (329)  and  (343) ,  this  has  been  deliberately  written  without  the  use  of 
the  bracket  conventions.  (In  Ref  4  the  wrong  sign  was  attached  to  all  the  4f 
terms. ) 

7  riMY -ORDER  ANALYSIS  YOR  J3  PERTURBATIONS 

7 . 1  Exact  aqpiatlons  for  rates  of  change  of  elements 

The  first-order  analysis  for  J3  is  very  similar  to  that  for 
(section  4),  the  main  differences  being  that  (as  for  all  the  odd  harmonics)  there 
are  no  secular  terms  in  the  solution  and  that  (as  for  all  with  f  >  2) 

there  are  long-pe.tiod  perturbations  Induced  by  .  Our  starting  point  is  the 
set  of  exact  equations  for  the  variation  of  the  osculating  elements,  given  by 
substituting  (49)  into  the  planetary  equations,  (22)-(27).  Working  with  (33)  and 
(34),  rather  than  (26)  and  (27),  to  keep  the  analysis  (inevitably  more  laborious 
than  for  )  as  single  as  possible,  we  get 
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i  -  “221  (£)*  {5f(7«y^  +  6T3  -  etj)  +  12g(5«Cj  +  2Cj  -  Sec^)  }  ,  (360> 

12q® 


4  -  •—  d/  {  5«  <7eT5  ^  ^^2  ■  •^1’ 

48q 

+  129(5003  +  12C3  +  68C3  -  4Cq  -  3ec_.)  }  ,  (361) 


i  .  “21  (£)'  iity  +  49C3)  , 


(36?) 


Q  -  ^  (f/  {5£03  +  (4  -  ISf)*^}, 


49- a 


(363) 


y  ”  d)^  {  5£(70O5  +  200^  +  8003  -  40^  +  eOj) 


4809' 


+  129(6083  +  I2S3  +8083  +  48^  +  308^3)  }  (364) 


and 


2Hns 


,2  'r 
3q 


(?)'  <5f<r3  +  12983)  . 


(365) 


Zntagratlon  of  thaaa  aquations  is  qultu  straightforward,  on  changing  the 
Ind070nd0nt  varlabla  from  t  to  v  (or  rathsr  v  )  as  in  section  4.  From  the 
integration  of  each  (  we  obtain,  in  the  notation  of  (72)  and  (73) ,  the  required 
expressions  for  (a.ip  and  (3  .  Results  for  (O,  O,  i  ,  and  hence  M  ,  are 
available  at  once  from  those  for  a,  Q,  y  and  p  ,  and  the  complete  set  of  (3 
and  (3, Ip  are  developed  in  sections  7.2  to  7.5.  Formulae  for  r3,  b3  and  W3 
are  obtained  in  section  7.7,  based  on  the  expressions  derived  for  V3  and  113 
in  section  7.6.  In  Part  3  we  ma)ce  use  of  the  Cl, Ip  our  study  of  the  long¬ 
term  variation  of  the  mean  orbital  elements,  as  perturbed  by  Jj  and  J3  ,  with 
singularity  very  much  to  ths  fore. 

7 . 2  Pectusbatlon  in  a 

As  in  sections  4.1  and  5.1,  the  best  way  to  evaluats  a3  is  by  the  special 
method  that  appeals  to  the  exact  energy  constant  a'  .  On  this  basis,  (52) 
leads,  for  the  potential  given  by  (49)  alone,  to 
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s/a*  “  1  +  -^  Ha<l  +  e  cos  v)*q“^  (SfOj  +  12gs^) 


(366) 


Thus 


*3  “  tSfOj  +  i2gSj)  . 


(367) 


By  ccnbinlng  successive  powers  of  p/r  (  -  1  -i-  e  cos  v  )  with  O3  and 
Si  ,  we  car.r  as  in  section  4.1,  write  down  progressively  more  complicated 
alteimative  expressions  for  ai  .  These  expressions,  all  of  them  necessarily 

(368) 

(369) 

*3  “  •;^aq‘*s(p/r)  {Sf  (e^Oj  +  6e^<i5  +  3et4.1)o^  +  4(2,3)03 

+  6e^Oj  *■  s^Op) 

+  12g(e^s^  +  6e^S3  +  3e[4,l)S2  +  4l2,3)Sj  +  36[4,1)Sq 

+  6e^3_j  +  e^s_3)  } 

. .-(370) 

and 

*3  “  {  5f  (e^O^  +  8e^Og  +  46^(6,1)03  +  2(8,24,3)03 

+  86(4,3)0^  +  4e^(6,l)Oj  +  Be^Oj  +  eVj) 

+  12g(eS5  +  8e^s^  +  46^(6,1)63  8«I<>3J»2  +  2(8,24,3)83 

+  8e(4,3)Sj3  +  4e^(6,l}s_3  +  8e^s_3  +  e^s_3)  }  . 

.  (371) 


sysnetrlc  in  both  their  O  coiqsonents  and  their  s  con^onents,  are 
Sj  -  "JJ**!’**  (p/r)^  {  Sf  (eo^  +  2O3  +  eo^)  +  12g(e33  +  283  +  es^)  }  , 

63  “  ^aq'^s  (p/r)^  {  5f  (e^Oj  ■(  4eO^  +  2(2,1)03  * 

+  12g(e^S3  +  4eaj  +  2(2,1)83  +  4eSjj  +  e^s_3)  }  , 
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W«  shall  lisa  (369)  for  th«  devslopnent  of  M3  via  J3  ,  (3^0)  for  the  develop¬ 
ment  of  rs  ,  and  (371)  for  the  development  of  P3  . 

The  final  expression,  (371),  can  also  of  course  be  developed  from  (360)  by 
the  general  method  now  to  be  used  for  03,  13,  etc. 

7 . 3  VcrtudhatlOBs  in  e,  1,  p  and  pe^ 

From  (361),  with  the  independent  variable  changed  to  v  by  use  of  (38),  we 
get 

^  “  Tm  «a{5f  (7e\7  +  48e^76  +  4et27,81y^  +  16[5,9)t^  +  42e[4,l]y^ 

+  16[1,6]y  +  4e[3,4]Y  -  e\  ,  ) 

2  1 

+  12g(5e^C5  +  32e^c^  +  4eC17,4)C3  +  16(3, 4)03  +  14e[4,l]Cj 

-  leq^Cg  -  28ec_3  -  ISe^c^j  -  36^0,3)  }  . 

. (372) 

The  term  in  cq  leads  to  the  long-period  variation  specified  by 

*3,<P  “  "  “  '  <273) 

whilst  the  rest  of  (372)  yields,  on  integration, 

®3  ”  lir  ®  (se^o,  +  40e^O5  +  4e[27,8)05  +  20[5,9]O^  +  70e[4,l]a3 

+  40(1,6)02  20e(3,4]Oj  +  406^0^  +  5e^O_j) 

+  4g(3e^Sj  +  24e^8^  +  4eI17,4)S3  +  24(3, 4)S2  +  42e(4,l]Sj 

+  24(1,6)S(,  +  84es.j  +  24e^a.2  +  30^3.3)  } ; 

. (374) 

the  coefficients  of  Oq  and  sq  in  (374)  have  l>een  chosen  to  ma)ce  rj  ,  given 
by  (408),  as  simple  as  possible. 
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Fxon  (362) ,  in  th«  suw  way,  w«'  get 

^  ^  He  {  5f  +  2t2,l)tj  +  4«|j  +  e^^) 

+  4g{e*C3  +  4eCj  +  2[2,l)e.  +  4eC(j  +  e^c_j)  }  .  (375) 


The  term  in  cq  leads  to 

i,  Ip  -  egc  cos  <0  ,  (376) 

whilst  the  rest  o£  (375)  yields 

ij  -  e{f  (ae^Oj  +  15e<T^  +  10(2,1)03  +  SOeO^  +  156^03) 

+  4g(e^S3  +  OeSj  +  6t2,llSj  -  3e^s_j)  -  66(4,-15)8^}  .  (377) 


The  coefficient  of  sq  in  (377)  has  been  chosen  to  malce  ba  ,  given  by  (410),  as 
single  as  possible;  the  fact  that  it  is  a  multiple  of  4  -  ISf  ,  rather  than 
4  -  Sf  as  for  the  other  s  coefficients,  is  bound  up  with  the  fact  that  in 
Q3  ,  to  be  given  by  (385),  the  coefficient  of  cq  is  a  multiple  of  4  -  5f 
whilst  the  other  c  coefficients  are  multiples  of  4  -  15f  .  (See  also  Part  2.) 

From  (373),  using  (28),  we  get 

Pj  lp  -  2pegs  cos  (0  ,  (378) 


and  from  (371)  and  (374)  we  get 

P3  -  P8{f  (3e*Oj  +  15eO^  +  10(2, IjOj  +  aOeO^  +  15e^Oj) 

+  4g(e^83  +  SeSj  +  6(2,1)83  +  18eSg  -  Se^s  ^)  }  .  (379) 
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Mu  now  haw*  our  uaual  chock  in  toraa  of  tho  quantity  pe^  that  ia  an  abaoluto 
conatant  of  tho  aotion.  Iquationa  (376)  and  (378)  yiold 

<Pc^)3.<p  -  0  .  (380) 

whilat  (377)  and  (379)  yiold 

(pc*)3  -  Mpoac^(8  -  15f)aj,  ;  (381) 

(381)  ia  'conatant'  (to  ordor  rero),  which  conpletas  tho  chock,  non-zoro  becauae 
tho  conatanta  in  03  and  i:  aro  not  choaen  by  concomitant  criteria. 

Uainq  (28)  in  rovorao,  wo  can  got  an  exproaaion  for  03  that  is  aomowhat 
moro  compact  than  (374)  but,  liko  (173),  la  unfortunate  in  carrying  tho  factor 
o~^  ;  tho  oxproaalon,  based  on  (367)  and  (379),  is 

O3  -  ^  o'^o{4(p/r)'(5f03  +  lags^)  -  q^f  (30*05  +  l5oO^  +  10(2,1)0, 

+  30oOj  +  15e*Oj) 

-  4q*g  (0*03  +  6OS3  +  6[2,l)Sj  +  18oSp  -  3e*s_j)  }  . 

. (382) 

7 . 4  PortudSationa  in  fi,  y  and  (O 

From  (363)  wa  got,  after  tho  change  of  integration  variable, 

~  Hca"^  {sf  (0*05  +  4oo^  +  2(2,1103  +  40Oj  +  0*03) 

+  (4  -  15f)  (0*03  +  400^  +  2[2,l)0j  +  403p  +  0*3,3)  }  • 
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Th*  t«t«  In  «o  iMdc  to 


Q-  ,  -  lfc«(4  -  15f)c*  sin  •  , 

3f«P 


whilst  the  rest  of  (3B3)  yields 


Qj  -  ~  {^(aeSs  +  ISey^  +  lOta.llYj  +  30eTj  +  ISey^) 

+  (4  -  15f)  (e^Cj  +  SeCj  +  6[2,l)Cj  -  3e^c_j)  -  24090^} 


The  coefficient  of  cc  is  chosen,  in  conjunction  with  the  coefficient  of  sq  in 
(377),  to  make  ba  as  simple  as  possible;  we  have  already  remarked  on  the 
effective  swapping  of  the  coefficients  4  -  5f  and  4  -  15f  ,  as  between  (377) 
and  (385) ,  and  the  full  explanation  for  this  (in  the  context  of  the  general  ) 
is  reserved  for  Part  2. 

From  (364)  we  get,  after  the  change  of  integration  variable. 


^  -  -i-  He'^sfSf  <7eV  +  48eV  +  2e[54,ll)o.  +  16[5,6]o,  +  24e(4,ljo 
dv  192  \  I  t  D  4  3 


-  16[l,-3)Oj  -  2e(6,-5)Oj  +  e^(T_^) 

+  12g(5e^S5  +  32e^8^  +  2e(34,9)Sj  +  16[3,§)Sj,  +  24eL4,l]Sj 


+  16[1,4)8q  +  14e[2,l)s.j  +  IBe^s.^  +  Se^Vj)  } 


The  term  in  sq  leads  to 


Vj  Ip  “  e'^gs  (1  +  4e*)  sin  u  , 


whilst  the  rest  of  (386)  yields 


■  \  y. 

'>'P 

' 'i 
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p.  '■  r  ‘  • 

'4 1 


•  M'  ' 


X13 


f 

I 


i 


i 

\ 


I 

! 

! 

Vj  -  -  ^  (5«=V,+  2«£S4,11Jy^+  20t5,«J7^  +  46«t4,lly^  j 

-  40tl.-3JYj  -  10«[6,-5}Yj  -  40«Sj|  -  5«\j  ) 

+  4g(u\  +  24«*c^  ♦  2«I34,9}c3  +  24t3,51Cj  +  72e[4,lJCj 

i 

-  24c^  -  42«[2,l)c_j  -  24e*c_2  -  3e^c_3)  }; 

. (388) 


^h«  co«fficl«nts  of  Yo  and  cq  in  (388)  have  been  chosen  to  make  W3  ,  given 
)3y  (413)  r  as  simple  as  possible. 

Expressions  for  <U3,{p  and  U3  are  nou  lamadlate/  by  subtraction  of 
cl^3,4p  Va.fp  and  of  CO3  from  ^3  .  These  expressions,  which  are 

necessarily  more  conqpllcated  than  for  Va.tp  and  ^3  because  differing 
combinations  of  f  are  involved,  are 

•  “  T  [<^9  -  e^(4  -  35f  +  35f*)  ]  sin  U  (389) 

A 


and 


<“■1  "  -  (Se’fY,  +  40e*fT  +  2e[54f,-(6,-17))  y  +  20t5f,-3(l,-3)  )y^ 

-  40e(2(l,-3),(l,-2)Jy^  -  40(f,3(l,-2))Y^ 

-  10et6f,  (6,-11)  ly^  -  40e^fy^  -  5e^fy_^) 

+  2(6e^fgcj  +  48e*f9c^  +  etl36f9,-(8,-74,75)  JCj 

+  12{12fg,-(4,-39,40)]C2  -  12et2(4, -43,45),  (4,-31, 30)]Cj 

-  48g(f,-(l,-l)lCp  -  3eI56f9,-(8,-66,65)  Jc,j  -  48e*fgc.2 

-  6e^fgc.3)  } . 
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7 . 5  PaxtuidMtloaa  la  p,  a,  J,  ll  aad  L 

From  (365)  w«  gmt,  after  the  change  of  variable, 

Hqs{5f  (e^Oj  +  4eO^  ♦  2(2,l]Oj  +  aeO^  +  e^Oj) 

+  12g(e^a,  +  403^  +  2(2,  l)Sj  +  4e8g  +  e^s.j)  }  .  (391) 

The  ter-n  In  sq  leada  to 

P3  ,lp  “  Seqsg  sin  U  ,  (392) 

whilst  the  res.t  of  (391)  yields 

P3  ”  - -^qslsf  (30^7^  +  1567^  +  10(2,1)73  +  30073  +  ISe^^) 

+  12g(56^C3  +  30ec3  +  30(2,l)c.  +  laec^  -  15e^c_^)  };  (393) 

the  coefficients  of  70  and  cq  (the  former  (seing  zero)  have  been  chosen  to 
produce  the  appropriate  coefficients  (of  7o  and  cq  )  in  M3  given  by  (400) . 

Ebcpresslons  for  03, fp  and  O3  are  obtained,  in  view  of  (32),  by 
subtraction  of  qy3,<p  from  P3,tp  and  efifj  from  P3  ,  the  results  being 

O3  -  -  a'^qsg  (1  -  4e^)  sin  o)  (394) 

and 

O3  -  {  5f  (lSe\+  120e^7g  +  6e(54, -5)75  +  60(5, -2)7^ 

-  40e(4,5l7^  -  120(1,5)7^  -  30e(6,ll)7^  -  120e^7jj  -  15e^7.^) 

+  12g(l5e^05  +  120e^c^  +  10e(34,-7)c3  +  SeOq^c^ 

-  120e(4,5]C3  -  24(5,16)C(,  -  30e(14,-9)c.j  -  120e^c.3 

-  150^3)  }. 

. (395) 
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>'oc  ^3  **  pgof,— d  as  la  aaotioa  A.2,  using  tbs  aquation 

oosca^poading  to  <1S9)«  togathas  irith  (3C9)«  to  darlaa 

.  -X  qs  {5f(a‘o.  ♦  4a«.  ♦  2[2,l)e.  +  4ao,  +  a*a,) 
vM  l§  *  *  '  ^  * 

+  \2qW»^  ♦  408^  *■  2[2,l)s^  +  408^  +  a**_j)  }  .  (396) 


Zn  vlaw  of  (40),  tra  now  gat 

Ij  jp  -  -  Soqsg  sin  (a  (397) 

and  (on  integrating  t)ie  remaining  terms) 

^3  “  MO  (30*75  *  *  1012,1)73  +  30e7j  +  15e27i) 

+  l2g(Se*C3  +  lOeCj  +  30(2, l]Cj  +  12005  -  15e*c_j)  }  .  (398) 

These  results  could  have  been  obtained  directly  from  (392)  and  (393),  since 

-  T  P.  ,  and  I.  -  4*  P<  ixom  the  general  formula  (quoted  in 

3««P  •  3a«p  3  *3 

section  4.2)'  with  <  set  to  3;  the  coefficient  of  co  in  (398)  has  been  ta)cen 

to  conform  with  this  relation  between  P3  and  /3  .  (Taking  different  constants 
in  (3  and  *4*^3  *'°**1^  provide  greater  flexibility  than  we  need;  this  was 
tacitly  assumed  also  for  and  '^^Pj  section  4.) 

The  formulae  for  a  and  J  can  now  be  combined,  the  results  oi.'  the 
combination  being  that 

Mj  ,p  -  -  e"*q*sg  sin  a  (399) 

and 

M3  -  e'^qsff  (iSe^r,  +  120eV5  +  6e(54,l)75  +  60[5,1)7,  +  80eq^3 

-  120(1,2)7^  -  30ef6,5)7j  -  120e*7o  -  15e*7-i) 

+  12g(3e*C5  +  24e*c^  +  2e(34,-l)c3  +  72Cj  +  48eq*c. 

-  24C1,2)C5  -  6e(14,-3]c.^  -  24e*c.j  -  3e*o_3)  }  . 


TK  II0(« 


(400) 
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Tb*  perturbation  in  L  need  n<^t  detain  ua  long.  Thus,  the  ioraulae  for 
Ii3,{p  and  L3  follow  at  once  frt*  (.hose  for  p3,tp  and  ps  (or  {3,1^  and  /3  ) 
since  (from  the  general  fomula  of  section  4.2,  with  t  set  to  3) 


“3.«P 


i 


-  -  -S-  / 

3  '3.«p 


••  iaqsg  sin  «> 


(401) 


and 


L 


^  a  B 

a  ^3 


(402) 


there  (seing  little  point  in  quoting  I13  in  full. 

7 .  C  Short-period  perturbationa  in  w  and  u 

It  is  recalled,  from  (74),  that  the  (first-order)  J3  contribution  to 
is  given  by  H((|3  -f  Ci.tpn)  •  Also,  S-v  and  6u  may  be  expressed  in  terms  of 
by  means  of  (42) .  Thus,  from  03,  <1)3  and  M3  we  may  derive  expresslor.s  for 
V3  and  U3  ,  whilst  from  e3,|p,  tt3,«p  and  M3,jtp  w«.  derive  expressions  for 
the  quantities  that,  with  an  obvious  notation  referring  to  carry-over  effects  (cf 
section  6.1),  we  may  denote  by  v3co  and  03^0  >  We  have  no  need  for  a  quantity 
V3,(p  >  such  that  V3co  *■  V3,{p  m  ,  because  the  long-period  perturbations  are 
applied  (in  a  non-singular  manner)  to  the  orbital  elements,  and  not  to  derived 
quantities . 

From  (42)  it  follows  that  ^v  can  be  developed  as  the  sum  of  three 
components : 

(A)  2q"^  (S«  sin  v  +  eq'^  cos  v)  ; 

(B)  )4oq"*  (fe  sin  2v  +  eq"^  SH  [cos  2v  +  3) ) ;  and 

(C)  q'^  8V>  . 

Ne  derive  the  three  components  separately  for  both  the  1^3  and  the  C3,tp  ,  and 
then  combine. 

In  combining  the  appropriate  multiples  of  63  sin  v  and  M3  cos  v  to  form 
Component  A,  eight  pairs  of)  terms  cancel,  the  result  being 

- q’*S  [  f  { ISe^Tj  +  eOe^Tj  +  e(28,47)y^  -  90q^3  -  5e[28,-13]T3 

-  60e^7i  "  15e\o)  + 


* 
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+  12g{3«^c^  +  12«*e3  -  •t4,-19JCj  -  24q*Cj  -  3*[4,l)Cj 

-  12.*c.j  -  . 

Rather  not*  t«tna  cancel,  in  a  sinilar  nanner,  in  forming  Coapenent  B,  the  result 
being 

^  eq‘*8  [f{30*^j  +  4*t3?,-7lT^+ 15tl3,-3]r^  +  SSoq^Yj  "  15[1,9]72 

-  2  0otl,5]ri  -  30e\} 

+  24g{3*^c^  ♦  6*I3,-1JC3  +  3[8,-3}C3  +  ISeq^Cj  -  ISe^c^j 

-  6e[3,-l]c,j  -  36^2 }]  • 

Four  pairs  of  terms  cancel  in  the  combination  of  Components  A  and  B,  and  a  factor 
q^  emerges  to  cancel  the  q~^  ,  the  result  being 

^  s  [f  {25*^5  +  139*T^  +  12 [15,8)73  +  ^SSey^  +  lOOe^,  } 

+  48g{3*^C3  *  *  4t3,2)Cj  +  Sec^  -  3e*c,j}]  . 


Component  C  is  given  at  once  by  (400),  of  course.  It  then  follows,  from  the 
overall  combination,  that 

V3  -  e'^8  [f {15*^73  +  120e^7j  +  2*1162,31)75  +  2[150,169]74 

+  8*[55,14)73  -  10(12,-29)7j  -  10e[18,-5)7,  -  120e\ 

-  ISeYj} 

+  12g{3*^C5  +  24**c^  +  2*(34,11)C3  +  8[9,14)C3  +  leotO.Dc^ 

-  24C5  -  6*(14,l)c.j  -  24*^c.3  -  3*^3}  3  . 

. (403) 


ta  ssocs 


FroM  cj.Cp  Na.tp  ,  It  Follows  in  tho  saas  way  that  Component  A  o£  th« 

long-poriod  carry-over  is  -  2  a  g « si  ,  Coapononc  B  is  -  M  s  s  g  m  (S2  3so)  ,  and 
Coaponsnt  C  is  -s'^q^sgaso  .  Putting  thsss  together,  we  get 

Vj^^  -  -  »  e"*  8  ga  (e*Sj  +  <0  8^  +  [2,llSj,)  .  (404) 

For  ju  ,  we  8iaply  feed  in  the  effecto  of  (Oa  and  tt3,<p  .  Thus 

^  s'^  {  f  (2e^<9,-10)Y5  +  e(90,-101>Y^  +  4[5(6,-7),  <15,-16>]Y3 

+  5e(36,-43>Yj  +  lOe^  (9, -14) ) 

+  6  (e^(4, -15,10)03  26(12,-61,50)03 

+  2(6(4,-31,30),  (12,-85,80)  ICj  -  6e(4,-9,5)OQ 

-  3e^(4,-31,30)o_j)  } 
. (405) 

“sco  “  “  i"  +  lefgSj  +  e(8,-66,65)8g]  .  (406) 

7 . 7  Short-period  perturbations  in  r,  b  and  w 

We  oan  now  ooioplate  the  J3  analysis  formally,  by  derivation  of  the 
quantities  required  in  the  spherical-coordinate  representation,  viz  ra,  ba  and 
W3  (pure  short-period  terms),  together  with  raco,  baco  and  waco  (long-period 
carry-over  terms) . 

From  (95)  it  follows  that 

ra  ••  (r/a)  aa  -  (a  cos  v)  03  +  (aeq*^  sin  v)  M3  .  (407) 

Eight  (pairs  of)  terms  cancel  in  the  combination  of  the  terms  in  ea  and  M3  , 
the  combination  being 
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-  as  {f  (iSa^O^  +  90a*O5  +  a(176,49]0^  +  15(7,13)«3  +  Sat32,13]02 

+  90a^Oj  +  ISa^Oj) 

+  12g(3e^a^  +  IBe^Sj  +  •(32,13182  +  12(2, 3]8j  +  3«(16,-1)Sq 

+  ISe^s.j  +  }  . 


On  conbinlng,  further,  with  (r/alaa  ,  given  by  (370),  nine  more  tenea  disappear, 
with  the  terms  in  si,  s-i,  <Ti  and  o.i  vanishing  as  a  result  of  the  choice  of 
constants  in  (374)  and  (400),  and  a  factor  can  be  extracted  from  the  rest. 
The  result  is 

r^  -  ps{  f  (4eO^  +  ISdj  +  20e<T2)  +  48eg(s2  -  3S(,)  }  .  (408) 

From  83, Ip  and  Mj^ip  ,  given  by  (373)  and  (399),  we  have,  similarly, 

*3co  "  psgm  cos  u  .  (409) 

Next  we  require  ba  ,  and  it  follows  from  (96)  that 

ba  -  (sin  u)  ia  -  (s  cos  u)  Q3  .  (410) 

Condblnation  of  the  f-terms  in  (377)  and  (385)  is  iimediately  possible,  giving 
cf  (3e*C^  +  ISeCj  +  10(2,1)02  +  30eCj  +  ISe^C^)  . 


The  combination  of  the  g-terms  in  (377)  with  the  (4  -  15f)  term  in  (385)  is 
only  a  little  more  involved,  and  leads  to  two  sets  of  terms,  vix 

-  ~  c  (2  -  5f)  (e^  cos  2v  -  3e  cos  v  -  3(2,1]) 


V  v* 
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-  of  +  6^3  +  6t2,l]Cj  -  3«*Cq)  . 

Lastly,  tha  chosan  Sq  tarn  in  (377)  conl^inas  alth  tha  chosan  cq  tern  in  (385) 
to  giva 

-  -i-  ca  [siC^  +  2(2  -  5£)  cos  v]  , 

pracisaly  tha  combination  raquirad  to  cancal  tha  terms  in  Ci  and  cos  v 
already  generated.  Combination  of  tha  four  sets  of  terms  yields 

ba  -  -  c  {f  (2e^C^  +  ISeCa  +  20  [2  +  e^lc^  -  SOe^C^j) 

4  4(2  -  bf)  (e^  cos  2v  -  3  [2  +  e^] )  }  .  (411) 

From  ia.ip  and  ils.tp  ,  given  by  (376)  and  (334),  we  have,  similarly, 

bjco  “  -J-  ecm  [sfSj  4  2(2  -  5f)  sin  v]  .  (412) 

Finally,  we  requite  W3  and  w3eo  .  Adding  U3  ,  given  by  (405),  to 

CQ3  ,  with  fl3  given  by  (385) ,  we  have 

’'3  “  ■  m  ^  4  e^lYj  4  3567^  4  SOe^Yj ) 

-  24g(e^C3  "  "  2[l8  4  Ye^jc^  4  9e^c_j )  }.  (413) 

Equation  (413)  may  also  be  obtained  by  direct  addition  of  (403)  and  (388),  » 

representing  V3  and  y3  respectively.  Yet  another  derivation  is  possible 
without  using  either  U3  or  V3  ;  we  return  to  the  sum  of  the  'A'  and  'B' 
components  of  V3  derived  in  section  7.6;  if  to  this  sum  we  add  L3  ,  as 

expressed  by  (402)  via  P3  or  J3  ,  instead  of  q'^  M3  ,  then  we  have  W3  at 
once.  From  U3co  and  003, |p  ,  or  from  V3C0  and  yfj.tp  >  ot  most  directly  from 
tha  'A*  and  'B'  components  of  together  with  q'l  La.tp  ,  we  obtain, 

similarly. 
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-  -  )lsgn(*S2  +  4Sj  -  7«8p)  .  <414> 

It  will  be  obaetved,  froa  (408) >  (411)  and  (413),  that  rs,  03  and  W3 
contain  five,  aix  and  nine  terms  respectively.  Similarly,  from  (188) -(190),  ri, 
bi  and  wj  contain  two,  two  and  three  terms  respectively.  As  an  introduction 
to  Part  2,  we  now  give  the  number  of  terms  associated  with  the  general  , 
noting  that  the  quantities  generically  denoted  by  Cl  the  present  Report  will 
have  to  )3e  relabelled  C2  i''  ^Art  2,  where  the  analysis  is  taken  to  'first 
order*  only,  but  covers  every  t  .  Odd  and  even  values  of  8  have  to  be  dealt 
with  separately. 

For  odd  8  ,  the  number  of  terms  is  8(8  -  V2)  > 

8(8  -  +  V2  for  b|  ,  and  8^  for  \>i  .  For  even  8  ,  the  number  is 

8(8  -  ’a)  +  i  for  both  r^  and  b<  ,  and  8^-1  for  wj  .  That  the  last 

figure  is  8^  -  1  ,  rather  than  8  ,  is  somewhat  fortuitous,  being  due  to  the 
fact  that  (as  we  shall  see  in  Part  2)  the  coefficient  of  ,  ie  of 

ain{  (8  +  l)v  +  2<a],  which  arises  only  when  8  is  even,  always  vanishes.  Thus 
the  absence  of  an  S3  term  in  (190) ,  which  seemed  strange  when  the  analysis  was 
originally  performed,  was  not  a  fluke  occurring  just  for  J2  . 

8  COKCIitlSZOKS 

The  main  function  of  this  Report,  Part  1  of  a  projected  trilogy,  has  been 
to  provide  full  details  of  the  untruncated  second-order  orbital  theory.  Involving 
J2  and  J3  only,  of  which  the  outline  was  given  in  Ref  4 .  The  principal 
novelty  of  this  theory  lies  in  the  reduction  of  second-order  J2“i>iduced  pertur¬ 
bations  to  very  contact  expressions  in  a  special  system  of  spherical-polar  co¬ 
ordinates  based  on  a  mean  orbital  plane.  First-order  J3-lnduced  perturbations, 
regarded  as  second-order  in  the  overall  theory,  have  been  expressed  in  the  same 
way,  the  results  being  particular  cases  of  formulae  (for  the  general  )  that 
will  be  developed  in  Part  2.  (Postscript:  now  available  as  TR  89022.) 

The  special  coordinate  system  cannot  be  used  in  the  treatment  of  secular 
and  long-period  perturbations,  but  this  is  of  no  consequence  as  so  few  torms  are 
present  in  the  basic  expressions  for  perturbations  in  the  standard  elliptic 
elements.  The  J3  perturbations  in  these  elements  suffer  from  singularities, 
however,  and  the  Report  has  included  introductory  material  on  the  treatment  of 
these  singularities.  Additional  formulae  relating  to  the  long-term  evolution  of 
an  orbit,  subject  to  J2  and  J3  ,  will  appear  in  Part  3  of  the  trilogy,  which 
will  also  Include  numerical  results  from  an  assessment  of  the  overall  accuracy  of 
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th*  ■•cond-ordct  inodAl  for  motion  under  J2  end  J3  .  Over  long  periods  of 
time,  such  that  the  total  angle  described  by  the  satellite  (in  true  or  mean 
anomaly)  approaches  in  order  of  magnitude,  the  errors  in  the  model 

Inevitably  Increase  from  third  order  to  second  order,  ie  typically  (for  close- 
Barth  sativllites)  from  centimetres  to  decametres. 

The  trilogy  of  Reports  is  to  be  viewed  as  part  of  a  continuing  study.  The 
current  limitations  of  the  theory  (for  Earth  satellites)  may  be  summarized  as 
follows:  only  perturbations  due  to  the  geopotential  have  been  considered  (though 
lunieolar  perturbations  were  addressed  in  a  precursor^  to  the  present  theory) ; 
tesseral  harmonics  have  not  been  covered  (though  the  effects  of  J2,2  on  a  low- 
eccentricity  orbit  were  considered  in  Ref  S);  the  zonal  harmonics  have  been 
restricted  to  J2  *nd  ,  except  in  Part  2  of  the  trilogy;  the  analysis  has 
been  taken  to  second  order  only,  implying  only  first  order  in  J3  ;  and  the 
build-up  of  long-term  error,  as  indicated  in  the  last  paragraph,  cannot  be 
avoided.  In  spite  of  these  limitations,  however,  the  present  Report  serves  as  a 
significant  step  towards  the  goal  of  an  analytic  (or  semi -analytic)  orbit 
generator  that  is  much  more  efficient  than  numerical  generators  of  con^arable 
accuracy . 

To  facilitate  a  reference  back  to  the  formulae  that  summarize  the  formulae 
developed  for  the  generator  coveting  J2  and  J3  only,  we  conclude  the  Report 
by  recording  the  relevant  equation  numbers.  The  modified  version  of  Kepler's 
third  law  is  given  by  (299) .  Secular  variation  rates  ate  given  by  (150)  and 
(297)  for  Q  ,  and  (152)  and  (298)  for  O)  .  Long-period  variation  rates  ate 
given  by  (300) -(306)  for  J2*  effects;  for  J3  effects,  the  equations  are  (373), 
(376),  (384),  (387),  (389),  (392),  (394),  (397),  (399)  and  (401).  The  carry-over 
effects  on  coordinates  are  given  by  (307) -(309)  for  J2^  •  and  (409),  (412)  and 
(414)  for  J3  .  Last,  but  far  from  least,  the  short-period  effects  on  coordi¬ 
nates  ate  given  for  J2  (first  order)  by  (188)-(190);  for  J2^  ,  by  (320),  (343) 
and  (359);  and  for  J3  ,  by  (408),  (411)  and  (413). 
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latltuda-lllca  coordinate  of  (r,  b,  w)  ayatem 
coa  i  (alao:  obaolate  cylindrical  coordinate) 
coa(jv  +  (0) 
coa(jv  +  2(0) 

J  coa  kS  dr 

°  t 

defined  auch  that  n'  +  dn*  conatitutea  Jtec 

used  in  context  of  J2^  contributions  to  dl^/dv 

eccentricity 

eccentric  anomaly 

sin^i 

Merson's  functions,  defined  in  equations  (122) -(124) 


I  0r3(R/p)^ 
orbital  inclination 


j  arbitrary  integer  (usually  associated  with  v) 

geopotential  (zonal  harmonic)  coefficient  of  degree  { 

)c  arbitrary  positive  Integer  (usually  associated  with  (O) 

K  -i  J2(R/p)^ 

t  positive  integer  (suffix  for  J),  usually  2  or  3 

li  non-singular  quantity  such  that  L  »  M  qy 

m  V  -  M  (alao;  suffix,  in  r  in  section  1) 

''  M  mean  anomaly 

f  n  mean  motion 

n*  energy-based  fixed-mean  mean  motion 

H  modified  mean  motion  (ni,  in  Ref  4) 
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Le^endre-polyncffiial  function  of  degree  t 

V  (1  -  e*) 

used  as  factor  of  d{/dv  In  analysis 

geocentric  radius  vector,  coordinate  of  (r,  b,  w)  system 

obsolete  cylindrical  coordinate 

equatorial  radius  of  the  Earth  (6378.14  km) 

rotation  matrix  for  j***  axis  (j  -  1,  2,  3) 

derived  rotation  matrix 

sin  1  (of  Cl  and  sj,  Sj,  5|(  by  analogy) 

time  (measured  from  epoch) 

affix  for  matrix  transposition 

argument  of  latitude  (v  (o) 

obsolete  cylindrical  coordinate 

M  +  0) 

quantity  such  that  U'  -  M  +  qffl 

disturbing  function  W2  and  O3  as  particular  Uf) 
true  anomaly 

longltude-llke  coordinate  of  (r,  b,  w)  system 
w  +  ciitp 

q‘3  (p/r)2 

usual  (equator  based)  geocentric  coordinates 
geocentric  coordinates  based  on  the  mean  orbital  plane 


P  geocentric  latitude  (ie  declination) 

Yj  cos(jO  +  SS) 

Fj  cos(j\)  +  i&) 

total  short-period  perturbation  (in  0 
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A 
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A 
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A 

cofactor  (likewise)  in  C 
see  C 
e  sin  (0 
see  C 

arbitrary  angular  quantity 

Earth's  gravitational  constant  (398  600.4  km^/sZ) 
cofactors  of  K  and  K  in  Kepler's  third  law  (modified) 
e  cos  (0 
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